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Forging “new metals” has been routine 

procedure at Wyman-Gordon for 75 years. Contrasting the 
modern marvels of metallurgical development is the oldest 
metal employed by man—copper. Here you see the largest 
copper closed die forging ever produced—a Re-entry 

Shield weighing 1875 pounds (67-1/2” x 21-1/2”). Unexcelled 
know-how, with the availability of the most modern forging 


equipment, assures the ultimate in forging quality to meet 
the constantly expanding demands of progress— man’s quest 
for greater speeds and power in his unending exploration 
farther and farther into yesterday’s unknown. 
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To send the U.S. Pioneer more than 60,000 miles 
into interplanetary space, Space Technology Laboratories 
in seven months designed, developed, assembled, and 
tested an 88-foot combination of three integrated stages witha 
payload incorporating 36 separate ignition systems. 
stt’s Astrovehicles Laboratory focused on the payload itself and 
the sensitively related problems of propulsion, weight, 
and stability. These are in addition to the overall complexities 
of the structural configuration. 


Pioneer, setting new apogees in science and missilery, 

typifies the achievements STL is making in the advancement of 
space technology. Those who are able to contribute 

to and benefit from these developments are invited to 


consider joining our staff. 
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to the field of astronautics, by the current or 
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con by the importance of distributing the in- 
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Grove Powreactor Dome Regulator Model GH-408... 
50-6000 PSI Inlet...5-3000 PSI Reduced Pressure. 
Grove offers a wide variety of regulators for 

various pressure and flow conditions. 


The Control of High Pressure 


luids there’s no margin for error... 
his is no child’s play. Response must be instantaneous...accuracy absolute...dependability 
ositive. This is why for over a quarter of a century, Grove Regulators have been the No. 1 
hoice for the critical control of high pressure fluids. Here at Grove—unmatched experience 


d unique design/engineering skills that can work for you. See or call your Grove man today. 
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Phosphor bronze reticle (actual size) and space frequency transfer characteristics of circular aperture reticle. 


TARGET DISCRIMINATION 
IN INFRARED DETECTION SYSTEMS 


The pioneering field of infrared detection offers 
many challenging opportunities to scientists and 
engineers at Ramo-Wooldridge for advanced 
studies in the solution of target discrimination 
problems. Research is continually under way at 
Ramo-Wooldridge in the integrating of infrared 
detection devices with the latest electronic sys- 
tems techniques for enhanced target detection 
on the ground and in the air. 

The phosphor bronze reticle, or image chop- 
per, illustrated above was developed by Ramo- 
Wooldridge. It indicates a marked stride in space 
filtering discrimination concepts, and is used for 
target signal enhancement in guided missiles, 
anti-aircraft fire control and air collision warn- 
ing applications. 

The reticle is used in the focal plane of an 
infrared optical system and is rotated to chop 
the target image for the desired space filtering. 
It is also employed in time filtering, such as pulse 
length discrimination, or pulse bandwidth fil- 
tering. 

Space filtering is critical to infrared systems, 
because of its ability to improve the detection of 


objects located in the midst of background inter- 
ference. In a manner similar to that used in the 
modification of electronic waveforms by elec- 
trical filtering, space filtering enhances the two- 
dimensional space characteristics of a target. The 
size and features of the target are highlighted 
and the undesired background eliminated. 

Scientists and engineers with backgrounds in 
infrared systems—or any of the other important 
areas of research and development listed below 
—are invited to inquire about current opportuni- 
ties at Ramo-Wooldridge. 


Electronic reconnaissance and 
countermeasures systems 


Analog and digital computers 
Air navigation and traffic control 
Antisubmarine warfare 

Basic research 

Electronic language translation 
Information processing systems 


Advanced radio and wireline 
communications 


Missile electronics systems 


RAMO -WOOLDRIDGE 


P.O. BOX 90534 AIRPORT STATION * LOS ANGELES 45, CALIFORNIA 


a division of Thompson Ramo Wooldridge Inc. 
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Recent Advances in Gaseous Detonation 


R. A. GROSS 
Fairchild Engine Division 
Deer Park, N. Y. 


A. K. OPPENHEIM 
University of California 
Berkeley, Calif. 


A review of recent work in gaseous detonation is presented. Early work is briefly mentioned and 
treatises listed. Theoretical calculations of Chapman-Jouguet detonations are reviewed, com- 
pared and the ambiguity concerning the speed of sound in a reacting gas mixture discussed. Experi- 
mental Chapman-Jouguet measurements are reviewed. Recent studies of the interior of a deto- 
nation wave are presented. Standing detonation wave research, detonation limits, two-dimensional 
detonations, spectra, ionization and magnetohydrodynamic treatments are brought to the reader’s 
attention. A qualitative description of the development of a flame to a detonation is presented. 
Experimental observations are examined and recent theoretical attempts to explain these observa- 


tions are reviewed. 


Dr. Gross, Chief Research Engineer at the Fairchild 
Engine Division, is currently doing fundamental research 
in high temperature gas dynamics, supersonic combustion 
and plasma physics. He received his B.S. degree in 
Mechanical Engineering from the University of Pennsyl- 
vania (1949) and his Ph.D. in Engineering Sciences and 
Applied Physics from Harvard (1952), where he was a 
teaching fellow. Since 1954 he has been working in ad- 
vanced propulsion engineering at Fairchild, and during 
1959-1960 will be on leave of absence as a National Science 
Foundation Senior Post Doctoral Fellow. Dr. Gross has 
written numerous papers on advanced propulsion and 
problems associated with detonation and is a fellow of 
the American Rocket Society. 


Dr. Oppenheim is now Professor of Aeronautical Sciences 
at the University of California, Berkeley, where he has 
taught since 1950; and is also consultant for the Shell 
Development Company. He received the Dipl. Ing. in 
Aeronautical Engineering from the Warsaw Institute of 
Technology (1943), the Ph.D. from the University of 
London and the D.I.C. from Imperial College (1945). He 
has served on the faculties of City and Guilds College, 
London, 1945-1948, and Stanford University, 1948-1950. 
An active member of the American Rocket Society, Dr. 
Oppenheim is a past president of the Northern California 
Section and has recently been elected to the National 
Board of Directors. 


Steady-State Phenomena 


DETONATION is a wave in which an exothermal 
reaction takes place and which moves with supersonic 
velocity with respect to the undetonated (reactant) gas. It 
is characterized by this supersonic propagation velocity and 


Received Jan. 16, 1959. 
_ 1 This review was prepared with the support of the Air Force 
Office of Scientific Research, ARDC, Washington 25, D. C., 
under Contracts AF 49(638)-15 and AF 496(38)-166. The sec- 
tion Steady State Phenomena was written by R. A. Gross, and 
— Development of Detonation was prepared by A. K. Oppen- 
eim. 
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a large pressure and temperature increase across the wave. 
A detonation differs from a flame (deflagration) in that a 
flame moves with subsonic velocity. The propagation rate 
of a detonation is determined mainly by aerothermodynamic 
considerations, whereas flame speeds appear to be determined 
by a complex interaction between high temperature chemical 
kinetics, heat transfer and diffusion. 


History 


The phenomenon of detonation was discovered concurrently 
by Berthelot (1)? and Mallard and LeChatelier (2) in 1881. 
Chapman (3) and Jouguet (4,5) independently postulated 
that the steady-state detonation commonly observed in ex- 
periments is distinguished by the fact that the fluid flow 
immediately behind the detonation wave has a Mach number 
exactly equal to one relative to the wave. Such a deto- 
nation wave has subsequently been called a ‘“Chapman- 
Jouguet”’ detonation or simply “C-J”’ detonation. Follow- 
ing these pioneering works, the literature has carried a steady 
stream of research reports by chemists, physicists, mathe- 
maticians and engineers who became interested in this phe- 
nomenon. Reviews of this work can be obtained by referring 
to books or treatises that present comprehensive information 
on detonation (6 to 13). Work in this field transcends na- 
tional boundaries; important contributions having been made 
by French, English, German, U.S. and Russian scientists. 


Classification of Flows 


Combustion is distinguished by the numerous and often 
complex types of phenomena that occur in nature. To ob- 
tain a qualitative understanding of some, examine the solu- 
tions of the one-dimensional steady flow equations of motion 
where thermal energy is added to the flow. In Fig. 1 are 
shown the solutions of the equations representing conserva- 
tion of mass, momentum and energy. is the velocity (rela- 
tive to the wave) ahead of the detonation, V2 the product 
velocity and M, the approach Mach number. Details con- 
cerning these equations and their solutions can be found in 
(11, 13, 14). Important points in this figure are labeled, and 
their physical implications noted. Points 6, 7 and 10 are 


2? Numbers in parentheses indicate References at end of paper. 
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POINT | - DEFLAGRATION (SIMPLE SUBSONIC COMBUSTION) 
2 - EXPANSION SHOCK + COMBUSTION 

3 - EXPANSION SHOCK 

4 - COMPRESSION SHOCK + HEAT WITHDRAWN 
5 - COMPRESSION SHOCK 

6 - COMPRESSION SHOCK + COMBUSTION 

7 - COMBUSTION; NO SHOCK 

8 - NOTHING HAPPENS 

9 - HEAT WITHDRAWN, NO SHOCK 

10 - CHAPMAN - JOUGUET POINT 


Me 2 


Fig. 1 One-dimensional steady flow with heat addition 


detonations. Point 10 is the C-J point and represents the 
singular solution of the equations. Point 6 is a “strong” 
detonation and is characterized by the fact that the product 
flow is subsonic. Point 7 is a “weak” detonation and has 
the property of a supersonic product flow. A strong deto- 
nation has a larger pressure ratio across the wave than a weak 
detonation at the same value of M;. The C-J point is the 
coalescing of the strong and weak detonation branches and 
has a product Mach number of 1, or “choked” flow. There 
is a corresponding subsonic C-J point, although it does not 
usually bear this name. Between the two C-J points, for a 
given amount of heat addition, there are no steady flow solu- 
tions to the equations of motion. For a given quantity of 
heat addition, the C-J points represent a minimum super- 
sonic wave propagation rate and a maximum subsonic wave 
propagation rate, respectively. 


Detonation Wave Propagation Calculations 


The fundamental physical problem in combustion is the 
determination of the combustion wave propagation rate or 
speed and the attendant change in the properties of the 
gases. Since the early work of Chapman (3) and Jouguet 
(4), scientists have been able to predict C-J detonation ve- 
locities with fair precision. The equations of conservation 
of mass, momentum, energy, chemical equilibrium of the 
product flow and the fact that the product flow Mach num- 
ber is 1 constitute a determinant problem. These equations 
consist of a large simultaneous set of nonlinear algebraic 
equations. Because of the labor involved in obtaining solu- 
tions, the majority of published reports concern relatively 
simple two atom chemical systems, such as hydrogen-oxygen. 

Complete detonation calculations which include dissociation 
of a three or more atom system are given in (4,9,14,15,23 and 
24). All detonation computations have been concerned with 
the same system of equations, although the numerical methods 
of solution have varied. Progress has been made primarily 
in the use of improved thermodynamic data (14 to 23) 
and care in the proper use of the appropriate sound speed for 
the C-J point. Recently (14,23,24) electronic digital com- 
puters have been employed to solve more complex systems 
with high accuracy, and we can expect considerable future 
contributions in this field from such devices. 

Recently there has been considerable concern about the 
uniqueness of the definition of the C-J point. This involves 
the question of the velocity of sound in a reacting gaseous 
medium. Surprisingly, nearly al] theoretical calculations of 
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the C-J point contain a theoretical error that arises from the 
ambiguity of the speed of sound definition at the C-J plane. 
This point is amplified in the next section. 

To see how such present computations compare, consider 
the relatively simple case of the stoichiometric hydrogen- 
oxygen C-J detonation data of Table 1. Unfortunately, 
many authors describe the source of their thermodynamic 
data inadequately, and they do not give the chemical com- 
position of the product flow. Such data are presented, how- 
ever, in references such as (14,15 and 17). 


Speed of Sound in a Reacting Gaseous Medium 


The question of the appropriate sound speed in high te1- 
perature gas mixtures was raised by Brinkley and Richardson 
(25) and has been followed by an interesting literature debate 
recently summarized by Wood and Kirkwood (26). Tiie 
problem centers about the difference between an equilibrium 
sound speed and a frozen sound speed. Consider a gus 
whose equation of state is given by 


P = P(p, 8, Bi. . Bn) [t} 


where p is the pressure, p the density, s the entropy and 8; the 
mole fraction of the ith chemical constituent of the mixture. 
Then the equilibrium speed of sound a, is defined 


Op), ar =0 
_ op) (28: 
Op 8, Bi op; 8, By Op s, AF = 0 


i=1 
where 
AF, = Vim Mi 


AF, is the Gibbs tree energy; yu; is the chemical potential of 
substance 7, and »y,* is the stoichiometric coefficient of the 
ith substance in the ath chemical reaction. 

The term (Op/Op)s, g, is the frozen speed of sound, i.e., the 
sound speed in a gas mixture of constant chemical composi- 
tion. This is the sound speed familiar in low speed aerody- 
namics. The summation in Equation [2] corresponds to 
the contribution to the sound speed by the slight change in 
chemical composition associated with the change in pressure 
and temperature in the sound wave itself. For a fixed chemi- 
cal composition, the sound speed of an ideal gas is given by 


re) R 
Op m 


where 7’ is the temperature, y the ratio of specific heats and 
R the gas constant. When the fluid is a mixture of ideal 
gases 


ZBiCpi 
Co “4 


where C>; is the specific heat at constant pressure of the 7th 
species, and m is the molecular weight of the mixture. Most 
calculations of the C-J point, a priori, set the product flow 
Mach number equal to unity using Equation [3]. However, 
the classical C-J point, which is the singular point, has the 
product Mach number equal to unity only if the product 
Mach number is defined in terms of Equation [2]. To illus- 
trate this point consider Fig. 2. Here is shown the classical 
Hugoniot curve which is interesting because the coordinates 
are purely thermodynamic (velocities having been elimi- 
nated). The data shown here for a hydrogen-oxygen deto- 
nation were computed using the computer program of (14) and 


[4] 
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thermodynamic data of (24). Point B is the classical C-J 
point in which the equilibrium product Mach number is 
unity. It is the singular point, i.e., point of tangency on the 
product Hugoniot curve, and has the additional property 
that ds = 0 where the derivative of the entropy s is evaluated 
along the product Hugoniot curve (10). C is the point found 
in all theoretical calculations using Equation [3]. The value 
of V;, the detonation speed, the quantity most often experi- 
mentally measured, is quite similar between B and C, but 
the pressure, temperature and composition are significantly 
different. 

The frozen speed of sound at point B is, for most C-J 
stoichiometric fuel air mixtures (14), larger than the equilib- 
rium speed of sound and of the order of a,/a, ~ 0.96. 


Ex perimental Measurements of C-J Detonations 


shock and detonation tubes have been used to measure the 
veiocity of propagation of C-J detonations (17,27,29). Ex- 
perimental work in this area is difficult because of the very 
fast response instrumentation required. Most hydrocarbon- 
air stoichiometric mixtures under room conditions detonate 
with velocities of the order of 6000 ft per sec or a Mach number 
of about 5 (14). The use of piezoelectric, ionization and re- 
sistance gages have been the usual techniques employed. 
Kistiakowsky (32) and Duff (35,36) have succeeded in making 
dnsity ratio measurements across detonation waves using 
x-ray absorption techniques. The effect of initial pressure 
ald temperature upon detonation speed were experimentally 
examined by Cannon (33) and Hoelzer (34). 

The question concerning the product gas sound speed 
mentioned previously appears difficult to resolve experi- 
mentally. Relaxation effects (37) cause some doubt on the 
question of equilibrium. More important, however, is the 
fact that the flow field behind a detonation gives evidence of 
two-dimensional effects. A Schlieren photograph of such 


2 He + Op DISSOCIATED PRODUCTS 


FROZEN 


Vi #9323 FT/SEC 
V2 = 5247 FT/SEC 
atm. 

Tz =6595°R 

=5247 FT/SEC 
5053 FT/SEC 


C-J (equilibrium) Me =| \ ac 


Vi = 9316 FT/SEC \ 
Ve = 5065 FT/SEC 

p2 = 18.80 atm. \ 

Te = 6616 °R \ 

a, * 5260 FT/SEC 
a, * 5065 FT/SEC \ 


\ 
Pi = | atm. 
Ti = 537°R 


A 
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Fig. 2. Stoichiometric hydrogen-oxygen detonation 
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Table 1 Stoichiometric hydrogen-oxygen detonation 


a = 1.0 atm 
T, = 70F 2H: + O2 — dissociated products 
Measured m———Calculated 
Velocity Velocity Pressure Temper- 
Vi, Pr ature 
Source ft/sec ft/sec atm T2, F 
(14) Bae 9342 18.90 6175 
(i7)° 9268 9357 18.06 6160 
(100)? 9100 8961 17.22 5908 
(15)? 9249 9206 18.05 5989 
Eisen and 
Gross? 8 9316 18.80 6156 
(27) 8800 


* These data should be compared with the frozen speed 
of sound solution shown in Fig. 2. 

> In (14) some of the thermodynamic data used are out- 
dated. Consequently the calculation has been redone 
using the latest thermodynamic information available. 
This accounts for the difference between (14) and the 
more recent Eisen-Gross data given in Table 1. The 
sources of thermodynamic data used in this latest calcula- 
tion can be found in (24). 


effects was recently published by Fay (38) wherein weak 
shocks are seen originating from the wall, and he suggests 
that this is the effect of the viscous boundary layer and heat 
transfer. It has been known for some time that the diameter 
of the experimental tube has an effect upon the detonation 
velocity; elementary analyses have been offered by Kisti- 
akowsky (23, 30) and Zeldovich (39), and partially con- 
firmed by experiment (30). 

Again, it is interesting to compare some experimental results 
for the hydrogen-oxygen detonation. These are shown also 
in Table 1. It can be seen that the experimental results 
cluster near the computed values, but the theoretical values 
are slightly higher. It is suspected that this results from 
the aforementioned two-dimensional effects. 

Finally, some interest has been focused upon the flow behind 
the C-J plane by Taylor (40). A rarefaction wave propagates 
upstream behind the C-J plane so as to establish conditions 
which satisfy the boundary conditions from which the 
detonation originated. Kistiakowsky (41) has made some 
measurements which confirm the Taylor analysis. 


Interior of a Detonation Wave 


Although the detonation speed can with good precision be 
determined solely from aerothermodynamic considerations 
across the wave, considerable interest exists in the wave in- 
terior. In this region the chemical changes take place. Von 
Neumann (42), Doring (43), Zeldovich (39) and others 
analyzed this problem assuming that the detonation wave 
consists of a normal shock wave followed by a deflagration. 
The flame was thought to be initiated by the high temper- 
ature created by the shock wave. Unfortunately, many of 
these early works obscured the real problem by their failure 
to consider gradients within the wave. The concept of a 
shock followed by a deflagration implies the existence of a 
pressure peak immediately behind the shock and prior to 
the thermal acceleration of the flow to a Mach number of 1. 
This maximum pressure is sometimes referred to as the von 
Neumann “spike.” Kistiakowsky (44) and Duff (35) 
looked for this peak, and, although their measurements 
showed qualitative agreement with this theory, they failed 
to find the required maximum spike height. 
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Recently Hirschfelder et al. (45,46) have made a vigorous 
attack upon this problem. They solve numerically the dif- 
ferential equations in the wave interior under some simpli- 
fying chemical kinetic assumptions and obtain composition, 
temperature and pressure as a function of distance. The 
effects of viscosity, diffusion and thermal conductivity were 
included. In the examples Hirschfelder studied, there was a 
strong coupling between the reaction zone and the shock 
zone so that the solutions never exhibited the height of the 
von Neumann spike. 

The physical results of Hirschfelder’s analysis are summa- 
rized in Fig. 3. Here again is shown the Hugoniot curve. 
The loci of the von Neumann solution and Hirschfelder solu- 
tions are indicated. The von Neumann model is like a zero 
order solution to the more general problem treated by Hirsch- 
felder. 


Non-C-J Detonations 


Under certain conditions a detonation propagating along 
a tube will exhibit nonsteady effects ascribed mostly to spin. 
Gordon et al. (47) recently published some experimental data 
on conditions under which such spinning occurs. Fay (48) 
and Chu (49) have developed theories which predict spin fre- 
quency and pitch. These theories are based upon acoustic 
(linearized) treatments of the problem, in which the heat 
release is phased to reinforce the oscillation. The actual 
phenomena are nonlinear, and observed waves are often 
different and frequently complex. 

Steady strong or weak detonation waves have never, until 
recently, been observed either in nature’ or created in the 
laboratory. Although there are many statements in the 
literature claiming that steady non-C-J detonations do not 
exist, suggestions that they might exist are given by Reed 
(50), Hayes (51) and others (52,54). Gross (55) recently re- 
ported on the first laboratory observations of steady, strong 
detonations using hydrogen and air in a specially designed 
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Fig. 3 Interior of a detonation wave 


high temperature, steady flow, supersonic research tunnel, 
Repeated experiments and measurements have demonstrated 
the existence, stability and reproducibility of steady strong 
detonations. This research has opened up a large super- 
sonic flow regime in which steady combustion can be achieved. 


Standing Detonation Research 


Recently, efforts have been made to obtain a detonation 
wave that is stationary in laboratory coordinates. In 1956 
von Behrens (56) reported on the diffusion burning of par- 
tially burned rocket fuel with air in the flow field of a rocket 
exhaust. Nicholls et al. (57,58) have succeeded in burning 
a preheated mixture of hydrogen and air when it flows through 
the Mach reflected shock existing in open space downstream 
of an underexpanded nozzle. It is not clear whether this 
can be called a classical detonation, since light emission from 
apparent combustion took place for large distances down- 
stream of the wave. However, this represented a real stp 
forward in aerothermodynamic research. 

The strong detonations reported by Gross (55) are steady 
standing waves and have been held stationary in the labo- 
ratory for time periods of up to 1 hr; longer times could be ol- 
tained with increased fuel storage capacity. The waves aie 
less than one quarter of an inch thick, one-dimensional, and 
stationary in a Mach 3 premixed stream of hydrogen and air. 
Thermal ignition takes place through a normal shock wave 
when the stagnation temperature is above about 1370 I. 
Once the detonation has been established, unusual proper- 
ties, hitherto unreported, have been observed. Mixtures 
of hydrogen-air have been burned under conditions far 
leaner than previously reported ignition, flammability and 
detonation limits. Of considerable interest theoretically is 
the fact that once the detonation wave is established, the 
upstream stagnation temperature can be reduced greatly 
without apparent effect upon the combustion process in the 
wave. For example, the upstream stagnation temperature 
of the system has been reduced to the order of 500 F and the 
detonation maintained. The degree to which the stagnation 
temperature can be so reduced is a function of the fuel-air 
ratio. These experimental observations give strong support 
to the thesis that transport properties within a detonation 
wave are important (see Hirschfelder (45,46)) since the condi- 
tions behind a shock wave at these stagnation temperatures 
are far too low to produce thermal ignition. This work will 
be reported on in detail in the literature in the near future. 

As this work progresses we may see the first purposeful 
applications of detonation. Hypersonic detonative ramjets 
were first discussed by Roy (59) and recently by others (60 
to 63). Such developments may extend the Mach number 
of air-breathing propulsion systems. The control of deto- 
nation may also open up some interesting chemical production 
methods. 


Miscellaneous 


It is not always possible to generate a detonation. Data 
on detonation limits were reported by Fay (64) and Stein- 
berg (65) using shock tube methods. More recently Belles 
(66) proposed an analytical method for the prediction of the 
detonation limits of hydrogen-oxygen. This problem ap- 
pears closely tied to chemical kinetics, a subject which is still 
evolving. 

Some analytical aspects of steady detonation in two dimen- 
sions have been reported by Manson (67), Zeldovich (68), 
Sedov (69) and others. 

Detonation spectra have recently been studied by Wagner 
(70) and Gaydon (71). Fay (72) has made some shock tube 
ionization measurements behind detonation waves, and 
Chinitz et al. (24) have presented extensive calculations of 
the aerodynamic, chemical and electrical properties behind 
air shocks and H,-air and C,H-air detonations. 
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Fig. 4 Schlieren photographs of the development of detonation 


A study of the electromagnetic effects on exothermal 
waves has been made by Gross et al. (73,101). They study 
the interaction of a transverse magnetic field with a one- 
dimensional detonation. (The latter report contains an 
error in section B and is being revised.) 


Development of Detonation 


The problem of the onset of detonation has baffled the scien- 
tist for a long time. It was, in fact, this problem—brought 
to focus by accidents in coal mines and laboratories—that 
gave the start to the prolific research on the properties of the 
steady-state detonation described in the first part of this 
paper. 

The development of detonation is, by definition, an un- 
steady process. Understanding and exploration of it has 
been handicapped by the lack of instrumentation with suffi- 
ciently fast response characteristics to study the transient 
wave phenomena which control its progress. Such instru- 
mentation is now available. It was developed by improve- 
ments at first of Schlieren techniques combined with micro- 
second light sources and rotating drum or mirror cameras, 
and more recently of electronic apparatus which resulted 
from the progress in shock tube technology. The insight 
gained into the details of the phenomena gave rise to modern 
theories on the mechanism of the processes leading to the 
development of detonation. 


Experiment 


Fig. 4 represents the type of records obtained some 15 
years ago by Schmidt, Steinicke and Neubert (74). It con- 
sists of a series of Schlieren photographs of a flame acceler- 
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ating into detonation in a constant area tube. It was ob- 
tained by means of four Schlieren systems using intermittent 
spark gap light sources and a camera with four rotating drums. 
The record reveals not only the time-space trace (world- 
line) of the flame front, but also world-lines of shockfronts 
which sweep through the combustible mixture ahead of the 
flame. On the basis of Fig. 4 and similar observations, a 
phenomenological description of the process can be con- 
strued as follows. 

The laminar flame front accelerates at first due to the 
thermal expansion of the burned gases. At the same time 
the unburned mixture immediately ahead of the flame is set 
in motion by a series of Mach waves which are generated in 
the course of combustion. The laminar flame front becomes 
first curved, then wrinkled, and finally it breaks into a turbu- 
lent flame brush. In the case of ignition at a closed end, 
the wrinkling is initially due to the fact that the flame is 
pushed from behind by the expanding hot products and is 
consequently exposed to cold gas flow at an increased rate. 
In the meantime the Mach waves coalesce into shock fronts, 
and the flow ahead of the flame induced by their action reaches 
a supercritical velocity and also becomes turbulent, intensi- 
fying still more the process of flame acceleration. 

The breakdown of the laminar flame front into a cellular 
structure, which occurs under the gentle conditions prevailing 
in the vicinity of the flammability limit, is a well-known 
phenomenon. Under the more intense conditions of an ex- 
plosive mixture, this process is so rapid and the wrinkling so 
fast that the combustion zone itself becomes a strong source 
of turbulence and a generator of finite pressure waves. 

Thus, in rapid succession, the reaction zone is subjected 
first to the effects of turbulence and then to effects of elevated 
pressure which increase the density, and consequently the 
volumetric rate of energy release and of high temperature; 
this produces dissociation and, in some cases, even partial 
ionization of the unburned mixture, yielding radicals which 
increase the speed of the combustion reaction. 

The phenomenological description of the development of 
detonation is due to the observation of many investigators. 
Besides Schmidt et al., the more recent contributions worth 
mentioning are those of Turin and Huebler (75), Evans et al. 
(76), Bollinger and Edse (77), Brinkley and Lewis (78), Cook 
(79), Laffitte and Bouchet (80), Martin (81,82) and Kistia- 
kowsky (31). The attention given this problem in the USSR 
is especially significant, as evidenced by the numerous articles 
of Zeldovich (83), Shchelkin (84), Sokolik (85), the coverage 
given it in the recent text of Zeldovich and Kompaneets 
(86), and the very recent experimental studies of Bazhenova 
et al. (87) and Salamandra et al. (88). 

A more sophisticated study than that of flames in pipes 
concerns the problem of spherical explosions. Here especially 
worth noting are the recent experimental studies of Manson 
(89), Freiwald and Ude (90) and Zeldovich, Kogarko and 
Simonov (91). 

The understanding of the process is inadequate yet for the 
establishment of a predictive theory. That is, having been 
given all the initial conditions, we are still unable to derive 
from the theory a specific description of the subsequent de- 
velopment of the process. The state of our knowledge is 
such that the most prolific effort is still spent on the rational- 
ization of experimental results and on the performance of 
more elucidating experiments. 

The study of the development of detonation is concerned 
primarily with wave interaction phenomena—the most im- 
portant feature of unsteady flow processes. 

The most idealized model of a combustion system develop- 
ing into detonation consists of a shock wave which is followed 
by a flame front. In considering such an unsteady double 
discontinuity, Oppenheim (92,93) pointed out that the states 
of the products immediately behind the combustion front 
are not represented by the classical Hugoniot curve, and 
consequently that the classical theory is insufficient for the 
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study of the unsteady phenomena associated with deto- 
nation, such as its stability. For this purpose, he developed 
the theory of the Q-curve, that is, the locus of states attained 
behind a double discontinuity system under conditions similar 
to those of the Hugoniot curve (constant heat release per 
unit mass and given pressure and density of reactants), with 
the additional stipulation that the products of combustion 
move with respect to the flame front at a velocity equal to the 
local velocity of sound. 

The model provides the basis for an idealized description 
of a continuous, dynamically compatible transition from de- 
flagration to detonation. Of particular interest in this con- 
nection is an intermediate state reached by the products when 
they move at the same velocity as the undisturbed reactant 
mixture. At lower temperatures and pressures, the products 
move in the opposite direction, at higher temperatures and 
pressures, in the same direction as the shock-flame system. One 
of the significant results of the analysis is the conclusion that, 
in order for the detonation to be established in a dynamically 
compatible manner as a result of the coalescence of the com- 
bustion front with the shock just before coalescing, both 
fronts must move at a higher velocity than the final steady- 
state value, and the intermediate pressure between the two 
fronts must be considerably higher than the final one of the 
upper Chapman-Jouguet state. 

The algebra of the plane double discontinuity system has 
been explored recently by Troshin (94), while Adams and 
Pack (95) presented a discussion of the dynamic compatibility 
which governs the coexistence of shocks, flame fronts and 
rarefaction waves in a nonsteady combustion system which 
develop into detonation. Popov (96) used the method of 
characteristics to rationalize the experimental records of the 
development of detonation in closed tubes, and, in particular, 
to demonstrate that the observed fluctuations of the flame 
front are due to reflections of waves from the closed ends. 

Especially noteworthy is the thorough algebraic and arith- 
metic analysis of the double discontinuity system presented 
by Taylor and Tankin (97). In their exposition, they re- 
strict their attention to only two points of the Q-curve: The 
state corresponding to the case in which the velocity of the 
products is the same as that of the undisturbed reactant mix- 
ture, and the upper Chapman-Jouguet state which—they 
confirm—has to be exceeded in the course of the development 
of the process as a necessary condition for the flame to catch 
up with the shock wave. They also analyze the case of dy- 
namic compatibility of a spherical deflagration following a 
spherical shock wave and obtain some interesting specific 
solutions for such a double discontinuity system. 

In a recent report based on lectures given by von K4rmén 
at the Sorbonne, Millan (98) describes the development of the 
process following the account of Zeldovich (83). His ap- 
proach to the problem is essentially descriptive, but special 
emphasis is given to the action of the shock front preceding 
the accelerating combustion wave. A more detailed discus- 
sion of the transition from normal combustion to detonation is 
in section 18 of the text of Zeldovich and Kompaneets (86). 

The wave interaction phenomena occurring during the de- 
velopment of detonation which have been observed by 
Schmidt et al. (Fig. 4) were analyzed recently by Oppenheim 
and Stern (99) who rationalize the results and demonstrate 
how they could be exploited for the study of reaction kinetics 
of accelerating flames. 
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Performance Calculations 
for Hybrid Nuclear-Chemical 
Rocket Propulsion Systems 


LEON GREEN Jr.' 
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Preliminary performance calculations for propulsion systems in which water and light metal 
propellant components are vaporized by nuclear heat prior to chemical reaction indicate that the 
hybrid system can remove limitations due to high heats of vaporization on the capacity of light 
metals to generate combustion products with high available enthalpy, and avoid the limitations 
placed on gas temperatures achievable by pure heat transfer systems by material strength require- 
ments. The hybrid system can thus reconcile high specific impulse performance with high pro- 
pellant density and avoid the use of cryogenic propellants. A semiquantitative comparison of the 
relative cost and feasibility of achieving a given flight performance with various propellant systems 
indicates that the hybrid water-metal systems compare favorably with high energy chemical and 
pure heat transfer systems. A qualitative discussion of the engineering difficulties faced by the 
various nuclear systems indicates that, by virtue of avoiding the material protection problem, de- 
velopment of a low-cost hybrid water-magnesium system requires no major technological advances 
beyond today’s achievement, although systems employing lithium promise higher performance. 
It is concluded that the indicated performance achievable with hybrid systems is attractive enough 
soe to warrant more detailed consideration of such systems for possible application to ambitious 

missions. 


HE CONSTANT demand for increased rocket perform- 

ance has prompted several evaluations of the ultimate 
specific impulse which can be derived from chemical propul- 
sion systems, notably those of Carter (1),3 Zwicky (2) and 
Tormey (3). These studies generally agree that the upper 
limit of specific impulse values achievable by the ambient 
temperature reaction of molecular chemical species will be 
about 400 Ibf sec/Ibm,‘ if not lower. To attain impulses 
beyond this value, use of nuclear heating, free radicals or 
still more exotic methods would be required, and the limita- 
tions of all such schemes have not yet been fully assessed. 
The purpose of the present paper is to compare briefly some 
limitations of the two methods of thrust generation which 
are practical in today’s technology, namely, chemical and 
nuclear heating of a working fluid, and to suggest that hybrid 


promoted by the presence of monatomic products in the ex- 
haust. Tormey suggested the use of the light metals lithium, 
beryllium and boron as propellant components because of 
their low atomic weight, and Zwicky recommended their use 
(together with that of magnesium, aluminum and silicon) 
because of their high heat of reaction with oxygen or fluorine. 
Carter, however, noted that these extremely high heat re- 
leases per unit weight cannot be fully exploited in conven- 
tional rocket systems, because of the high heat of phase change 
required for conversion into the vapor state. This fact is 
illustrated by Table 1, which presents pertinent physical 
property data for lithium, magnesium and aluminum for 
comparison with those of more conventional propellant 
components. 


nuclear-chemical propulsion systems may deserve more de- t 
tailed consideration than they have hitherto received. I 
Table 1 Pertinent physical properties of selected I 
Performance Limitations of Chemical Systems propellant components ‘ 
It is recognized that to achieve high specific impulse a ; ae Approx V 
chemical propellant system should have a high energy con- _ Melting Boiling Le, dm, Lr cost, r 
tent (to yield a high stagnation enthalpy), a low molecular Fluid point,C point,C cal/g g/cc cal/ee $/Ib 
weight and a high effective value of specific heat ratio y = H, —259 —253 108 0.07 7.6 1.5% n 
C»/¢» for the nozzle expansion process.6 High y-values are F, —-2233 -187 41.1 1.14 4.8 4.0 n 
Presented at the ARS Spring Meeting, Dallas, Texas, March NH; —78 —33 327 0.69 226 0.15 6 
1958. H.0 0 100 539 1.00 517 0 
1 Head, Applied Mechanics, Corporate Research Division. N-H 1.4 113 ; Pp 
Member ARS. oH, 335 1.01 337 25 
2 Corporate Research Advisor. Member ARS. Li 179 1317 4680 0.51 2380 8.0° a 
3’ Numbers in parentheses indicate References at end of paper. Mg 650 1103 1137 1.57 1790 = 0.36 
‘ All specific impulse epee in this paper are referred to a Al 660 2450 3050 2.2 6700 0.26 % 
pressure expansion ratio of 34 to 1. a 
* ‘Yn the basis of an interpretation of the idealized specific im- * Hypothetical value based a arbitrary assumption b 
ptuse equation in which the y-value for the combustion condition that present cost of about $15/lb in laboratory scale pro- ~ 
(introduced by the perfect gas approximation for enthalpy in duction by NBS can be reduced by about a factor of 10 in 
terms of temperature) was not distinguished from the effective quantity production. P 
y-value governing the expansion process, Tormey states impre- > In large quantity production. eI 
cisely that a low y-value is desirable. of 
ARS JournaL N 


180 


| 
| 
: 
] 
\ 
€ 


It is apparent that if the latent heat of vaporization were 
supplied by some external source, the high heat of reaction 
of the light metals could be manifested in the production of 
superheated combustion products with high values of en- 
thalpy available for expansion work, and thus exploited for 
propulsion. This observation immediately suggests the use 
of a nuclear reactor for purpose of propellant vaporization 
and superheating to a vapor pressure level sufficient for 
thrust chamber operation. In this fashion, one energy con- 
tent limitation on light metal chemical propellants could be 
removed. 


Performance Limitations of Nuclear Systems 


The use of a nuclear reactor to heat a working fluid for 
rocket propulsion purposes has been widely discussed, and a 
description of a conceptual design for such a heat transfer 
rocket, together with a good bibliography of the open liter- 
ature on nuclear rocket propulsion, has recently been pro- 
vided by Rosenblum, Rinehart and Thompson (4). In the 
presently practical case wherein heat is generated by a fission 
process within solid fuel elements, the principal factor limiting 
achievable performance is the temperature which can be 
tolerated by the fuel element without structural failure and 
at which chemical attack by the propellant is not too severe. 
A closely allied limitation results from the problem of heat 
removal from the fuel element by a gas of low molecular 
weight without suffering either too great a difference between 
the maximum solid temperature and the bulk gas temper- 
ature, incurring a performance sacrifice, or too great a pres- 
sure drop in the core, incurring a penalty in the form of in- 
creased pumping power necessary and increased turbo- 


_ machinery weight. 


The design outlined by Rosenblum et al. is the often 
mentioned type employing fuel elements fabricated from 
graphite, a moderator material of exceptional strength vs. 
temperature behavior in that its short-time breaking strength 
reaches a maximum value at a temperature of about 2500 C 
(5), and which even at temperatures of incandescence is not 
susceptible to erosion by high velocity flows of chemically 
inert fluids free of particulate matter (6). At temperatures 
of this order, however, graphite is subject to deformation by 
creep, and no experimental information is available concern- 
ing the degree to which the rate of the thermally-activated, 
diffusional processes which are manifested by creep may be 
increased by a high neutron flux® or under the bombardment 
by relatively massive and highly charged fission fragments 
which would be experienced by fuel-loaded material. In any 
event, it is doubtful that a graphite fuel-element structure 
can withstand temperatures of more than about 2400 C for 
even a relatively short time, without suffering creep deforma- 
tion in excess of acceptable design limits. When the tem- 
perature drop in the solid and gas film is considered (as will 
be discussed later), it appears that a bulk gas temperature of 
about 4500 R (about 2230 C) represents the most optimistic 
value which can be envisaged in conjunction with a graphite 
reactor. 

Thus the concept of a reactor used as a vaporizer of a light 
metal chemical propellant again suggests itself. By such 
means, the temperature limitation imposed by the reactor 
core material could be avoided, the core operated at a more 
conservative temperature at which an increased variety of 
possible cladding materials is usable, and combustion flame 
temperatures far in excess of the stagnation temperatures 
achievable by a solid-to-gas heat transfer system could be 
realized. In contrast to the case where heat generated within 
a solid must be removed at high rates and the material must 
be operated at its maximum usable temperature, high tem- 

6 An order-of-magnitude estimate by G. Schoek (J. Appl. 
Phys., vol. 29, 1958, p. 112, indicates that the creep rate of 


graphite at temperature 7’,/2 would be increased by a factor 
of 2 at a flux of about 10! fast neutrons/cm? sec. 


Marcu 1959 


ATOMIC PERCENTAGE 


20 40 460 70 80 90 
T 

4 1200 
600 

4 1000 
500 4 

4 800 
400 
300 +600 

200 400 
100 


Mg 10 20 30 40 50 60 70 80 90 Li 
WEIGHT PERCENTAGE 


Fig. 1 Liquidus curve for lithium-magnesium (from ‘‘Liquid 
Metals Handbook’’) 


peratures generated by gas phase reactions can be contained 
in thrust chambers of conventional materials by employing 
efficient wall cooling methods, such as porous wall tran- 
spiration. 


Chemical Performance of Vaporized Metal-Water 
Systems 


The foregoing discussion has suggested the possibility that 
the principal performance limitation of high energy, light 
metal chemical propellant systems (i.e., low available en- 
thalpy of combustion products due to high heat of vapori- 
zation for the metals) and of nuclear heating systems (i.e., 
attainable gas temperatures limited by material strength 
considerations) could be obviated in a hybrid nuclear- 
chemical system in which the propellant components are 
vaporized by nuclear heat and then reacted chemically to 
yield high enthalpy combustion products. In order to in- 
vestigate the potential performance achievable by such a 
scheme, sample performance calculations for the systems 
lithium-water (with an excess of lithium as working fluid) and 
magnesium-water (with an excess of water as working fluid) 
were made. In addition, a calculation for a mixed lithium- 
magnesium fuel system was included, since the two metals 
form a continuous solid solution of variable melting point 
over a wide range of composition, as shown in Fig. 1. In 
this system, magnesium and water were in stoichiometric 
proportions for reaction, and lithium served as the principal 
working fluid. Water was arbitrarily selected as an oxidizer 
in preference to fluorine or oxygen, because of its obvious 
cost and availability advantages, because it provides some free 
hydrogen in the exhaust because of its noncryogenic nature, 
and because its dissociation permits the existence of high 
enthalpy combustion products without excessively high flame 
temperatures. If these considerations are not governing, the 
substitution of fluorine would yield higher flame temperatures 
and possibly higher specific impulse values. 

As a comparative basis for the calculations, it was assumed 
that both the metal and water were heated to a temperature 
of 1500 C (an arbitrary figure, 1000 C below the maximum 
strength point of graphite) before reaction. (Actually it 
would not be desirable to superheat water any higher than is 
necessary to obtain a sufficient pressure level.) The calcu- 
lations also assume that the reaction products are expanded 
isentropically over a pressure ratio of 34 to 1, so that the 
results might be compared with published performance data 
for various chemical systems under expansion from 500 to 
14.7 psia. Assumption of such a relatively low expansion 
ratio actually tends to penalize the water-metal propellant 
system in comparison to some of the other systems, since a 


181 


| 

al 
he 

ht 

ns 

he 
ns 
nd 
| 
1e 
es 
zh 
us 

m, 
of 
ise 
n) 
n- 
ge 
4 
is 
al 
or 

nt 

L 

A 


more complete expansion, such as would prevail at high alti- 
tude conditions, would permit greater recovery of the energy 
in any dissociated water or hydrogen fraction for useful 
nozzle work. 

In this latter regard, Tormey (3) states that only the trans- 
lational portion of the energy in rocket combustion gases is 
available for mechanical work, and that all energy present in 
the form of rotational and vibrational modes, as well as that 
due to dissociation and ionization, is not available. This 
statement is equivalent to saying that in a rocket nozzle, 
expansion time is so short that all products behave as mon- 
atomic gases, with a “frozen’’ specific heat ratio of 1.67 during 
expansion. For sufficiently rapid expansions (very small 
nozzles with very high expansion ratios), this statement may 
be accurate. In actual rocket nozzles for motors of thrust 
levels of 1000 Ibf and greater, it may be inaccurate. Effective 
exhaust velocities may be about 10,000 ft/sec, and nozzle 
lengths will vary from a few inches for 1000-lb thrust to a few 
feet at 100,0001-b thrust. Assuming an average velocity one- 
half the effective exhaust value, nozzle lengths of 1 to 2 ft 
correspond to total residence time in the nozzle of 2 to 4 X 10~4 
sec (7,8,9). Relaxation times for rotational and vibra- 
tional energies are of the order of 10~" sec at high temper- 
atures. Specific reaction rates under rocket nozzle condi- 
tions are such that substantially complete reaction (reassoci- 
ation) occurs in about 10-8 to 10~ sec, and ionization rates 
and the reverse are of the same order (10,11). Under these 
conditions, then, there is ample time for these other forms 
of energy to contribute to the work done by the gases in the 
nozzle. Altman and Carter (11) state that reassociation 
during expansion may effect increases in the exhaust velocity 
as high as 10 per cent for the Oo-H» and F.-He systems, anda 
similar effect would be expected for water-metal systems 
with excess water. 

There is at present considerable discussion concerning the 
participation of these various nontranslational energies in the 
expansion process (12). It should, however, be kept in mind 
that although relaxation times for vibrational energy, for 
instance, of a pure gas of homopolar molecules may be rela- 
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Fig. 2 Performance data for system water-lithium 


tively long, these times may be greatly shortened by the 
presence of other molecules, particularly if they are hetero- 
polar, such as water. The question of contribution of non- 
translational energy to increase of velocity in a nozzle must be 
regarded as debatable, at least insofar as the completeness 
of this energy conversion is concerned. Present evidence 
provided by performance tests with conventional propellants 
indicates the conversion is at least a large fraction of what 
would be predicted on a complete equilibrium basis. If, as 
is proposed for some spaceflight programs, combustion pres- 
sures are dropped to near normal atmosphere, and pressure 
expansion ratios are increased to 100 or more, the question 
would become much more serious. This is so because of the 
larger fraction of nontranslational energy, due to vibration an:| 
dissociation, and to the predictable decrease of conversion 0: 
this energy, due to lower density during expansion. 

The calculations also assumed that thermal and kineti: 
equilibrium is maintained during the expansion process. 
Dissociation of metal oxides was neglected. The thermody. 
namic data from (13 and 14) were used where applicable 
The heat content and entropy values of magnesium oxid 
(assumed solid over the entire temperature range) wer 
based on specific heat data given by Kelley (15). The heai 
content and entropy of lithium oxide were calculated by on 
of the authors (JMC). In this calculation, specific heats o: 
the vapor were estimated from those of steam, using differ- 
ences in force constants and in atomic masses to estimat 
vibrational specific heats. The heat of vaporization wa: 
estimated from scanty vapor pressure data, and by compari- 
son with the entropy of vaporization for other substances 
boiling at about the same temperature. In view of uncer- 
tainties, no heat recoverable by condensation of lithium 
oxide was assumed in calculating performance. 

The results of the calculations are presented in Figs. 2, 3 
and 4. From these curves, typical performance points at 
one mixture ratio for each system were selected for com- 
parison with the high energy chemical and nuclear systems 
included in Table 2. In this selection, no attempt was made 
to optimize the water-metal mixture ratio by means other 
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& Table 2 Propellant system comparison (pressure expansion ratio 34 to 1) 
Type > Chemical —Nuclear-Chemical— — Nuclear ~ 
F- H.0O- H.,O- H,0- 

Propellant system —O),-H, F.-H.-— NH; N2H, Li Li-Mg Mg NH; 

mixture ratio, 8.0 3.5 1.9 4.5 3.0 2.0 0.60 0.39 1.9 
(wt. oxidizer/wt. fuel) 

chamber stagnation 6630 4960 8990 5460 7730 8200 7040 6080 6300 4500 4500 
temperature, 7’, (R) 

molecular weight in chamber, 16 9.0 18 8.9 19 19 11.1 22 2135 2 8.5 
M (mol) 

effective specific heat ratio in 1.23 «641.26 1.34 1.33 1.33 1.33 1.21 1.138 1.05 1.4 1.3 
expansion, 7 

bulk propellant p(g/cec) 0.43 0.26 0.75 0.32 1.16 1.30 0.625 0.842 1.15 0.07 0.69 
density, (Ibm /ft*) 2.8 16.7 46.8 20 72.3 81.0 39.0 52.5 71.5 4.4 43.0 

mass specific impulse, 317 364 338 373 305 = 316 423 362 278 722 351 
(Ibf see/Ibm) 

vol. specific (Ibf sec/Ibm X g/ec) 136 94.6 25.3 119 354 = 411 264 305 318 50.5 242 
impulse (10? lbf sec/ft?) 8.50 6.07 15.8 7.46 22.1 25.6 16.5 19.0 19.9 3.18 15.1 

mass flow rate for thrust of 1580 =1370 1480 1340 1640 1580 1180 =61740 692 1420 
500,000 Ibf (lbm/sec) 

vol. flow rate for thrust 58.8 82.2 31.7 67.0 22.6 19.5 30.3 26.3 22.0 157 = 33.1 
of 500,000 lbf (ft*/sec) 

propellant cost rate for above 334 990 5740 4760 4980 5540 5900 = 3920 216 1040 213 
thrust ($/sec) 

reactor power to heat propellants 9:2 7.4 4.3 15 8.2 
(10° btu/sec) 

propellant wt. fraction for 0.91 0.88 0.90 0.87 0.92 0.91 0.84 0.88 0.94 0.66 0.89 
burnout vel. of 25,000 fps 

feasibility parameter(¢ 0.043 0.036 0.082 0.045 0.10 0.13 0.11 0.12 0.075 0.036 0.084 


than visual inspection. The data for pure chemical systems system H.O-Li-Mg shows an intermediate value of 362 Ibf 


were obtained from (3, 16 and 17). 

It may be seen from Table 2 that, in terms of mass specific 
impulse, the H,O-Li system appears highly attractive, its 
value of 423 lbf sec/Ibm surpassed only by the nuclear heated 
Hz system at 722 Ibf sec/Ibm. The system H,O-Mg ranks 
lowest on this scale with a value of 278 Ibf sec/Ibm, and the 
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Fig. 4 Performance data for system water-lithium-magnesium 
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sec/lbm, which surpasses the nuclear NH; system and com- 
pares favorably with the O-H, and F.-H: chemical systems. 
In terms of volume specific impulse, the order shifts consider- 
ably, with the high density F.-N2H, and F.-NH; chemical 
systems heading the list at 25,600 and 22,100 lbf sec/ft’, 
respectively, followed by the hybrid system H,O-Mg at 19,900 
Ibf sec/ft’. On the volume specific impulse basis, all the 
low density systems employing substantial fractions of ele- 
mental hydrogen compare unfavorably. 

In view of the approximations noted above, necessitated by 
the general meagerness of high temperature data for metals 
and oxides, it is not claimed that the performance data pre- 
sented are the results of rigorous calculations. Nevertheless, 
they are believed accurate enough to permit the following 
semiquantitative comparisons of the feasibility of achieving 
a given flight performance. 


Comparison of Chemical, Nuclear and Hybrid 
Systems for a Given Flight Performance 


As illustrated by the trend of the above data, the relative 
attractiveness of a propulsion system cannot be assessed by 
consideration of specific impulse values alone; flight per- 
formance and cost data are required for this purpose. 
For a truly quantitative comparison of flight performances, 
detailed mission analyses and preliminary vehicle design 
studies are necessary, but interesting semiquantitative com- 
parisons can be made by assuming arbitrary levels of thrust 
and burnout velocity. Since it is evident a priori that 
nuclear powered vehicles will compete with chemically 
powered ones only for ambitious missions requiring high 
thrust levels, a thrust of 500,000 lbf and a burnout velocity 
of 25,000 ft/sec were selected for comparison purposes. 


Relative Propellant Cost 


Specification of the thrust level determined the values of 
propellant mass flow rate and approximate cost rate compared 
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in Table 2. For the systems employing nuclear heating, the 
approximate reactor power levels required to achieve the 
specified temperatures (1500 C for the hybrid systems and 
2230 C for the pure nuclear systems) are also determined. 
Together with knowledge of the mass specific impulse, speci- 
fication of the burnout velocity fixes the required propellant 
mass fraction of the loaded vehicle, as determined by the 
simplified expression of Fig. 5. This expression neglects 
drag and gravity effects, but some compensation for the error 
is provided by use of the specific impulse values for a 34 to 1 
pressure expansion ratio, which will be improved in expansion 
under actual high altitude flight conditions. 

Comparison of the flow rate figures does not lead to any 
new observations but does serve to illustrate the relatively 
large volumetric rates required for the hydrogen systems, 
especially the nuclear heat transfer system, which will in turn 
require turbopumps larger than those needed for the higher 
density propellants. Also illustrated is the large variation 
in reactor power requirement among the various nuclear 
systems (lowest for the water-magnesium hybrid, highest for 
the hydrogen heat transfer system), which will be reflected 
as a significant variation in reactor weight and fuel investment. 

The relative propellant cost rate figures indicate that the 
hybrid water-magnesium and pure nuclear ammonia systems 
are the most economical. Although these figures are rela- 
tively firm, there is considerable uncertainty about the costs 
of the systems using hydrogen, which are based upon a crudely 
extrapolated cost value? as indicated in Table 1. Most 
nebulous of all are the costs of the hybrid lithium systems, 
since the lithium cost figure of Table 1 is for the commercial 
material, 92.5 per cent Li? with a thermal neutron capture 
cross section of 0.033 barn, and 7.5 per cent Li® with a ther- 
mal cross section of 945 barns. Any application of lithium 
as a reactor coolant would thus require either the use of at 
least partially purified Li? with an increased per-lb cost, or 
of a reactor with a fast or epithermal neutron energy spec- 
trum. It is possible that the use of completely purified Li’ 
may not be desirable, since the presence of Li® in the stream 
(if tolerable from control and fuel investment standpoints) 
can provide both a mechanism for direct energy transfer to 
the fluid and a source of a low molecular weight exhaust 
components by the slow neutron (n, a) reaction 


+ on! iH? + 2He* 1] 


7 An estimated cost of about $0.55 per Ib at the liquefaction 
plant has recently been quoted by the Linde Co., but the extrap- 
olated value of $1.50 per lb used here is probably more represen- 
tative of the cost at the firing site. 
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The magnitude of these possibly beneficial effects has not been 
estimated, but will probably not be great. 


Relative Engineering Feasibility 


Despite the considerable difference between the derived 
values of Table 2, it is not possible to evaluate the relative 
attractiveness or feasibility of the various propellant systems 
for the selected mission without an additional criterion. 
The question which needs to be answered is ‘‘How feasible is 
it to build a vehicle of the required propellant mass fraction 
using propellants of the density indicated?”’ Obviously, any 
precise answer would require a detailed preliminary design, 
but it has been suggested that a rough figure of merit may he 
derivable from simple, qualitative considerations (18). By 
definition of propellant mass fraction £ = m,/(m, + m,), 
where m, = mass of propellant, and m, = mass of dry 


vehicle 
= mt 1) 


Since we are seeking a basis for gross comparison only, the 
admittedly gross approximation® that the mass of the pro- 
pellant tank constitutes a certain constant proportion of the 
dry vehicle mass seems admissible 


mM, & me [3} 


Introduction of a second rough approximation that, for a 
given mass of propellant, the mass of the propellant tank 
varies inversely as the mass density of the propellant 


m:<1/p [4} 
leads to the result that 


1 
= constant [5} 


It should again be emphasized that this result permits only 
gross comparisons in view of the crudeness of the approxima- 
tions of Equations [3 and 4], but even these comparisons 
appear significant. 

In order to use Equation [5], it is necessary to evaluate the 
constant involved. Since only a rough rule is desired, it 
appears legitimate to use the statement of Sutton (16) that 
“Tt requires extremely careful design to exceed a mass ratio 
of 0.80; mass ratios approaching 0.90 appear to be the prob- 
able practical limit.” Although this statement might not be 
rigorously exact, it provides a useful basis for comparison. 
It is presumed that Sutton’s judgment applies to single-stage 
vehicles using propellants with a specific gravity of the order 
of unity, in which case a maximum propellant fraction of 
0.90 would correspond to a value of 0.11 for the constant 
in Equation [5]. Values higher than this (i.e., low pro- 
pellant fraction coupled with high propellant density) may 
be considered easily feasible, whereas the achievement of 
lower values (i.e., high propellant fractions using low density 
propellants) would seem relatively formidable tasks. If, 
from a still more optimistic viewpoint, a value of ¢ = 0.95 
is considered the upper limit using propellants of unit spe- 
cific gravity, the corresponding value of the feasibility param- 
eter would be 0.052. In view of the highly approximate 
nature of this feasibility criterion, however, comparisons to 
more than one significant figure are not warranted. 

On the basis of this criterion, it is possible to estimate the 
relative feasibility of achieving the desired performance 
(25,000 ft/sec burnout velocity) using the various propellant 
systems listed in Table 2. It may be observed that the sys- 
tems fall into three categories. Those carrying a substantial 
fraction of their fuel in the form of elemental hydrogen (both 
O»-He systems, the low mixture ratio F:-H, system and the 
nuclear Hs system) show a feasibility parameter of the order 
of 0.04, a value which may be considered too low to warrant 


8 Equation [8] disregards the fact that thrust transmission 
structures do not scale directly with tank weight. 
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serious consideration. The chemical systems F.-NH; and 
F.-N2Hy, together with the hybrid systems H,O-Li and 
H.O-Li-Mg, yield a relatively attractive value of the order 
of 0.1, with a spread too small to permit a preference without 
more detailed consideration. The high mixture ratio F.-H» 
chemical system, the hybrid H,O-Mg system and the nuclear 
NH; system show an intermediate value of the order of 0.08, 
which may be interpreted as indicating marginal feasibility 
requiring evaluation by more detailed study. 

It must be emphasized that these gross observations apply 
only to the mission assumed. Ambitious spaceflight mis- 
sions would demand increased specific impulse performance, 
and the relative attractiveness of the nuclear Hy, system 
would be improved. 


Engineering Problems of Pure and Hybrid Nuclear 
Systems 


Although the present paper is primarily concerned with 
the relative performance of the various propellant systems 
considered, it may be briefly noted that the water-metal 
hybrid systems can avoid some of the engineering problems 
encountered in the pure heat transfer systems. 


Materials 


It was emphasized earlier that the material strength prob- 
lem which limits the gas temperature levels achievable by a 
pure heat transfer system using a graphite reactor may be 
obviated in the hybrid scheme. The problem of chemical 
attack by the propellant, and its prevention, was only noted 
in passing. It is well-known that ammonia seriously attacks 
graphite at temperatures of 1500 C and above, whereas hy- 
drogen is less reactive. The resistance of graphite to reduc- 
tion by hydrogen can be improved by refractory carbide 
coatings (19), but protection against ammonia attack is more 
difficult. 

In the hybrid systems, the use of lithium also involves a 
protection problem, but there are indications that the refrac- 
tory metals molybdenum, columbium and tantalum show 
good resistance up to 1000 C (20), and their use at still 
higher temperatures may be possible. No data on the re- 
sistance of tungsten, rhenium or high melting intermetallic 
compounds to lithium attack is known to the authors. In 
the use of magnesium with graphite there does not appear to 
be any serious corrosion problem (20). Water may be handled 
by limiting the temperature levels to relatively low values 
sufficient to give the required steam pressure, i.e., by using 
water as the rgflector coolant and the liquid metal for the 
high temperattire core coolant. 

It is worth noting that, in the event that the solid-temper- 
ature limit could be raised by use of some material other than 
graphite, the relative standing of the pure and hybrid systems 
here considered would not be affected. It appears that the 
only system which might show a relative improvement would 
be a mixed-metal hybrid system, which could then substitute 
aluminum (with an increased heat of combustion and increased 
mass density) for magnesium. Aluminum was not considered 
here because of the estimate (on a highly extrapolated basis) 
that temperatures of the order of 3000 C would be required 
to yield a usable high vapor pressure. As shown in Fig. 6, 
aluminum is also miscible in lithium. 


Heat Removal 


The data of Table 2 indicate that operation of nuclear sys- 
tems requires high rates of heat removal from the core, al- 
though this demand is reduced somewhat in the hybrid scheme. 
It is evident that the removal of heat at total rates of the or- 
der of 10’ Btu/sec (the order of 104 megawatts) requires 
either a highly efficient heat transfer mechanism, permitting 
high surface heat fluxes, or a highly extended heat transfer 
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Metals Handbook’’) 


surface for utilization of lower flux levels. In the case of 
purely gas-cooled reactors there is no choice but to use the 
extended heat transfer surface, owing to the relatively low 
heat transfer coefficients obtainable with gaseous coolants. 
The types of extended surface mentioned by Rosenblum et 
al. (4) included porous beds and arrays of closely-spaced 
flat plates or small holes, all of which present tedious fabrica- 
tion problems. Since the surface must be protected against 
chemical attack, the magnitude of the protection problem 
also increases with that of the heat transfer surface. 

This is not to say that high heat flux rates cannot be real- 
ized using gaseous coolants. The experiments of Durham, 
Neal and Newman (21) demonstrated that heat flux rates 
of the order of 10° Btu/ft? hr could be obtained in forced con- 
vection to helium flowing through straight, electrically- 
heated, g-in. ID graphite tubes. However, in typical cases 
the differences between the inner surface temperature and the 
bulk gas temperature were of the order of 10? F, corresponding 
to effective film coefficients of the order of 10* Btu/ft? hr F. 
Although better performance might be expected using hy- 
drogen or ammonia owing to dissociation effects, it appears 
that heat removal by gaseous forced convection alone in flow 
through straight channels involves a significant sacrifice in 
terms of high temperature drop between tolerable local values 
in the solid and usable bulk values in the gas, increased com- 
plexity, reduced compactness and increased weight. These 
drawbacks may be overcome by use of porous construction, 
but not without a considerable development effort. 

In contrast to the situation encountered using purely gase- 
ous coolants, the use of water and liquid metals permits 
order-of-magnitude improvements in the heat transfer condi- 
tion. The maximum heat flux rates achievable with boiling 
water (22) and boiling liquid metals (23) are also of the order 
of 106 Btu/ft? hr, but the driving temperature difference 
required for this flux is only about 10 to 20 F. As an ex- 
ample of more efficient heat removal, approximately 90 per 
cent of the power required for heating the propellants in the 
H,0-Li system of Table 2 is expended in liquid phase heating 
and vaporization of the propellants. 


Component Design 


The most important design and handling simplification 
made possible by the hybrid scheme is that cryogenic fluids 
are avoided. The hybrid systems would involve “hot’’ 
vehicles in which the fuel would be maintained in a liquid 
state by electrical heaters within the fuel tank, which could 
be constructed from thin gage metal by conventional fabrica- 
tion techniques. Tank insulation at the temperatures in- 
volved presents no problem beyond the scope of everyday 
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technology, and might indeed be easier than insulation of a 
liquid hydrogen tank. 

Of the various nuclear systems, the only one which would 
appear to involve significant advances in pump development 
is the hydrogen heat transfer system, due to the low density 
and cryogenic nature of the working fluid. The low density 
is not only reflected in a greatly increased volumetric ca- 
pacity requirement relative to the other propellants, as illus- 
trated in Table 2, but also in high pumping head require- 
ments. If for a given flow this head is high enough to re- 
quire more than one centrifugal stage, loss considerations 
may indicate the use of an axial-flow design, thus involving 
increased complexity and weight (24). 

The problems of radiation heating of the propellant and of 
radiation damage to solids (especially to close tolerance turbo- 
pump components at cryogenic temperatures), such as were 
mentioned by Rosenblum et al. (4) in discussion of a hydrogen 
reactor system, would be mitigated in the case of a hybrid 
system, not only by the reduced power requirement, but also 
primarily by the fact that the propellant pump and manifold 
system operates at a temperature high enough to permit 
rapid annealing of any radiation damage, and thus could be 
used to shield more sensitive components. Finally, the 
availability of water as a propellant component presents 
opportunities for various design simplifications employing 
water not only as a thrust chamber and nozzle coolant, but 
also for reflection, shielding and possibly moderating pur- 
poses, and also steam-turbine operation. 


Conclusion 


The foregoing discussion has presented preliminary per- 
formance calculations for chemical propellant systems con- 
sisting of water and light metals, assuming sufficient preheat 
by an external source, such as a nuclear reactor, to vaporize 
the propellants and thus permit the generation of products 
with high available enthalpy in the ensuing combustion. 
The possible advantages over conventional chemical and 
nuclear propulsion systems to be derived by use of the hybrid 
nuclear-chemical system include: 

1 The performance limitation imposed on light metal fuels 
by high latent heats of vaporization is removed by the nuclear 
preheating. 

2 The performance limitation imposed on pure nuclear 
heat transfer systems by mechanical strength properties of 
the core structure is avoided by operation of the core material 
at reduced temperature levels. 

3 By means of this dual removal of limitations, the hybrid 
system can achieve high specific impulse values without using 
low density working fluids. 

4 Problems associated with the use and handling of 
cryogenic fluids are eliminated. 

5 Reactor power requirements are relaxed by the utiliza- 
tion of chemical energy to supplement the nuclear energy. 

6 Use of propellants with improved heat transfer char- 


acteristics and reduced radiation sensitivity can permit 
simplified thermal and mechanical core and engine design. 

A semiquantitative comparison of several chemical, nuclear 
and hybrid systems was presented, using a simplified criterion 
for roughly estimating the relative engineering feasibility of 
employing the various systems to achieve a given flight per- 
formance, and it was suggested that the hybrid systems appear 
potentially attractive enough to warrant more detailed con- 
sideration. The low-cost hybrid system water-magnesium 
appears to involve no development advances beyond today’s 
achievement, but use of the higher performance lithium fue! 
would involve a problem in material protection equivalent to 
those involved in the use of hydrogen or ammonia in a pure 
heat transfer system, as well as increased propellant costs 
comparable to those of fluorine systems. If the corrosion 
problem is tractable, still higher performance fuels, such as a 
slurry of lithium hydride in lithium metal, may be envisaged, 
but consideration of such system elaborations is beyond the 
scope of the present prefatory discussion. 
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Magnetohydrodynamics and 


Aerodynamic Heating 


RUDOLF X. MEYER! 


Space Technology Laboratories, Inc. 
Los Angeles, Calif. 


The basic equations and some fundamental concepts of magnetoaerodynamies are discussed. 
As an application, the problem of the flow near the stagnation point of a body of revolution is re- 


viewed and an exact solution is given. 


Data are presented for the heat transfer coefficient at the 


stagnation point, and for its gradient in the stagnation point region. The magnetic field strength 
required to accomplish an appreciable reduction of aerodynamic heating in hypersonic flight is 
discussed, for the case in which ionization is due to thermal motion. Alternatively, methods 
which involve electrical breakdown of the air are considered. 


T HYPERSONIC speeds, exceeding a flight Mach num- 
ber of approximately 15, the air behind a normal shock 
reaches a temperature which is sufficient to ionize the air to 
an appreciable extent. It is then possible to alter the char- 
acteristics of the flow by means of externally applied mag- 
netie fields (1 to 8).2. Due to the motion of the ionized gas 
(plasma) through the magnetic lines of force, currents are 


induced in the gas, which in turn give rise to a force 7 X B, 

per unit volume, acting on the gas. In this expression, 7 is 


the current density and B the magnetic induction. The 
plasma currents in turn give rise to a magnetic field which 


adds to the field produced by external currents. B represents 
the total field. The theoretical analysis of the ensuing flow 
therefore requires in general a solution of Maxwell’s equations 
for the electromagnetic field, coupled with the equations of 
fluid mechanics containing the additional force term. The 
equation of energy is also modified, in a manner which will be 
discussed. 

Confining the present considerations to the stationary flow 
of a fluid of constant properties, the usual equations of mag- 
netohydrodynamics take the form 


curl H = j [1] 
curl FE = 0 [2] 
div B =0 [3] 


where H and E are the magnetic and electric field strengths, 
respectively. Maxwell’s equations are supplemented by the 
constitutive equation 


B = pH [4] 
and Ohm’s law for moving, isotropic media 
j = o(E B) [5] 
where 
= magnetic permeability 
o = conductivity 
v = velocity 
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Furthermore 


= 0 [6] 


since a fluid of constant density is assumed. Supplementing 
the Navier-Stokes equation by the magnetic force term 


Dv 
= ~eradp t+ + B [7] 


p = density 
p = pressure 
= viscosity 


Equations [1 to 7] suffice to determine, for instance, the 
velocity and the current distribution in the fluid. For the 
purpose of computing the distribution of temperature and 
rates of heat transfer, consideration must be given to the 
equation of energy* 


DT 
— 8 
[8] 
where 
c¢ = specific heat per unit mass 
T = temperature 
k = thermal conductivity 
and 
Ov; Ov; Ov; 
Ox; ” (2 + 1) 


The terms on the left of Equation [8] represent the rate of 
heat removed per unit volume by convection and by heat 
conduction, respectively. The terms on the right are the 
rates of heat added through Joule heating and through 
viscous dissipation. 


Stagnation Point Flow 


As an illustration, the rotationally symmetric stagnation 
point flow is considered. The corresponding two-dimensional 
case was treated by Neuringer and Mcllroy (10). In the 
present analysis, a rotationally symmetric magnetic field is 
assumed, which vanishes at a large distance forward of the 
stagnation point and has a normal component on the surface, 

3 A full discussion of the equation of energy in magnetohydro- 
dynamics, applicable to more general conditions than needed 
here, is due to Goldstein and to Chu (9). 
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which is constant. It follows from simple symmetry argu- 
ments that a solution exists, for which 


jr =j. =0 
By=0 

vg = 0 

E, = E, = E, =0 


where r, z and ¢ indicate the corresponding components in 
cylindrical coordinates. It is seen that the current flows 


in circular loops about the z-axis. The magnetic force 7 X B 


° 
Fig. 1 The functions F(¢), G(¢) and their derivatives, for 
a= B=1 
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Fig. 3 Function F;’(¢;) 


is in the meridional plane, and its direction is such as to im- 
pede the motion of the fluid. 
The trial solution 


v, = rf'(z) 

= —2f(z) 

B, rg'(z) 

B, —29(z) 
p—p =h(z) + rk(z) 
T —T = Xz) + 


where f, g, h, k, land m are functions of the coordinate z alone, 
and where p and T are the pressure and temperature at tlic 
stagnation point, satisfies Equations [1 to 9] identically, 
resulting in a set of differential equations for these functions. 
Full details are given elsewhere (11). 

The following boundary conditions are imposed on the sur- 
face, i.e., at z = 0: 

1 The normal component of the total magnetic field ‘s 
constant, i.e., B,(0) = B, where B is a given constant. 

2 The velocity vanishes, i.e., v,(0) = v,(0) = 0. 

3 The surface temperature is constant, i.e. T(0) = T. 

Ata large distance from the surface, the following conditions 
are imposed: 

4 The magnetic field vanishes. In particular, it is re- 
quired that B,( ©) = 0. 

5 The velocity and pressure distributions approach the 
corresponding distributions for an inviscous, nonmagnetic 
stagnation point flow. In particular, v,(©@) = ar, where the 
constant a is determined by the characteristics of the flow 
outside the stagnation point region. 

6 The temperature is independent of the radius. In 
particular, = [where is a given constant. 

Numerical results were obtained by means of a digital com- 
puter. Introducing the dimensionless quantities 


= [12] 

F(g) = f(z) [13] 

GS) = — Fale) [14] 


typical distributions of velocity and magnetic induction are 
plotted in Fig. 1. These results depend upon two dimension- 
less parameters namely 


a = [15] 

B = pon/p [16] 

In aerodynamic applications, a and @ are invariably very 
small, typically of the order of 10-6. The ratio 


_ 


2a 4pa 
which is a measure for the ratio of the magnetic forces to the 
inertial forces, however, is generally of order 1. Letting 
a and £ tend to zero gives rise to a boundary layer problem 
(11). The velocity and the magnetic field distributions then 
depend only upon the single parameter y. Figs. 2 and 3 rep- 
resent the two velocity components outside this boundary 
layer ((; = ¢, PF: =F). Similarly, Figs. 4 and 5 give the 
velocity distribution in the boundary layer, where {2 = 
B-1/2¢ and F,=B-1/2P, 

As one recognizes from the velocity distributions, an im- 
portant effect of the magnetic field consists in displacing the 
flow in the direction of the positive z-axis. The thickness 
of this magnetic cushion, defined by 


f= lim [z — f(z)/a] [18] 


is plotted in Fig. 7. 
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Defining the Nusselt number 
1/2 
(2) [19} 
K \pa 


where A(r) is the heat transfer coefficient, one obtains 


a’r2 
eT T) 
The coefficients Nu and Nu have been computed for a 
Prandtl number of 0.70 and are given in Fig. 8. With in- 
creasing magnetic field strength, the rate of heat transfer at 
the stagnation point is seen to decrease, a result which is in 
qualitative agreement with the computations by Neuringer 
and Mellroy (10) for the two-dimensional case. In the non- 
magnetic case which corresponds to y = 0, the result agrees 
with data computed by Cohen and Reshotko (12). 

Nu jnereases with increasing field strength, indicating 
an increased slope of the heat transfer rate considered as a 
function of the radius. In applications to aerodynamic 
heating, this is a detrimental effect. The crossover point, 
at which the advantage of a reduced heat transfer at the stag- 
nation point is canceled by an increased heat transfer at some 
distance from this point, depends upon the details of the 
application. 


Nu(r) = Nu + [20] 


Required Magnetic Field Strength 


It is seen that y must be large (of the order of 100), if an 
appreciable reduction in heat transfer at the stagnation point 
is to be achieved. The same does not need to be true, how- 
ever, for configurations of the magnetic field differing from 
the one presently considered. The contribution of the plasma 
currents to the total magnetic field is substantial, as seen 
from the curves in Fig. 6 where the magnetic field outside 
the boundary layer is plotted in terms of G,(¢1) = G(¢). 

The required field strength depends critically upon the 
conductivity of the medium. The conductivity of air in 
thermodynamic equilibrium has been estimated for temper- 
atures up to 24,000 K (Fig. 9). The quotient p/pp in this 
figure is the ratio of the density to the density at normal tem- 
perature and pressure. These data have been computed 
from the electron collision cross sections (13) and the caleu- 
lated equilibrium composition of air (14). Details are given 
in an Appendix to this paper. It must be emphasized that 
in view of the uncertainty in the presently available data 
on the appropriate cross sections, the degree of accuracy 
is small. 

In order to permit a quick orientation, the conductivity of 
air behind a normal shock is plotted in Fig. 10 directly as a 
function ot the flight velocity V and of the altitude H above 
sea level. Fig. 11 is a plot of the magnitude of the magnetic 
induction B which is required to make the ratio of the mag- 
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Fig. 10 Estimated electrical conductivity of air behind a normal 
shock 


netic forces to the inertial forces, i.e., ¢B*1/pu, equal to unity. 
Here, o, p and u are the conductivity, density and velocity 
behind a normal shock in air (taking into account dissociation 
and ionization). / is a characteristic length of the aerody- 
namic body, and was taken as 1m. With these assumptions, a 
definite field strength is associated with any given velocity 
and altitude. The curve labeled by the equation w, = », is 
discussed in the Appendix. 

It appears from these data that large field strengths are 
required if the heat transfer at the stagnation point is to be 
reduced substantially at altitudes for which aerodynamic 
heating is critical. Different configurations of the magnetic 
field may be more advantageous. 

In some cases, the temperature of the electron gas may be 
many times greater than the temperature of the gas, corre- 
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Fig. 11 Required magnitude of magnetic field for oB*l/pu = 1 
and / = 1m, behind a normal shock in air 
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Fig. 12 Model used in glow discharge experiments in a super- 
sonic stream 


sponding to an effective increase of the conductivity. This is 
the case with microwave breakdown and glow discharges, 
for instance. The device shown in Fig. 12 contains a large 
number of electrodes of alternating polarity, which are im- 
bedded in an aerodynamic surface. It has been found pos- 
sible to maintain a glow discharge in a supersonic airstream, 
thereby ionizing the air in the boundary layer. A typical 
current voltage characteristic is given in Fig. 13. 


Appendix 


The electrical conductivity of air in thermodynamic 
equilibrium at temperatures up to about 7000 K is known 
with considerable certainty from shock tube experiments 
reported by Lamb and Lin (13). At these temperatures, the 
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contribution of the positive ions to the total electron collision 
probability is still relatively small, but becomes increasingly 
important at higher temperatures. 

It appears that no detailed experimental data are avail- 
able at present for temperatures in excess of 7000 K. On the 
other hand, many of the contemplated applications of mag- 
netoaerodynamics involve considerably higher temperatures. 
In this paper, estimates of the conductivity up to 24,000 K 
are computed. This upper limit was decided upon because 
it coincides with the limit of Gilmore’s computations of the 
equilibrium composition of air (14). 

The usual practice in computing the conductivity of par- 
tially ionized gases is to add the Maxwell-averaged total 
electron collision cross sections (weighted with the number 
density) of the neutral molecules to a corresponding equiva- 
lent cross section of the ions, the latter being computed on the 
basis of a theory which applies to fully ionized gases. The 
degree of approximation involved in this assumption is not 
known. In the case of air, a further element of uncertainty 
is introduced by the fact that the cross sections of some of 
the neutral species, notably N, are not known with certainty. 

The conductivity o is expressed by 


ne? 


As 
mC. 


Ne = number density of free electrons 

e = electronic charge 

Me = electronic mass 

n; = number density of the species j 

Q; = Maxwell-averaged total electron collision cross section 
of the species 


C; = J a the mean speed of the electrons 
TM, 


k = Boltzmann’s constant 
T = temperature 


Thermodynamic equilibrium with an electron temperature 
equal to the gas temperature is assumed. 

The summation in Equation [A-1] is extended over those 
species, neutrals as well as ions, which contribute appre- 
ciably to the total cross section. An inspection of the 
equilibrium composition of air in the range from 6000 to 
24,000 K and a density ratio p/p (po = standard density) 
varying from 10~* to 10 indicates that the neutrals of major 
importance are No, N and O (14). In view of the approxi- 
mate nature of the computation, the inclusion of O., NO 
and other even less frequent species would not be warranted. 
The cross sections for Nz, N and O above an electron tem- 
perature of 6000 K are believed to be quite similar and nearly 
independent of the energy in the range considered. (At 
much higher temperatures, the effect of the neutrals becomes 
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Fig. 13 Typical current-voltage characteristic for glow discharge 
in supersonic air stream. Copper electrodes, parallel to free- 
stream velocity. Cathode area = 3 cm? 
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negligible.) The computations can be based therefore on a 
common cross section of the neutrals, Q,, which was taken as 
0.80 X 107 cm?. In view of the uncertainty in the cross- 
sectional data of N and O, a more detailed computation would 
not be justified. 

The sum over all cross sections, in Equation [A-1], can 


therefore be expressed as 


"hia [A-2] 
NA, 0 +N + 0)Q, + 

0 


where n, and n; are the number densities of neutrals and 
ions, respectively. 4, is the number density of atoms” 
at standard condition, i.e., n4,0 = 5.388 X 10%cm~*. N;, N, 
O and e~ are the number of N-molecules, N-atoms, O-atoms 
and free electrons, respectively, per “air atom.” The values 
computed by Gilmore have been used for these fractions. 

Within the limits of temperature and density considered 
here, the effect of doubly ionized ions on the conductivity 
amounts to less than 8 per cent, and is neglected. Similarly, 
the effect of O~ is negligible. Consequently, for a neutral 
plasma 

= NM = NA, [A-3] 
Po 

a result which was already utilized in the derivation of Equa- 
tion [A-2]. 

Spitzer and Hirm’s (15) results for the conductivity of a 
fully ionized gas, if written in terms of an equivalent cross 
section Q; for the positive ions, give 


(=) In 


A-4 
4kT YE 


for singly ionized ions, where 


Me kT 1/2 
| 


and yz = 0.582. C, is the rms electron velocity 


1/2 
C, = 
Me 


The results of a computation of o, based on Equations 
[A-1 to A-4] are given in Fig. 9. It is noted that the de- 
pendence upon density is relatively weak, and that its gradient 
reverses within the considered temperature interval. Above 
24,000 K, the effect of the second ionization becomes in- 
creasingly important. 

At the lower end of the temperature interval, the computed 
conductivity agrees with the experimental data of Lamb and 
Lin (13). The computed values fall inside the experimental 
scatter. However, this agreement is to be expected, since 
cross-sectional values were used which are derived from these 
data. 

Gilmore’s results for the thermodynamic variables of state 
of the air behind a normal shock resulting from hypersonic 
velocities have been used to calculate the conductivity for 
different flight velocities V and altitudes H above sea level 
(Fig. 10). For comparison, the velocity of a satellite in a 
circular orbit of 7 X 108 m radius and the escape velocity 
have also been indicated. 

In many applications of magnetoacrodynamics, the param- 
eter oB*l/pu plays an important role since it is a measure for 
the ratio of the magnetic forces to the inertial forces. In 
general, this parameter must be at least of order 1, if the 
magnetic field is to produce an appreciable effect. 

In order to permit a quick estimate of the typical field 
strength required in magnetoaerodynamics (assuming thermo- 
dynamic equilibrium), Fig. 11 has been drawn. In this 
figure, the required field strength is given in the case where 
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the value of this parameter is 1, and for a length 1 = 1m. 
From conservation of mass 


pu = paoV 


where po = p«(H) is the freestream density, and the compu- 
tation of B as a function of H and V is straightforward. 

Equation [A-1] gives o correctly only in the case of weak 
fields. Otherwise, the conductivity is no longer described 
by a scalar quantity. This occurs if the cyclotron frequency 
w, of the electrons is comparable or larger than their collision 
frequency ».. From 


Be 
Oo, =— 
Me 
and 
nee? 
MVe 


it follows from equating w, and »,, that 


Nee 


B’ [A-5] 
where B’ can be considered as a rough estimate of the maxi- 
mum field strength B for which Equation [A-1] is still valid. 
The locus of points in the H-V diagram for which this 
limiting condition exists, is indicated in the figure. To the left 
of this curve, w, > ve. 
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Compression Sensitivity of 


Monopropellants 


G. A. MEAD' 


Air Reduction Company, Inc. 
Murray Hill, N. J. 


A method for evaluating the sensitivity of liquid monopropellants to initiation of decomposition 
in the presence of rapidly compressed gas bubbles has been developed. The mechanical design of 
the equipment and the conditions of experimentation simulate the environment which can ¢ ccur 
in actual propellant systems. During the test, a gas bubble in contact with a liquid monofuel sample 
is compressed rapidly by a gas driven piston. The mass and velocity of the piston are analogous to 
the mass and velocity of a liquid propellant flowing in an actual propellant system against a dead 
end. It is postulated that the significant parameter to be associated with decomposition initiation 
by this type of stimulus is the minimum compressive energy per unit gas bubble volume. A 
theoretical analysis and experimental data for several monofuels are given in support of this hy- 
pothesis. The process of data reduction is described. 


ECAUSE monopropellants are necessarily capable of 

exothermic decomposition, a certain amount of risk is 
associated with their handling and use. To minimize this 
hazard, conditions which might result in inadvertent and un- 
controlled decompositions must be recognized and avoided. 
One possible cause of the unintentional initiation of mono- 
propellant decomposition is the rapid compression of small 
gas bubbles in contact with liquid fuel. Rapid compression 
ean result from mechanical shock to containers of fuel, or 
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from rapid closing of valves in propellant lines containing 
entrained gas bubbles. Such phenomena are thought to be 
responsible for a number of accidental explosions. Therefore, 
a method of estimating the relative sensitivity of various 
monofuels to this kind of initiation, or, of quantitatively 
measuring the severity of conditions necessary to cause 
initiation, should have considerable practical value. 

The work described herein is an extension of an investigation 
begun by the Aerojet-General Corp. (1).2 Similar equipment 
has been employed throughout. This equipment was de- 


2 Numbers in parentheses indicate References at end of paper. 
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signed to simulate conditions which might occur in an actual 
propellant line, and incorporates a device for applying pres- 
sure very rapidly to a gas bubble in contact with a liquid pro- 
pellant sample. In the previous work, results were reported 
in terms of the rate of pressure rise necessary to cause ex- 
plosion without a consideration of other thermodynamic 
parameters. In the present study, an attempt has been made 
to explain the significance of rate of pressure rise, relating it 
to other physical conditions of the test in order to estimate the 
energy input to the sample. 

The relative sensitivity of several monopropellants has 
heen measured, and a quantitative scale of sensitivity rating 
has been developed. It is thought that the method of evalua- 
tion should be useful for estimating the sensitivity of experi- 
mental monofuels and of explosive liquids in general. 


Equipment and Operation 


The tester used in this investigation is identical in critical 
dimensions and similar in construction to that first developed 
»y the Aerojet-General Corp. (1). During the course of this 
nvestigation, a few changes have been made to facilitate 
yperation and to improve documentation of the phenomena 
vecurring. A schematic diagram of the test equipment is 
shown in Fig. 1, and a photograph of the disassembled tester 
is shown in Fig. 2. 

The test sample, consisting of a gas bubble in contact with 
liquid monofuel, is subjected to rapid compression by a gas 
driven piston. Driving gas is contained in a high pressure 
tank and is admitted to the driving cylinder by the rapid 
opening of a solenoid valve. Since the area of the piston in 
the driving cylinder is 20 times the area of the piston in the 
sample chamber, a pneumatic multiplication of 20 results. 
Maximum driving gas pressure is limited by the strength of 
the tester body, and with present equipment is about 1000 
psi, giving a maximum static pressure in the sample cylinder 
of 20,000 psi. Rate of pressurization can be varied by adjust- 
ing pressure in the driving-gas tank, by varying the size of an 
orifice in the gas line between the driving-gas tank and the 
driving cylinder, and by changing the working medium. For 
example, the use of helium as a driving gas instead of nitrogen, 
with other conditions being the same, increases the rate of 
pressure rise in the driving cylinder by a factor of 2.8. 

The test chamber itself is formed by the small end of the 
piston, the walls of the cylinder in which it moves and a burst 
diaphragm held in place by a bolted retainer (Fig. 2). The 
burst diaphragm used with recording instrumentation is a 0.1- 
in. thick heat-treated stainless steel disk which deforms elas- 
tically up to at least 20,000 psi. Since the bursting pressure 
of the steel diaphragm is quite high, there is some danger of 
damage to the tester when an explosion occurs. Conse- 
quently, aluminum diaphragms are usually substituted for 
steel when testing without the use of recording instrumenta- 
tion. 

The total sample chamber volume is about 1.3 ml. Samples 
from about 0.4 to 1.1 ml liquid volume and from 0.2 to 0.9 ml 
bubble volume can be tested. The relative sizes of the two 
phases in the sample chamber are limited by the accuracy of 
the liquid phase volume measurement, since bubble volume is 
estimated by the difference. 

The tester is jacketed to provide heating or cooling of the 
sample, if desired. Piston motion is detected by a slide wire 
potentiometer. The sliding contacts are fixed to a shaft which 
is an extension of the driving piston rod. The linear resistor is 
mounted on the tester body and is stationary with respect to 
the piston. The variable resistance forms one leg of a re- 
sistance bridge. Unbalance of the bridge is amplified elec- 
tronically, and the resulting signal causes deflection of a 
galvanometer. A mirror mounted on the galvanometer fila- 
ment, by deflecting a light beam, allows photographic record- 
ing. The rated frequency response of the galvanometer is 
1300 eps, so that the maximum piston velocities encountered 
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Fig. 1 Rapid compression test equipment 


Fig. 2 Rapid compression tester—disassembled view 


(1) Body, (2), insert, (3) piston-plunger assembly, (4) strain gage 
bonded on diaphragm, (5) retainer, (6) flange, (7) cover and (8) 
linear motion potentiometer. 


were well within the range of the instrument. 

Instantaneous pressure in the test chamber is sensed by a 
strain gage cemented to the burst diaphragm. Again, the 
strain gage forms one leg of a resistance bridge, the unbalance 
of which is recorded by a system similar to that used to record 
piston motion. 


Analysis of Tester Action 


It is generally believed that high temperatures resulting 
from nearly adiabatic compression are the cause of ignition in 
the rapid compression tester. The original reason for measur- 
ing both pressure rise and piston displacement in the test 
chamber was that simultaneous readings of these quantities 
would permit the evaluation of n in the equation 


Py Vy 


P = pressure 
V = volume 
T = temperature 
n = polytropic exponent 
0 = initial conditions 
f = final conditions 


Thus, if the initial gas conditions were known and the pres- 
sure and volume of the system could be measured at the time 
of explosion, a temperature at the instant of explosion could 
be calculated even though the compression were not adiabatic 
(n ¥ y, where y = specific heat ratio). This analysis has 
disregarded ignition delay phenomena, which, as will be 
shown later, can have a substantial effect. Further considera- 
tion of the compression process has shown that extremely ac- 
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Fig. 3 Typical oscillograph record of blank test on water 


curate measurement of piston displacement would be needed 
to make this calculation, because the volume of compressed 
gas in the test chamber becomes extremely small at relatively 
low pressures, and large pressure changes are produced by 
very small changes in piston displacement. At higher pres- 
sures, in the range where explosions usually occur, the dif- 
ference between piston displacements for adiabatic and iso- 
thermal conditions is practically negligible on the scale used 
to measure total piston travel. 

Originally, it was anticipated that smooth pressure traces 
would be produced, on which a point of initiation could be 
easily located. A series of such points could be obtained, and 
the initiation pressure extrapolated to infinite rate of pressure 
rise, at which point adiabatic conditions would obtain and the 
specific heat ratio of the gas used to estimate the‘ temperature 
at the instant of explosion, ignition delays again being disre- 
garded. It was soon observed, however, that very rapid 
pressure oscillations usually occurred at the start of the pres- 
sure-rise cycles. These oscillations could have*been antici- 
pated, since the confined gases on each side of the-piston could 
be expected to act essentially like springs. A typical oscillo- 
graph record of a blank test on water is shown in Fig. 3. The 
magnitude of the rate of pressure rise in the first oscillation 
may be as great as 10’ psi per sec, and it is not directly related 
to the average rate of pressure rise, which would be equal to 20 
times the rate of pressure rise in the driving chamber. Clearly, 
if an explosion did occur during the first half of pressure rise, 
the pressure and rate of pressure rise at that point could not 
be measured with any accuracy. It was concluded that heavy 
damping of the system would be necessary to produce smooth 
pressure traces, but if such damping were incorporated, ap- 
preciable deviation from the real situation would result. 

To assist in explaining the action of the tester, an expression 
was derived for the rate of pressure rise on the driving side of 
the piston in terms of the physical constants of the equipment 
and driving gas. It was assumed that flow from the high 
pressure reservoir through the gas release valve and flow re- 
stricting orifice was adiabatic, and that the temperature of the 
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gas in the driving chamber remained the same as that in the 
reservoir throughout the process. Also, it was assumed that 
the gas pressure in the reservoir was constant, and that 
driving-chamber volume remained constant. The gas was 
assumed to be ideal. 


Then 
n RoT 
P = 2 
[2] 
_ dng Ral _ 
dt dt V MV 
where 
Ng = number of moles of gas 
Ry = universal gas constant 
dP/dt = rate of pressure rise in driving chamber 
m = mass rate of gas flow into driving chamber 
M = molecular weight of driving gas 


Flow between the high pressure reservoir and the driving 
chamber will be at sonic velocity at the throat (the valve or 
orifice) until pressure in the driving chamber exceeds the 
critical fraction of pressure in the reservoir. For nitrogen at 
room temperature, the critical pressure ratio is 0.528. 


For sonic flow (2) 
A,Pogey 


[4] 
M 


n= 


where 


A, = throat area 

= reservoir pressure 

Jc = gravitational conversion factor 
Y = specific heat ratio 

7) = reservoir temperature 


It follows that 


dP = AP 2 y-1 [5] 
dt V M y¥+1 


Although all of the assumptions are not exact, the agree- 
ment between predicted and observed rates of pressure rise 
was found to be good. Pressure-rise rates in the test chamber 
were measured with the test chamber filled with water so that 
the piston did not move. Smooth pressure traces were ob- 
tained. The average difference between observed and cal- 
culated rates was less than 10 per cent, which was about the 
limit of accuracy in estimating the rate of pressure rise from 
the oscillograph records. The range of pressure-rise rates 
measured in the test chamber to check the derived expression 
was from 10,000 to 700,000 psi per sec. Because of the suc- 
cess of the correlation in this range, extrapolation to lower 
(4000 psi per sec) and higher (1,000,000 psi per sec) rates 
should introduce no serious errors. Plots of rate of pressure 
rise vs. driving gas pressure, using helium and nitrogen as 
working media, for a series of orifices of different areas were 
prepared and used in subsequent data reduction. 

If the pressure rise in the test chamber is rapid, the compres- 
sion will be nearly adiabatic and a rapid temperature rise will 
result. There is no conclusive evidence that high pressures 
or high rates of pressure rise can by themselves initiate ex- 
plosions, so that what probably is being tested is the sensi- 
tivity of propellant vapors to very rapid heating. It is 
probable that the maximum temperature reached during com- 
pression is the quantity of greatest importance in determining 
whether or not reaction occurs. For an ideal gas, temperature 
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reached in an adiabatic compression can be related to the 
energy going into the compression as 


PW 
oVo 
W= — 6 
[6] 


Since the compression is adiabatic and the gas is ideal 


(72) [7] 
Po T» 
and 
PoVo = NnghoTo [8] 
so that 
PoVo Ty NghoT’o ( 7) 
We —— (1-4) = - 9 
where W = work of compression (negative) or expansion 
(positive). 


It can be seen that, in an ideal system, the maximum tem- 
perature reached in a compression will be determined by the 
initial temperature, the specific heat ratio and the quantity of 
gas being compressed. Therefore, it seemed reasonable to 
expect that, if a definite minimum temperature would be re- 
quired for the initiation of a given fuel, the minimum energy 
input necessary to cause reaction in the rapid compression 
tester would be directly proportional to the initial gas phase 
volume of the sample, provided that the tester action actually 
approached adiabaticity. 

The energy input to the test chamber can be evaluated as 
exactly as the  PdV expression can be evaluated for the 
material in the test chamber. This was not practicable with 
the available instrumentation, as can be seen from examina- 
tion of Fig. 3. The same difficulties, in somewhat lesser de- 
gree, applied to the evaluation of  PdV on the driving side 
of the piston. 

To estimate the energy input to the sample, it was pro- 
posed that gas expanding from the high pressure reservoir 
imparts kinetic energy to the piston, and the kinetic energy 
of the piston is expended compressing gas in the test cham- 
ber. The analysis given by Lamb (3) for the compression of 
a spherical bubble surrounded by liquid appears to be relevant. 

The correctness of this picture depends on the fact that the 
opposing force of pressure in the test chamber is usually 
negligible until the piston has traveled over 90 per cent of the 
distance from its original position to the surface of the liquid 
sample, thus allowing the piston to acquire very nearly the 
velocity which it would attain if the sample chamber were 
open. 

An expression for the piston velocity was developed in terms 
of the rate of pressure rise on the driving side of the piston and 
the distance traveled by the piston. The observed result, that 
rate of pressure rise in the driving chamber was constant 
throughout the first half of the pressure rise, was used in the 
expression, since this means that the rate of change of ac- 
celeration of the piston will be constant. Denoting this con- 
stant by a 


[10] 


instantaneous acceleration 

rate of pressure rise on large end of piston 
area of large end of piston 

mass of piston 


ou 


integrating 
a = S = at + a 
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(neglecting gravity, a) = 0) 


velocity = v = S atdt ace + u% [12] 
z 


(since initial velocity is zero. v) = 0) 


at® 
=—+X 
6 + Xo [13] 


distance traveled = X = dt 
Combining Equations [10, 12 and 13] and inserting the 
necessary conversion factors and physical constants leads to 
the following expression for piston velocity 


2/, 
v= 14.9 (2) x! [14] 


where v is in in. per sec; dP/dt is in psi per sec, and X is in 
inches. 

Actual piston velocities were estimated from measurements 
on oscillograph traces and were found to average about 0.66 
of the predicted rates. Inserting the empirical correction 


factor gave 
dP\ _?/s 
v = 10.0 (=) X [15] 


Average deviation from the empirical equation over a range 
of measured velocities from 40 to 110 in. per sec was some- 
thing less than 10 per cent. Although the correction is based 
on a relatively small number of tests because of difficulties 
encountered in obtaining consistently legible records, it is 
thought that the available results are sufficient to verify the 
prediction that the piston velocity should be proportional to 
the 1/3 power of the rate of rise of driving pressure and the 
2/3 power of the distance traveled, and that any difference 
between calculated energy input and actual energy input is 
consistent, if not small. 

A model for the reaction initiation and propagation in the 
compression tester process has been suggested. Temperature 
in the gas bubble is raised by compression; the heat generated 
tends to evaporate liquid from the bubble surface, which 
tends to lower the bubble temperature. Simultaneously, the 
decrease in bubble volume tends to condense vapor at the in- 
terface. From an examination of vapor pressure vs. tempera- 
ture characteristics of several monopropellants, it appears that 
the net effect during this stage of the process will almost 
always be evaporation accompanied by a temperature rise, if 
most of the material in the gas phase is insoluble in liquid. 
This occurs when the gas consists of a mixture of fuel vapor 
and air. However, if the vapor pressure of the liquid were 
only a weak function of temperature, or if the gas were 
soluble in the liquid, the net effect could be condensation. 
Also, during compression, heat is being lost to the mass of the 
liquid sample and to parts of the equipment in contact with 
the bubble. If the compression is fast enough, these losses will 
be negligible. 

When the temperature becomes high enough, fuel in the 
gas phase decomposes, generating heat. Some of this heat 
can be lost to the liquid sample and tester parts, and some can 
go into evaporating more liquid. If enough heat is generated 
in the bubble, the reaction will become self-sustaining. 

The situation in the rapid compression tester is complicated 
by rapid oscillations in pressure, with corresponding oscilla- 
tions in temperature. However, these are the conditions to 
which an entrained bubble in an actual propellant line would 
be subjected by the sudden closing of a valve, in either a 
pressurized system or a centrifugally-pumped system. There- 
fore, it appears that the sensitivity of a fuel as measured by 
the adiabatic compression tester should be a good indication 
of the sensitivity of the same material in use. 
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Fig. 4 Compression sensitivity of n-propyl nitrate 


Test Procedure and Results 


Because a prohibitive amount of time and effort would 
have been required to obtain oscillograph records of any large 
number of tests, the correlations discussed above were used 
to evaluate energy input to the fuel samples. For the reasons 
previously discussed, it was thought that the energy required 
to initiate liquid phase explosion in a given fuel would be 
directly proportional to the gas phase volume of the sample 
if the tester action actually approached adiabaticity, so two 
or more bubble volumes were tested for fuels which could be 
made to decompose. Only two types of results were con- 
sidered: Positive, in which essentially complete decomposi- 
tion of the sample occurred; and negative, in which the fuel 
sample remained apparently unaffected by the test. Un- 
doubtedly, partial decompositions or vapor phase ignitions 
which did not propagate were obtained occasionally, but re- 
sults of this type could not be identified with any certainty. 

Because other characteristics of n-propyl nitrate have been 

ac well documented, it was thought that the results of intensive 
a testing of this material would be of practical interest. The 
presence of air in the bubble, because of the high specific heat 
ratio and the possibility of ignition by oxidation, is thought to 
represent the most severe condition which might be en- 
countered in practice. 

Bubble volumes of 0.2, 0.4, 0.6 and 0.8 ml were investigated. 
A “reversal” technique was employed, starting in the approxi- 
mate range determined by preliminary tests, increasing or de- 


creasing driving pressure in uniform increments according to 
the type of result obtained first, until a reversal in result oc- 
curred. The direction of pressure change was then reversed, 
and the procedure continued until the desired number of tests 
had been completed. A total of 28 tests was made for each of 
the four bubble volumes. The results are shown in Fig. 4. The 
average of positive results (after starting a “reversal’’ series) 
is indicated by “X”’; the average of negative results is indi- 
cated by ‘O” for each bubble volume. The accuracy of the 
slope is thought to be +15 per cent. The numerical value is 
6.7 kg-cm/ml. At 0.8 ml, positive results were occasionally 
obtained at energy inputs as low as 3.2 kg-cm, and negative 
results as high as 5.9 kg-cm, so it appears that in order to ob- 
tain a reliable value of the parameter, minimum ignition 
energy per unit bubble volume, it is necessary either to run 
several tests at one bubble volume, or a few tests at severa| 
bubble volumes. 

Satisfactory records of explosions were obtained while test- 
ing n-propyl nitrate. Shown in Fig. 5 is a typical record of an 
explosion caused by an energy input of 6.2 kg-cm to a sample 
consisting of an 0.9 ml air bubble in contact with 0.55 ml oi 
n-propyl! nitrate liquid. It is in agreement with the average 
value of a fairly large number of tests, and there is a clear in- 
dication from the piston displacement trace that rapid reac- 
tion started during the first pressure surge, as can be seen by 
comparison with Fig. 6 which is a record of a test made on 
water under the same conditions. Fig. 7 shows the effect of 
an energy input of 4.0 kg-cm to a similar sample of n-propy! 
nitrate and air bubble, while Fig. 8 shows the result of a blank 
test on water with the same conditions. From the cusp in the 
pressure trace and the high rebound of the piston, it appears 
that reaction started, but quenched rapidly. Fig. 9 is the 
record of an explosion caused by an energy input of 3.2 kg-cm 
to another sample consisting of 0.55 ml of liquid n-propy] 
nitrate in contact with an 0.8 ml air bubble. The result ap- 
pears unusual for several reasons. First, the explosion was 
caused by an energy considerably less than the average re- 
quired for these conditions. Second, there is no evidence of 
initiation of reaction at the end of the downstroke of the pis- 
ton. Why the explosion occurred is not known. It can only 
be repeated that the process of ignition is complex and, ad- 
mittedly, not very well understood. 

Nitromethane was tested in the presence of 0.2, 0.4 and 0.8 
ml air bubbles. At the two smaller volumes, only enough 
tests were made to locate the first reversal in result. At the 
0.8 ml bubble volume, a reversal series of 10 tests was made. 
The average obtained in this series was reasonably well in 
line with the results at the lower bubble volumes. Since in 
normal operation the direction of piston motion is straight 
down toward a stationary liquid surface, it was thought that 
a difference in results might be observed if the piston were 
driven upward, accelerating the liquid, and possibly agitating 
the liquid surface enough to form a large number of small 
bubbles. No significant difference was noticed. Results of 
the two types of test are shown in Fig. 10. The sensitivity of 
nitromethane in the presence of air measured by either type 
of test would be 10.4 + 1.7 kg-cm/ml. 


DISPLACEMENT 


Fig. 5 Record of explosion Fig. 6 Blank test record 
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Fig. 7 Record of probable 


: se Fig. 8 Blank test record 
partial decomposition 
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Fig. 10 Sensitivity of nitromethane 
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Fig. 11 Sensitivity of 60 per cent ethyl nitrate-40 per cent 
propyl nitrate 


The results obtained for a mixture of 60 per cent ethy] ni- 
trate and 40 per cent propyl nitrate in contact with air 
bubbles of various size are shown in Fig. 11. Also shown are 
results for the same material with fuel vapor bubbles of various 
size. This mixture was the most sensitive of any of the ma- 
terials tested, and was also the only one which could be made 
to decompose in the absence of air in the vapor space. 

The sensitivity of ethyl-propyl nitrate (and of the other 
fuels which could be ignited in the presence of air) is clearly 
much lower in the absence of air. Several explanations are 
possible for the decreased sensitivity. Among these are the 
absence of oxygen and the high specific heat of the fuel vapor 
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relative to air. Also, vapor in the presence of a liquid can 
condense very rapidly. For example, the nitromethane vapor 
in a 1-ml bubble, in contact with a 1-cm? surface of liquid 
nitromethane at 20 C could condense in a time on the order of 
10-4 sec. This would mean, in the compression tests with 
liquid in the presence of vapor only, that the fuel vapor is 
being condensed instead of being compressed. This condition, 
of course, would raise the temperature of the liquid surface, 
but the temperature rise would be very small compared to 
that produced in compressing a bubble of insoluble gas with 
the same amount of energy. 

Methylacetylene was made to decompose at bubble volumes 
of 0.8 and 0.4 ml. For these tests, the liquid sample was 
loaded at about —40 C and fired at about —30 C. Cooling 
was effected by circulating trichloroethylene through the 
tester jacket and a coil immersed in a dry ice-trichloroethylene 
bath. Negative results were confirmed by the presence of un- 
affected liquid remaining when the chamber was opened while 
still cold. Decompositions were indicated by the presence of 
a compacted carbon pellet in the chamber, and also by pres- 
sure in the chamber. The results are shown in Fig. 12. Only 
enough tests were made to locate the first reversal in result. 
Since maximum energy input possible with this equipment is 
not a linear function of bubble volume, there is no special 
significance attached to the absence of any positive result at 
0.2 ml. A test with an 0.8 ml bubble, in which the sample was 
sealed and allowed to warm to room temperature before 
firing, gave a negative result at maximum energy input. 
The vapor pressure of methylacetylene is about 90 psig at 
room temperature. 

Hydrogen peroxide, hydrazine, unsymmetrical dimethy] 
hydrazine and ethylene oxide (tested both at —20 C and at 
room temperature) could not be made to decompose by test 
conditions of maximum severity, i.e., an energy input of 115 
kg-cm to an 0.8 ml air bubble, or 144 kg-em/ml. As far as is 
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Fig. 12 Sensitivity of methylacetylene at —30 C 
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known, no accidental explosions which could be blamed on 
rapid compression have occurred in the handling or use of 
these materials. 


Conclusions 


The rapid compression test, in its present state of develop- 
ment, can be used to rate numerically the sensitivity of 
liquids to ignition in the presence of rapidly compressed air 
bubbles. The parameter on which the scale of rating is based 
is the minimum compressive energy per unit bubble volume 
necessary to cause gas phase ignition followed by propagation 
to the liquid. The sensitivity scale is thought to be valid and 
useful because of the simple relationship between the indicated 
sensitivities of different materials, i.e., if one material is 
shown to be five times as sensitive as another, it can be ignited 
by one-fifth as much compressive energy if other conditions 
are the same. Although the action of the tester is not simple, 
it appears to simulate closely conditions which might occur in 
actual propellant systems. For each of the four monopropel- 
lants which have been made to ignite or decompose in the 
presence of air bubbles, it has been found that the sensitivity 
parameter is constant over the range of bubble volumes (0.2 
to 0.8 ml) tested, within the limits of accuracy of the test. 

For the materials tested, the sensitivity ratings in the 
presence of air bubbles are as follows 


Monofuel Sensitivity, kg-em/ml 
60 per cent ethyl nitrate 
40 per cent propyl 
n-propy! nitrate 
nitromethane 17 
methylacetylene 86 + 12 
hydrogen peroxide >144 
hydrazine >144 
unsymmetrical dimethyl hydrazine >144 
ethylene oxide >144 


All of the fuels which could be ignited in the presence of air 
bubbles were substantially less sensitive in the presence of 
fuel-vapor bubbles. It is believed therefore that the tests in 
the presence of air bubbles represent the most severe condi- 
tions which might be met in practice. 

Because of the reasonable constancy of the sensitivity 
parameter over the range of bubble volumes tested, and be- 


cause of the larger range of energy inputs possible with the 
larger bubble volumes, it is suggested that the sensitivity of 
experimental monofuels which were available only in small 
quantity could be evaluated by a series of tests at the largest 
bubble volume practicable. Almost certainly, because of the 
rapid increase of surface to volume ratio with decreasing bub- 
ble volume, resulting in increased rate of heat loss, there will 
be some minimum bubble volume below which the indicated 
sensitivity parameter will increase very rapidly with decreas- 
ing bubble volume. Although this ‘critical’? bubble volume 
was apparently less than 0.2 ml for the monofuels which could 
be made to react in the tester, with the possible exception of 
the ethyl-propy] nitrate mixture in the presence of fuel-vapor 
bubbles, the possibility of testing experimental materials at 
bubble volumes below “critical’’ should be minimized. It is 
thought that a series of 20 tests, requiring 10 or 15 ml of ma- 
terial, would be ample to establish a reliable sensitivity rating. 

There is an evident need for the definition by example of 
the potentially hazardous range of sensitivity by the selection 


’ of two standard fuels, one of which would mark the dangerous 


limit of sensitivity and one of which would mark the safe limit 
of sensitivity. Until such a range has been established, the 
usefulness of the test results will be restricted, and any need 
for extending the range of the tester will be uncertain. 
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Isothermal Compressibility of 
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The methods for computing the isothermal compressibility for RP-1 and liquid oxygen are pre- 
sented. Liquid oxygen is a single substance, and RP-1 is a complex petroleum mixture, hence dif- 
ferent methods of approach were used in each case. The accuracy of the calculated result is on the 
order of the experimentally derived result. 


HE LIQUID propellant engine system generally con- 

sists of one or more tanks to store the propellants, a feed 
mechanism for forcing the liquids into the thrust chamber, a 
power source for furnishing the energy required by the feed 
mechanism and a control device for regulating the propellant 
flow rates. In order to arrive at a satisfactory feed mechanism 
design, correct power supply and accurate mass flow rate, it 
is necessary to know the density of propellants under actual 
operating conditions. One of the present trends of combus- 
tion chamber design in the development of high thrust rocket 
engines is to increase the chamber pressure. Therefore, the 
study of compressibility of liquid propellants is one of growing 
importance. 

In general, the compressibility of liquid is small, a few 
hundred-thousandths of the initial volume (1 X 1075) per 
atm. The difference in volume between that at 1 atm pres- 
sure and that under 1000 psi pressure is rather small, generally 
smaller than the accuracy with which a flowmeter can measure 
today. Since most rocket engines are now operated at a 
pressure of less than 1000 psia, it is doubtful that the pressure 
exerted on the propellants has any influence on the measure- 
ment of the specific impulse of the thrust chamber. However, 
the compressibility of the liquid propellants may have its ef- 
fects on combustion. Kretschmar (7)? has shown that the 
compressibility of propellants may be related to the pressure 
transients in the injector system, which, in turn, may couple 
with natural frequencies of the combustion chamber in such 
a way that instability and rough combustion result. Since 
one of the most commonly used propellant combinations con- 
sists of liquid oxygen and RP-1, a knowledge of their isother- 
mal compressibilities would be of interest to engine develop- 
ment work. This paper attempts to make a thorough evalua- 
tion of the quantitative relation of densities for these propel- 
lants under pressurized conditions, and to arrive at a method 
by which their compressibilities can be determined. 


Isothermal Compressibility of Liquid Oxygen 


Partington (10) made a comprehensive review of the litera- 
ture on the compressibility of liquids; no simple and decisive 
rule could be found to generalize the calculation of com- 
pressibility for liquids, even for a pure substance. 

There are two types of compressibility for liquids, isother- 
mal and adiabatic, as defined by the following equations 


1 (a 
(5), 
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1 ov 
v Op s 
where 
K; = isothermal compressibility 
K, = adiabatic compressibility 
s = entropy 
p = pressure 
v= volume 


There are more adiabatic than isothermal compressibility 
data in the literature. In cases where isothermal compressi- 
bilities are needed, some method to convert adiabatic com- 
pressibility into isothermal compressibility must be de- 
veloped. Such conversion may sometimes be difficult or even 
impossible without sacrificing accuracy. 

As shown by thermodynamics 


Ki 
[3] 
Tva? 
Crp = 
P K; [4] 
From Equations [3 and 4] 
Tva? 
K; — K. = — 5 
[5] 
where 
a = coefficient of liquid volume expansion 
C, = molar heat capacity at constant pressure 
C, = molar heat capacity at constant volume 
v= molar volume 
T = absolute temperature 


y = ratio of C,/Co 


Direct application of the above relations to compute KA; 
from Ka is rather difficult. y cannot be calculated from C, 
and C,, because C, of the liquid is almost impossible to meas- 
ure, and, in turn, is calculated from known values of C, with 
Equation [4]. 

In general, a for many liquids is known for different tem- 
peratures at 1 atm pressure. However, since a is a function 
of both temperature and pressure, values of @ at constant 
temperature will vary with pressure. For example, at zero 
deg C, a for mercury is 181 X 10~¢ at zero pressure, while a 
is 152 X 10-8 at 7000 kg/cm*. The difference due to pressure 
is around 20 per cent. In addition, C, and v are also functions 
of both temperature and pressure. For the above mentioned 
reasons, it is rather difficult to compute K; and K, without 
some approximation, unless the whole set of data on a, Cp and 
v is known. 
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aor 
F gm/em? em*/gm K-*x 
90 —297.7 1.142 0.8757 4.21 1.772 
85 —306.7 1.167 0.8569 4.15 1.722 
80 —315.7 1.191 0.8396 4.03 1.624 
75 —324.7 1.215 0.8230 3.94 1.552 
70 —333.7 1.239 0.8071 3.82 1.459 
65 —342.7 1.263 0.7918 3.65 1.332 
60.5 —350.9 1.282 0.7800 3.46 1.197 


Table 1 Calculation of K; for liquid oxygen 


T a?, Ta?/C>5, Ki, 
K-! x Cy, Cop, cm?/cal X in.?/lb X 
10-3 cal/gm/K cal/em?/K 10-3 10-6 
1.595 0.4059 0.4635 3.441 12.9 
1.464 0.4039 0.4714 3.106 11.6 
1.299 0.4019 0.4787 2.714 10 4 
1.164 0.4008 0.4870 2.390 9.30 
1.021 0.3989 0.4942 2.065 8.27 
0.8658 0.3972 0.5017 1.726 7.23 
0.7242 0.3972 0.5092 1.422 6.30 


In the case of liquid oxygen, only two point data (8, 12) of 
a can be found from the literature. Because these two values 
come from two different sources and only one value lies in the 
temperature range of interest, neither is adopted in this com- 
putation. a values in the whole temperature range and 1 atm 
pressure can be derived from the density data of Baley 
(1). The specific volume of the liquid oxygen is plotted vs. 
temperature on large-scale coordinate paper. Tangent lines 
are then drawn at various temperatures. In order to avoid 
errors in drawing these tangent lines, the slopes of the tan- 
gents are plotted vs. temperatures to smooth out the dif- 
ferences in drawing tangents. These smoother slopes (0v/Ot), 
at one temperature multiplied by the density of that particular 
temperature will give the value of a at that temperature. a@ 
values derived this way, in addition to the density data from 
Baley (1) and heat capacity at constant pressure data from 
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Giauque (3) are the only known data which can be utilized 
for calculating K; by Equation [5]. All these values are 
given at 1 atm. 

In order to use these data at 1 atm for computation at all 
pressures, some assumptions must be made. These in turn 
must be proved sound, and the results of calculations based on 
these assumptions must be checked with two known experi- 
mental values of isothermal compressibility of liquid oxygen. 

The logical approach to this problem is to assume (va?/C,) 
to be a constant at a definite temperature disregarding the 
change of pressures. The accuracy of this assumption can be 
proved in the following way: Since the compressibility of 
liquid oxygen is in the range of 1 X 10~*, the fluctuation of v 
due to change in pressure is rather small, less than 1 per cent 
from that at 1 atm, even at 1000 psia. As to the change of a 
due to change of pressure, the information can be estimated 
from Hougen and Watson’s chart “Thermal Expansion and 
Compression of Liquids’ (4). In the temperature range of 
— 297.7 to —350.9 F (90 to 60 K), the range of reduced tem- 
perature 7’, is from 0.396 to 0.583. At 7, = 0.396, the w 
value’ read from the chart is 0.1420 at 1460 psia and 0.1415 at 
zero pressure. The ratio of these two numbers is around 
1.0035. Even at the high temperature limit of the range, the 
corresponding values of w at 1460 and 0 psia are 0.1270 and 
0.1258. This ratio is 1.0095. The “Heat Capacity” chart of 
(4) also shows that the difference between C, at 1462 psia and 
C, at zero psia is around one-fifth of 1 per cent, in the tem- 
perature range of interest. Based on this information, one 
may be sure that the ratio a?/pC, is a constant with a maxi- 
mum error of less than 1 per cent on pressure change. The 
soundness of the above assumption can be further proved by 
comparing two calculated values with two experimental 
values 


Temperature, Experimental K;, Calculated K;,, 
F in.?/lb in.?/Ib 

— 299.7 1.25 <X 10-5 1.29 X 10-5 

— 345.5 0.75 X 1075 0.72 X 107-5 


Based on the above assumption, results calculated are listed 
in Table 1 and plotted in Fig. 1. For convenience of engineers 
in the field, a chart has been made which shows the precise 
density of liquid oxygen at various pressures with temperature 
and density as coordinates. 


Isothermal Compressibility for RP-1 


RP-1 is a special kind of kerosene used in rocket engines; 
its characteristics are regulated by military specification 
MIL-F-25576A(USAF), April 8, 1957. In the literature, 
there is some information about the compressibility of kero- 
sene (13) in general and about the jet fuels (7) in particular. 
However, such data can only provide the range of compressi- 


Fig. 1 Compressibility of liquid oxygen 


3 = expansion factor of liquid at P, and T,. 
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bility for RP-1, not the data which can be used for calculating 
the density of particular samples. Because RP-1 is a mixture 
not a pure substance, and the military specification gives 
limited information about RP-1, it is difficult to obtain the 
data which can be used for estimation even when supported 
with experiments on a limited number of samples. In such a 
case, the only way to obtain the necessary data is correlation. 

There are three generalized correlation methods in the litera- 
ture which can possibly be applied here. The ‘Thermal Ex- 
pansion and Compressibility of Liquids” chart of (4) probably 
is the most widely known. It covers a wide range of pressures 
and temperatures, and can be utilized for estimation up to 
critical temperatures of a liquid and up to a pressure five times 
the critical pressure of a liquid. This chart can also be easily 
used for interpolation. On the other hand, it has many de- 
fects. First, it is drawn from a reduced temperature of 1.0 to 
0.5 which then has been extrapolated down to 0.4. Informa- 
tion for liquid at temperatures lower than reduced tempera- 
ture 0.4 cannot be found from the chart. Secondly, the use of 
the chart requires knowledge of the critical temperature and 
pressure of the liquid under investigation. However, critical 
properties are not always known and are generally inaccurate 
for mixtures. This chart inherits the defects of all generalized 
charts, namely, inaccuracy. It is good for general estimating 
and engineering purposes, but not for precision work. 

Jessup’s correlation (6) of compressibility and thermal ex- 
pansion of petroleum published in 1930 is probably the best 
source of information which can be applied here. This 
method resulted from tests of 14 representative samples of 
petroleum oils from various sources over the pressure range of 
0 to 50 kg/cm? (gage) and the temperature range of 0 to 
300 C. 

Jessup found from the result obtained on these samples that 
the compressibility and thermal expansion of two samples of 
the same specific gravity, but from different sources, differed 
by more than 30 per cent at high temperatures; whereas oils 
of the same specific gravity and the same viscosity had the 
same compressibility and thermal expansion within rather 
narrow limits. 

The parameter used for Jessup’s correlation is the volume 
modulus 


1/(dV/ log 1000 s) 


= the specific gravity of oils at 60 F 
s = the kinematic viscosity at 100 F 


As he claimed, the maximum deviations of an observed point 
from the compressibility curve vs. the independent variable 
is around 0.1 per cent in volume at 50 kg/cm’. After com- 
paring his correlation with the experimental results of Zeit- 
fuchs (14) and Bearce and Peffer (2), he found that the agree- 
ment between his correlation and experiment was good. 

Jessup’s method was studied by Thiele and Kay (11). They 
reviewed the literature published before 1933 containing data 
on the expansion of liquid petroleum products at low tem- 
perature. They thought these experimental data adequate 
for application in ordinary engineering estimating purposes. 
However, these data cannot be considered sufficiently ac- 
curate in precision work. Jessup’s method was more suitable 
for extrapolation and interpolation, since it made the volume 
modulus relation at a given temperature a straight line. 
However, extrapolation to high values of modulus was not 
found satisfactory, even at 302 F, since there is considerable 
curvature at high modulus. At 167 F, the curves are straight 
for wide ranges of values, and provide a reasonable means for 
extrapolation. 

Thiele and Kay attempted to use the law of corresponding 
states to extend Jessup’s data to high temperatures. They 
presented a final drawing to show the density as a fraction of 
the density at 60 F and 1 atm, for gage pressures of 0, 400 and 
800 Ib per in.?, and for values of square of Jessup’s modulus up 
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to 2.4. Thiele and Kay also presented an alignment chart 
from which Jessup’s modulus can be quickly found. 

In 1946, Jacobson, Ambrosius, Dashiell and Crawford (5) 
presented a new method of correlation for compressibility of 
liquid hydrocarbons. The method used by Jacobson et al. is 
to plot API‘ vs. compressibility with temperature parameters. 
However, this method can be criticized. API is similar to 
density which, as shown by Jessup long ago, is not the only 
determined function of compressibility. Furthermore, the 
chart presented by Jacobson et al. is far from perfect. The 
only complete set of data used to determine the characteristics 
of the curve is at one temperature, 100 F. The other two 
basic curves at 80 and 30 F used to determine the space of 
temperature curves were based on calculated, not experi- 
mental, values. The reliability of such calculated values is 
doubtful. Despite the fact that Jacobson’s original chart is 
plotted on a contracted API scale, the points are still scattered 
around the curves. For this reason, the method of Jacobson 
et al. cannot be considered to give accuracy. 

Based on the above statement and discussion, the author 
concluded that the only method which can be adopted for 
estimating the density of RP-1 accurately is Jessup’s method 
and its modified version by Thiele and Kay. 

According to the military specification MIL-F-25576A 
(USAF), the density range of RP-1 is 0.801 to 0.815(60/60). 
The density range of RP-1 at other temperatures can be 
estimated by Hougen and Watson’s “Thermal Expansion and 
Compression of Liquids” chart. Because here only the range 
of density is required, and not the density of the particular 
sample, the accuracy of the chart is adequate. The maximum 
and minimum densities at 100 F estimated by this method 
are 0.800 and 0.785, respectively. 

The specification only gives the maximum viscosity of RP-1 
at —30 F as 13.9 centipoise. The maximum viscosity at other 
temperatures can be estimated with decane as a reference 
material by using the equation, log u = a log P + 6 (9). wis 
viscosity, P the vapor pressure, and a and 6 are constants. 
The maximum viscosity of RP-1 at 100 F estimated by this 
method is around 2.5 centipoise. Since the specification does 
not give the minimum viscosity, it can only be estimated from 
experience in handling RP-1 samples. Experience shows that 
for 100 F, a good guess for the minimum viscosity is 1.0 centi- 
poise. From the maximum and minimum viscosity and 
density values of RP-1 at 100 F, the maximum and minimum 
kinematic viscosities of RP-1 at 100 F were calculated to be 
3.185 and 1.250 centistokes. From the value of kinematic 
viscosity at 100 F and specific gravity at 60 F, the maximum 
and minimum “volume modulus” of Jessup were calculated to 
be 0.709 and 0.655. 

The narrow range of “volume modulus” of RP-1 further 
simplifies the application of Jessup’s method to RP-1. The 
above mentioned method fails at high temperatures when the 
compressibility-volume modulus straight line is extended to 
high volume modulus. However, this is not the case for RP-1 
within the narrow limit of volume modulus between 0.655 and 
0.709. The straight line relation can be sure to hold even at 
high temperatures for RP-1. In this case, Thiele and Kay’s 
modifications do not apply. 

The only chance for heating RP-1 in the rocket engine be- 
fore it is atomized from the injector is in the passages of the 
regenerative cooling system of the thrust chamber. Because 
of the short period of passage through the cooling jacket, the 
final temperature coming to the injector is never too high. 
This is generally below 572 F, the highest temperature which 
Jessup’s experiment covered. The highest pressure used in 
his experiment covers the pressure limit in the present-day 
rocket engine. The lower temperature limit of his data is zero 
deg C. In order to take care of the cooling of RP-1 by liquid 
oxygen in certain rocket system arrangements, the lower 
temperature of his data is extended to —40 F. 

4 Abbreviation of American Petroleum Institute gravity, related 
to specific gravity by deg API= [(141.5)/(spec. gr. at 60/60)]-131.5. 
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Fig. 2 Density of liquid oxygen under pressure 
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Table 2 Average compressibility of RP-1 at various 
temperatures under pressures up to 500 kg/cm? 


Compressibility < 108, in.?/lb, temperatures, F 
Volume 
modulus 32 86 167 302 437 572 


0.64 3.78 4.29 5.27 7.80 11.80 18.60 
0.66 3.82 4.34 5.38 8.02 12.15 19.45 
0.68 3.86 4.41 5.48 8.25 12.55 20.21 
0 70 3.89 4.46 5.60 8.46 12.95 21.00 
0.72 3.92 4.54 5.71 8.69 13.35 21.80 


From Jessup’s experimental data, mean compressibility vs 
volume modulus at various temperatures can be plotted. 
The mean compressibility in the range of volume modulu- 
between 0.64 and 0.72 is found for every 0.02 unit o/ 
volume modulus. For engineering application, the unit: 
of compressibility are converted into in.*/lb from cm?/kg, anc 
temperature units are converted into deg F from deg C. The 
result is summarized in Table 2. 

The data in Table 2 are plotted in Fig. 3. Fig. 3 can furthe: 
simplify the use of Jessup’s method for test stand application. 
At temperatures below 100 F the difference in compressibilities 
within the limit range of volume modulus is within the ac- 
curacy of his experiment. Therefore, it is accurate enough 
to take the average compressibility within the limit range of 
volume modulus at one temperature to represent the com- 
pressibility of all samples of RP-1 at that particular tempera- 
ture. Fig. 4 is a chart by which the volume correction for 
specific volume of RP-1 based on the specific gravity measure- 
ment at the measuring temperature at 1 atm can be found. 
It is the author’s opinion that any chart presented in 
this paper is limited only to correction, not the final density 
or specific volume. Because the correction is rather small com- 
pared with the specific gravity or specific volume itself, the 
wide range of specific gravity variation of RP-1 itself and its 
variation at various temperatures prevents the accurate 
readings of the specific gravity directly from a chart of reason- 
able size. However, the specific gravity of RP-1 under pres- 
sure can be easily calculated by the equation 
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Fig. 3 Isothermal compressibility of RP-1 


Fig. 4 Volume corrections for RP-1 under pressure, below 
100 F 
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S;,p = sp. gr. of RP-1 at temperature ¢ F and pressure p 
Ib/in.? 

S:1 = sp. gr. of RP-1 at temperature ¢ F and 1 atm 

v. = volume correction from Fig. 4. in X in.?/Ib 


7, = temperature correction from Fig. 4 at the measuring 
temperature 
= pressure of RP-1 in Ib/in.? 


The calculation can be faster with the help of a table of re- 
ciprocals of numbers, and the product v, X T, X p can be 
calculated with a slide rule without affecting the accuracy of 
the final result. 

At a temperature higher than 100 F, the simplified method 
above cannot be used without sacrificing a certain amount of 
accuracy, because of the increasing effect of Jessup’s modulus 
on compressibility with increasing temperature. For ex- 
aniple, at 100 F, the compressibility change within the range 
of Jessup’s modulus is from 4.42 to 4.70 X 10~* in.?/Ib, while 
at 550 F, the compressibility change within the same range of 
Jessup’s modulus is 17.3 to 20.0 X 10-6 in.?/lb. However, 
the utilization of the compressibility of RP-1 at high tem- 
perature is not for daily routine work at the test stand, but 
rather for application in designing when one has more time 
to calculate for the purpose of accuracy. The following equa- 
tion can be used in this case to calculate the actual specific 
gravity after the compressibility has been found from Fig. 3 


1 
— CX p) 


Si.» [7] 
where C is the compressibility in in.2/lb X 10~* from Fig. 3 at 
measuring temperature ¢ of RP-1 with its characteristic 
“volume modulus.” Density of RP-1 at 60 F and kinematic 
viscosity at 100 F are the additional data required for the 
calculation of the density of RP-1 at measuring temperature ¢ 
under pressure p. 


Conclusions 

Liquid oxygen is a single substance, while RP-1 is a complex 
petroleum mixture. This paper presents methods to calculate 
the isothermal compressibility in both cases. The same 
technique presented in this paper can be applied to investigate 
the isothermal compressibility for other propellants. The 
present approach serves to reduce the amount of experimental 
work through the use of data from published literature. How- 
ever, such published data should be subjected to analysis in 
order to insure that the accuracy of the calculated result is 
comparable to that derived from experiments. 
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An Attitude Control System for 


Space Vehicles 


WALTER HAEUSSERMANN' 


Army Ballistic Missile Agency 
Huntsville, Ala. 


An attitude control system, applicable to any size space vehicle, is proposed which uses expulsion 
charges to compensate for initial disturbances. Final fine control is exerted by a low power flywheel 
reaction-type servo system which makes use of a specially designed d-c servv motor with an imped- 
ance-type commutator and transistor control. The d-c servo motor is able to work in vacuum. 
It is shown how the control system can be stabilized without direct derivatives of the attitude in- 
formation; therefore, the system also supports stabilization when manually controlled. Some 
possibilities are mentioned to derive attitude information for space and Earth fixed reference sys- 


tems. 


HE ATTITUDE control system for a space vehicle must 
be designed to meet requirements (1,2)? which differ so 
much from those of the control system for the booster stages 
that independent control systems will be preferable for both 
applications. 
Presented at the ARS Semi-Annual Meeting, Los Angeles, 
Calif., June 9-12, 1958. 
1 Director Guidance and Control Laboratory, Development 


Operations Division. Member ARS. 
? Numbers in parentheses indicate References at end of paper. 
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Some of the major requirements for a three-axis attitude 
control system for a space vehicle are: 
Required features: 
1 Ability to compensate for the relatively high initial 
angular impulse resulting from the last separation process. 
2 Automatic control as well as assistance to manual re- 
orientation. 
3 Continuous operational duty must be possible. 
4 Control system must have a lifetime at least equal to 
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Fig. 1 Functional diagram of a one-axis control system 


the lifetime of the space vehicle as determined by its main 
missions. 

5 Necessary power should be as small as possible and 
be drawn on small satellites from solar cells delivering d-c, on 
larger vehicles from nuclear power supplies or solar cells. 

Desirable features: 

1 Control system should be able to work in vacuum. 

2 The principles of control should be applicable to small 
and large space vehicles. 

3 The system should be able to withstand accelerations 
during ascent of at least 6 g for manned satellites, and high 
accelerations during ascent as occurring on instrumented 
satellites. 


Proposed Attitude Control System 


Several methods of exerting control torques are possible. 
One which appears most attractive has been selected in the 
following description. It uses the inertial reaction torque 
principle. 

For the attitude control of each of the three mutually per- 
pendicular axes, a separate control system is proposed, since 
such a solution gives minimum cross coupling effects between 
the three individual torque producing systems. 

The principle of the scheme of such a one-axis control 
system is shown in the functional diagram Fig. 1. 

The attitude control signal as derived from a sun sensor, an 
infrared horizon sensor, a supervised gyro, or from a manual 
control, is amplified in a control computer and accelerates a 
d-c motor with a flywheel. The stability of this system can be 
improved if either the first derivative of the input attitude 
signal or the flywheel velocity is fed back. 

The requirement to overcome the initial unpredictable 
angular impulse, caused by separation from the booster, can 
only be met by the described angular reaction system if the 
integrated reaction torque or the corresponding flywheel rpm 
change is tolerable. In order to avoid unfavorable design 
limitations and large constant friction losses of the system, it 
is proposed that the initial angular impulse be almost com- 
pletely compensated by an expulsion or thrust device which is 
controlled by the flywheel speed. By proper staggering of the 
magnitude of the angular impulses and by selection of the 
number of steps, a considerable reduction in the torque re- 
quirements, and therefore in the overall flywheel system 
parameters, can be gained. The lower limit for the design 
parameters is finally given by factors, such as the minimum 
desired loop frequency, system sensitivity, transient disturb- 
ances due to moving parts or bodies in the space vehicle, etc. 

The proposed expulsion scheme has been selected to achieve 
simplicity for the sequencer and for the control of the expulsion 
charges. The magnitudes, or “impulse factors,” of the im- 
pulses are given in Fig. 2 with the time sequence in which they 
must be fired. Two sets of N expulsion charges are provided 
for plus and minus angular impulse directions. Regardless 
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Fig. 2 Magnitude and sequence of expulsion charges 


of the magnitude of the initial disturbance which must be /e- 
duced, N pulses will always be fired in the sequence shown on 
the time scale of Fig. 2. The magnitude of each individual 
charge is twice the magnitude of the following charge, and its 
sign will be determined by the direction of the flywheel rot:- 
tion. 

Thus the indicated impulse factors are multiples of the last 
charge impulse, which may be called “unit impulse.” This 
unit impulse must be small enough to diminish the angulur 
momentum of the space vehicle to a remaining value which 
can be transferred to the flywheel. Therefore, the maximu:n 
angular impulse, which finally must be compensated by an 
average speed of the flywheel, is equal to the unit impulse. 

The design reduction factor, which may be defined as the 
ratio of maximum angular impulse of the space vehicle to tlie 
impulse which can be compensated by the flywheel alone, is 
also indicated in Fig. 2. It is equal to the ratio of maximum 
allowable initial angular impulse to the unit impulse. 

It can be expected that the resulting flywheel speed will be 
reduced after some time due to eddy current damping of the 
flywheel if the space vehicle travels through magnetic fields. 
This effect should be considered in the design and layout of 
the flywheel system. 

The flywheel-motor arrangement shown in Fig. 1 demands 
special design features to fulfill the necessary and desirable re- 
quirements; in particular, the d-c motor and the d-c generator 
must be designed in such a way that they will work without 
mechanical commutators. In order to explain the proposed 
modifications, the motor design will be discussed first. It is 
obvious that a d-c motor can have its armature winding in 
the stationary part and its field excitation in the rotor, if the 
armature current is conducted through sliprings to rotating 
brushes which transfer the current to the stationary commu- 
tator, which is connected to the stationary armature winding. 
A brushless solution is obtained by replacing the “slipring- 
rotating-brushes-stationary-commutator” system with an 
“impedance-commutator-transistor switch” system in which 
the impedance commutator controls the transistor-switching 
due to its air gap variations, which are a function of the stator 
to rotor positioning. An oscillator is necessary to excite the 
impedance commutator; however, this control circuit re- 
quires practically no power because it only controls the opera- 
tion of the transistors which perform the high power transfer 
to the stationary armature. 

Fig. 3a shows a cross section of the impedance-commutator- 
motor assembly; the motor uses permanent magnets for field 
excitation in order to avoid sliprings, excitation power and 
heat dissipation problems. Fig. 3b indicates the air gap 
variation of the impedance commutator developed along its 
circumference. 

Fig. 4 gives the schematic diagram for a one-axis control 
system. The proposed impedance commutator is, in princi- 
ple, a generator with an amplitude modulated output. The 
demodulated output voltage gives the three-phase signal for 
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the control gate of the armature motor current. The armature 
carries the motor winding and a generator winding in the same 
slots. The a-c voltage, which is induced in the generator wind- 
ing, is proportional to the speed of the flywheel. The rectifica- 
tion of this voltage is controlled by the demodulated output 
of the impedance commutator which also controls the com- 
mutation of the motor current. The rectified voltage of the 
generator winding provides the velocity feedback signal used 
for damping and also furnishes information to trigger the se- 
quencer for the expulsion device through a controlled gate. 

The derivation in the next paragraph shows that the system 
is stable due to the effect of the induced emf of the motor. 
Therefore, in the case of some calculated examples with 
suitable numerical system data, no velocity feedback signal 
had to be applied to stabilize the control system. 


Control System Characteristics 


The following derivations are based on the assumption of 
an idealized control system, the components of which do not 
have any time lag, friction, saturation or other nonlinearities. 
Furthermore, coupling terms of the three mutually perpen- 
dicular control axes are considered to be negligible. The equa- 
tion of torques about the axis to be controlled gives 


+ he =0 1] 


J; = moment of inertia to be controlled 

I, = moment of inertia of rotor (including flywheel, etc.) 
@ = angular displacement of space vehicle, input signal* 
a@ = angular displacement of rotor® 

Dots written above variables denote time derivatives 


The equilibrium of voltages in the armature circuit of the d-c 
servo motor gives 


a) + — & + khd = 0 [2] 
where 
f(¢,a@) = control voltage applied to the armature of the 
servo motor 
= ad + a(h — &) according to Fig. 4 [3] 
e(@ — &) = emf of servo motor 
kho = ohmic voltage drop in armature circuit due to 


control current proportional to its correspond- 
ing torque lid 


The characteristic equation of the system of Equations [1, 2 


and 3] is 
s+ 0 [4] 


In order to describe the system completely, it should be 
mentioned that the characteristic equation actually contains 
a factor of s?, which has been canceled out on Equation [4]. 
The two poles at s = 0 allow for two initial conditions 
with respect to a and &. The initial displacement a of the 
flywheel is of no concern, whereas the initial value of the fly- 
wheel speed @ is zero in order to have symmetrical starting 
conditions for the operating range of the control system. 
Thus, the stability of the flywheel control system is sufficiently 
described by Equation [4]. 

Equation [4] indicates that the system behaves like a 
damped pendulum with the undamped natural frequency 


[4a] 


1 
and the relative damping ratio 


late 1 
2 Vkao/T; (- [4b] 


3 Measured in space direction fixed system, as desired by control 
mission. 
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Fig. 3a Impedance-commutator-motor 
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Fig. 3b Air gap variation of impedance-commutator-motor 
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Fig. 4 Schematic diagram for a one-axis control system 


The described motor works most efficiently like any d-c 
motor with a proportional control system. However, the 
available power sources are constant voltage sources. Thus, 
highly efficient control inverters and amplifiers have to be de- 
veloped. Most promising is the use of pulse modulated 
transistor and controlled rectifier-type amplifiers. 

The use of control characteristics with a properly selected 
dead zone or an on-off behavior will further improve the ef- 
ficiency of the control circuit, because no power is required 
during the dead zone interval. 

Preliminary design results are very favorable. The weight 
of all the flywheels for a three-axis control system contributes 
less than 0.5 per cent to the weight of the space vehicle. It 
can be expected that the total weight contribution, which in- 
cludes the weight of the control circuitry, the servo motors, 
and the expulsion devices (however, no weight of the 
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Table 1 Feasible methods for attitude sensing 


Medium to be sensed 


From 


Attitude sensor 


space fixed direction 


direction fixed to center of emit- 
ting body 

direction fixed to center of body 
shadowing cosmic rays 

magnetic field vector direction 


visible light 
visible light, infrared light 
cosmic rays (if no radiation belt 


present) 
geomagnetic field (application 


stellar bodies 
sun, planet, moon 
planet, moon 


Earth 


light direction sensitive cell 
(5,6) 

light and infrared direction sen- 
sitive cell (5,6) 

cosmic ray sensor (3) 


magnetic compass 


questionable due to el. currents 


in ionosphere) 


direction fixed to gravitational gravitational field 


direction determined by 
space vehicle to Earth trans- 
mitter or receiver line 


field vector 

direction fixed normal to orbit (Angular velocity of satellite due to attitude control di- rate gyro 
plane rection fixed to center of orbited body) 

relative velocity direction, ac- gas atmosphere angle-of-attack indicator, ac- 
celeration celerometer 

determined by signal, reference RF signals Earth RF receiver, either on Earth or 


orbited body satellite body or pendulum (4) 


in space vehicle (distance for 
attitude sensing is limited by 
length of base line) 


attitude sensing devices and the power supply) can be kept 
within 3 to 5 per cent of the vehicle weight. 

A motor with impedance commutator has been built and 
works satisfactorily. Analog computer studies have been 
carried out to determine optimum conditions for the combina- 
tion of the expulsion and the flywheel systems. The results 
are as expected and confirm the feasibility of the proposed 
attitude control system. 
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Fig. 5 Control system test table 


Attitude Sensing 


In order to attitude control a body in space, attitude angles 
must be sensed with respect to a selected reference system. 
Such a reference coordinate system can be established by 
means: 

1 Within the space vehicle. For this purpose inertial de- 
vices alone are available, which maintain a space direction 
fixed system. 

2 Outside the space vehicle, such as by any radiating 
field, potential field or field with matter. 

Due to limitations in the performance of instruments, the 
accuracy of inertial devices is not sufficient to retain their meas- 
uring system for a long time in space fixed position. In order 
to avoid drift of the measuring system, it is mandatory to 
supervise it as mentioned under item 2 above. Thus the iner- 
tial system will be of benefit if interruptions of the supervisory 
information will occur as due to shadowing effects, etc., or if 
by averaging of a noisy supervisory signal, better attitude in- 
formation can be obtained. 

Information sources outside the space vehicle can be used to 
establish space fixed reference systems, celestial body (includ- 
ing the Earth) aligned systems, or a combination of both. 
Some feasible attitude sensors for different outside informa- 
tion sources with the resulting type of coordinate direction 
are given in Table 1. 

The most favorable combination for a three-axis attitude 
sensing system can only be selected with the space vehicle 
mission, its path, the resulting accuracy requirements, weight, 
power and volume allowances in proper consideration. 

For certain space vehicle missions, it may be possible to use 
a two-axis attitude control system, and a limited angular rate 
around the third axis, the roll axis, can be tolerated, or a 
specified angular rate may be required. Then it becomes 
necessary to use an input signal to the roll axis control system 
which depends on the spin rate. In this case the attitude 
sensor will be replaced by an angular rate sensor, e.g., a rate 
gyroscope if a small angular speed is requested, or an ac- 
celerometer if a specified angular rate is demanded. The 
feedback loop, which is shown in Fig. 4, producing the angular 
rate damping terms becomes unnecessary. 


Appendix Simulation of the Attitude Control 
Loop 


It is desirable to investigate first the complete one-axis sys- 
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tem under the influence of expected initial disturbances and 
with respect to its overall behavior due to inherent and built-in 
nonlinearities. The frictionless support of the simulated 
moment of inertia of the space vehicle can be accomplished 
by an air bearing as indicated in Fig. 5. The mounting table 
represents the space vehicle and carries the complete one-axis 
control system, including batteries. The attitude signal can 
be derived from an optical pickup, such as in real application 
or from a simulated inductive pickup. The initial disturbance 
or the first angular impulse may be introduced by an air jet. 

Simulation of the three-axis control system behavior can be 
achieved with the complete control system mounted on the 
table, which is now supported by a spherical bearing without 
the cylindrical guide (see Fig. 5). The simulating device will 
serve to study coupling effects if large attitude angles occur, 


such as during the period of alignment to a reference system 
or during the period of changing from one reference system to 
another. 
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Application of Astronomical 
Perturbation Techniques to the Return 


of Space Vehicles 


ROBERT M. L. BAKER Jr.” 


Aeronutronic Systems, Ine. 
Newport Beach, Calif. 


An application of astronomical perturbation techniques is made to the return of space vehicles 
from satellite, lunar and interplanetary missions. Such techniques increase the accuracy of the 
computations and shorten the computational times required. The proposed variation-of-parameters 
perturbation method is illustrated by applying it first to the standard problem of a gravitational 
orbit perturbed by drag and, subsequently, by applying it to the inverse problem of a gravity-free 
drag orbit perturbed by gravity. The perturbation method is employed specifically in the analysis 
of the return from space voyages, and it is concluded from other comparative calculations that the 
perturbation computations are 3 to 10 times faster and are 11 to 165 times more accurate than the 


nonperturbation computations. 


HE MANY integration steps that are required in a 
numerical integration lead to the accumulation of end 
figure error. The usual re-entry trajectory calculations, i.e., 
step-by-step numerical integration of the total accelerations 
acting on the vehicle—known to astronomers as Cowell’s 
method—have been especially subject to this difficulty. 
Furthermore, the many integration steps involve too much 
calculation time. In the case of perturbation methods, one 
integrates only the departures from a reference orbit. Thus, 
only a few steps are necessary in those phases of the vehicle’s 
flight when the perturbative forces are small instead of the 
very many steps required when all the forces are included in 
the integrand. (It has been found that, in this case at least, 
the advantage of having fewer steps more than compensates 
for the disadvantages of using the perturbation formulas 
which are not always as simple or symmetrical as the non- 
perturbation formulas, see (11).%) 
The reference orbit that is to be perturbed is usually chosen 
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as the orbit that best depicts the gross motion of an object. 
All forces not taken into account in this gross motion are then 
categorized as perturbative forces. It is usually advisable to 
utilize as the reference orbit, an orbit that can be integrated 
analytically in closed form at the outset without requiring 
a step-by-step numerical procedure. 

These considerations immediately led astronomers to the 
utilization of the two-body Keplerian orbit as the reference 
orbit for the study of planetary perturbations. All forces 
exerted by drag, other planets, etc., are then incorporated as 
perturbative forces that tend to change the orbital parame- 
ters, such as eccentricity, semimajor axis, etc., of the two- 
body orbit. 

Comparative computations of the problem discussed in the 
fifth section of this paper have demonstrated that, in this par- 
ticular case, the perturbation method requires one-fifth to one- 
thirtieth fewer integration steps than does Cowell’s method, 


and, hence, proceeds 3 to 10 times faster. It should be noted 


also that the accumulation of end figure error cannot be over- 
come by simply increasing the number of digits carried in the 
computation, because such a procedure merely conceals the 
inherent round-off errors of any tabular data included in the 
integrand (e.g., atmospheric density): As pointed out by 
Dirke Brouwer (1) the end figure error of a double integration 
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increases in proportion to the three halves power of the number 
of integration steps; hence, the perturbation techniques pro- 
posed here should provide answers with between 5’ and 30° 
or between about 11 and 165 times more accuracy than would 
the commonly employed nonperturbation techniques. 

The need for extremely accurate and rapid computations of 
re-entry trajectories has arisen, in part, out of the require- 
ments for return from voyages into space, e.g., when the re- 
covery of a human being is involved, a fast, accurate orbit 
determination is mandatory. Even in other re-entry prob- 
lems, the computational efficiency usually afforded by per- 
turbation techniques (which have often been employed ad- 
vantageously by astronomers (9)) should also be taken ad- 
vantage of by missile designers. 


Perturbation Analysis 


In essence, the following perturbation analysis utilizes the 
technique of the variation of parameters,* which involves the 
continuous correction (“rectification”) of a reference orbit 
such that at every instant the reference orbit is tangent to the 
true orbit. This technique can be best illustrated by ‘first 
applying it to the standard problem of a purely gravitational 
two-body orbit perturbed by drag, and subsequently by ap- 
plying the technique to the inverse problem of a rectilinear 
gravity-free, drag orbit perturbed by gravity. Fig. 1 and the 
following two subsections illustrate a comparison of these two 
applications. As the vehicle’s path is perturbed, its motion is 
represented by a series of osculating or continuously varying 
reference orbits. Perturbative forces cause a change in the 
parameters that specify the orbit and, consequently, cause a 
progression from one osculating orbit to the next. 

The notation and formulation of the following subsections 
are based largely upon the work of Herrick (2). 


Drag Perturbation of a Drag-Free, Gravity Orbit 


If we denote the perturbative derivative by the grave sym- 
bol (*), and define 


§ = the speed of the object tangential to its trajectory (in 
terms of surface circular-satellite speed) 

the radial distance of the object from the geocenter (in 
Earth radii) 

a = the semimajor axis of the osculating ellipse 

&* = the perturbative acceleration tangential to the trajectory 

of the object (e.g., resulting from drag) 
a’ = the perturbation in the semimajor axis 


ll 


then the perturbative derivative’ of the well-known vis-viva or 
energy integral 


2 1 
[1] 


becomes 
= /a*® [2] 
From Equation [2] we find that 
a’ = [3] 
‘ Instead of a variation-of-parameters technique, one could also 


profitably employ an “Encke” perturbation technique (see (8)). 
5T.e., if f is a function of the ane Near variable 7, and cer- 


tain variable parameters Po,... Pi. .., then 
of of dP; 
dr or 5 > oP; dr 


Symbolically, the terms of this equation are identified by df/dr = 
_+f*, where the parameters P; vary only because of perturba- 
tive forces. See chapter 17 of (2). 


GRAVITY -PERTURBED 
RECTILINEAR ORBIT 


DRAG-PERTURBED 
ELLIPTICAL ORBIT 


INITIAL ORBITS 
PERTURBED OSCULATING ORBITS ——— — 
Fig. 1 Variation-of-parameters technique 


Taking the derivative of the semi-lactus rectum or parameter 
p where 


p = a(1 — e?) [4] 
we have 
= — e?) — 2aee* [5] 


where e = the eccentricity of the osculating ellipse. 
By the law of areas (or the principle of the conservation of 
angular momentum) 


p = rX(rv)? [6] 


where v = the angular speed of the true anomaly v. (For 
geocentric orbits » is the angle between the object and the 
perigee measured, at the geocenter.) 

Again perturbative differentiation yields 


p» = [7] 


but, since drag is a perturbation that is tangential to the 
trajectory and, furthermore, since it is the only perturbative 
force to be considered, we find that 


rv’ = cos @ = 8 rv/s [8] 
where @ is the angle of elevation, i.e., the angle between the 
tangent to the trajectory and local horizontal plane. 

Hence, Equation [7] becomes 


= 2ps* /8 [9] 
If we employ the equation of the osculating ellipse 
p/r =1+ecosv [10] 
and solve for e’, then we find that 


23* 
es = (e + cos v) [11] 


Equation [11] agrees with previous computations made by 
Moulton (3) and Herrick (2). (Moulton’s more complex 
equation reduces to [11] if Eqs. [1 and 10] of this paper are 
substituted into the fifth equation of problem 4, p. 405 of (3).) 
Similar expressions for v‘, etc., may be obtained in a straight- 
forward manner. 


Perturbations of a Gravity-Free, Drag Orbit 


In the case of re-entry orbits, the drag forces predominate, 
and it is therefore advisable to discard the two-body orbit and 
substitute in its stead a purely gravity-free, drag orbit.6 Be- 
cause all the drag forces are directed opposite to the velocity 
vector, the orbit will be reetilinear. The orbital parameters, 
perigee distance’ and the longitude of perigee @ vary only 


6 The changeover from the gravitational to the drag reference 
orbit should take place when the drag accelerations exceed the 
gravitational accelerations. 

7 The distance from the center of the Earth to the point of 
nearest approach of the orbit to the Earth. 
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because of gravitational perturbative forces. Integrations of 
these parameters yield a complete history of the trajectory 
(see (4)). Lift forces, normal to the trajectory, also can be 
easily included in the new perturbation scheme. 

It is necessary to consider the rectilinear orbit only with 
respect to a spherically symmetrical Earth; therefore, the 
one angle @ and one distance z are the only parameters needed 
to specify the orbit. (Of course, considerations of the aspheric- 
ity of the Earth, rotation and periodic fluctuations in the 
Earth’s atmosphere, and gravitational anomalies can be in- 
cluded also in the perturbation analysis.) Another quantity 
(e.z., y or v) is required to give the position of the object on 
this orbit. 

Consider the coordinate system shown in Fig. 2. The 
rather unconventional labeling of the coordinate axes is em- 
ployed in order to preserve the direction of the z-axis toward 
perigee. 

Let the grave symbol (*‘) also denote perturbative variation 
of a quantity with respect to time in which the variation re- 
su ts from gravity. Thus, the radial and transverse perturba- 


tive accelerations are 
= —1/r? rv» = 0 [12] 


Because the osculating orbit is by definition the one along 
which the vehicle is instantaneously traveling, there exists no 
instantaneous change in the radial distance r and the true 
longitudinal of the vehicle 1; therefore, since 1 = v + @, we 
find that 


S=vS+6S=0 or = [13] 


Of course r? = 2? + y*;z = rcosv, and y = sinv. Accord- 
ingly, we find that 


zs = = [14] 


Through reference to Fig. 2, we observe that * = rv tan v. 
Therefore, since rv’ = 0 


= (rv sec? [15] 


1 
(=) 4 [16] 


But cos v/rv = 1/y, whereas r cos v = x so that 


and, consequently 


L 
(see footnote 8). 
Thus 
= zy/yr* [18] 
But 
y y 


Since dr/dr = =x, and dx/dr = (dx/dy)(dy/dr) X 
(x — @), the perturbation derivatives with respect to y as the 
independent variable (denoted by’) are, similarly 

= + [20] 
and 

x’ = xy/(y + (21) 


Let us examine the behavior of Equations [17, 18 and 19] in 
two limiting cases. First, assume that the object is approach- 
ing the Earth along a rectilinear orbit through the geocenter and 
that there is no drag. In this case dy/d tr = —y*r’, while x 


§ It should be recognized that y + y® is the speed of the object 
with respect to the perturbed z-axis. 
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x 


FIXED ARBITRARY DIRECTION IN 
SPACE [i.e., SPECIFIES AN 
INERTIAL REFERENCE FRAME) 


Fig. 2. Coordinate system 


and, consequently, x‘ and ®* are zero. Second, assume that 
the object is circling the Earth in a drag-free environment. 
In the case of a circular satellite y = 1/+/r andz = r, there- 
fore, ® = 1/r*/*, and z= 0. From another point of view, the 
x-axis would be rotating around the Earth just fast enough to 
keep pace with the object, so that the speed relative to the 
perturbed axis is zero; hence, y + y\ = 0. But, as can be 
easily shown (see Eq. [14]), ys‘ = —z@*, therefore y = 2@S = 
1/-/r. All of the preceding results are in agreement with 
conventional theory. 


Inclusion of Lift in the Gravitational 
Perturbation Scheme 


In certain high speed re-entry problems, the vehicle may 
be expected to exhibit lift as well as drag. If we denote the 
lift acceleration (directed normally to the velocity of the ve- 
hicle) by L, then, from Fig. 3. 


*’ = Lcosv — 1/r? [22] 

ru’ = —Lsinv [23] 

Also tan v = rv/i, or i= rv tan v, so that differentiation yields 

= rv tan v + rv sec’. Substituting the previously 
determined values of ** and rv’, we obtain 


1 
Leosv — > —L sin v tan v + rvv’ sec? v [24] 
r 


Solving for v*, we find that 


cos? v 
L cosv — [25] 
v = = 
Tv 
or, in terms of x, y and r 
v = (L — 2/r*)/y [26] 
From Fig. 3 we observe that x = rcosv; therefore 
z\=rsinvs= —yv~ 
Hence 
y\r 
whereas 


Of course, y* remains unchanged. 
The perturbational derivations with respect to y as the in- 
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Fig. 3 Coordinate system for the incorporation of lift 
dependent variable (denoted by’ ) are 


1 x 
2 
y 
=—*__(— 29 
Formulation of the Drag-Perturbation Equations 


We have determined the perturbation variations in the 
parameters specifying the rectilinear gravity-free orbit in the 


dy 
dr 


preceding sections and it now remains to incorporate 
these into equations suitable for computation. Since all the 
motion of the vehicle is directed along the y-axis, the basic 
drag equation of motion is given by 


dy C pA opoV eg? 
dr 250mo “fal 
[30] 
A —D,oypay? — -— 
r3 
where 


Cp, = the reference value of the drag coefficient 
Ao = the initial projected frontal area of the vehicle 
mo = the initial mass of the vehicle 


= — 


D,? y 


y 


go = the acceleration of gravity at unit distance (979.81979 
cm/sec?) 

fo = thesea-level atmospheric density (1.225 X 10-3 gm/cm?) 

= p/p 

=Cd/Cr 

= 

a =A / Ay 

y = the speed of the vehicle with respect to atmosphere 
p, Cp, A and m represent, respectively, the true value of the 
atmospheric density, drag coefficient, projected frontal area 


and mass of the vehicle. 
Actually, y can be expected to be a function primarily of 
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(4 


0 B 
dy ( =) y (4 *) 


the three variables 7',, o and y (where 7’, is the surface tem- 
perature of the vehicle). These variables will determine the 
transitional and continuum variation of the drag coefficient. 
In certain re-entry problems in which the transitional 
variation and compressibility variation of Cp (and hence of 
y) oceur at different points along the trajectory, it is possible 
to separate y into two factors® 


= [31] 


(oc) specifies the variation of the drag coefficient during tle 
transition from free molecule to continuum flow, and ‘s 
strongly a function of the local molecular mean-free pat',, 
and consequently, is a function of the density ratio a (5). y(.') 
specifies the variation of Cp with the Mach number, aid 
hence is properly a function of y (assuming that the local spec d 
is nearly constant). 

If it is assumed that wa remains constant and equal to 
unity during re-entry’ and that the gross variation of oy(c) 
and y(y¥) (for supersonic speeds), respectively, can be repr:- 
sented by 


ay(c) => Co Cor? + [32] 
and 
¥y) =A+B/y +... [335] 


(where for Cpp = 0.92 and mean sonic speed = 320m/sev, 
A = 1 and B = 6.84 X 107 for a fit to the data of (12)) 


then 


B B 
y 


One can proceed by either integrating with respect to time 7, or 
distance, y. For nearly circular satellite orbits (in which dy/dr 
and y are close to zero) direct integration of Equations [17, 18 
and 34] with 7 as independent variable is advantageous. In 
this case, one integrates the first term in Equation [34] 
analytically (as per Eq. [42]—with P = 0) and integrates the 
second and third perturbational terms of Equation [384] 
numerically to obtain y; a second integration of dy/dr = 
+ — then required to obtain y. On the other 
hand, for either large 7? or small z? (i.e., for either high speeds, 
e.g., meteorites, or large angles of elevation, 6) the employ- 
ment of y as independent variable is advantageous. In this 
latter case one proceeds as follows 


y 


-(4+9)] 
| v5) (4 + —) | + 


where dy/dr = y + ys = y — x*/yr’, and we define 


2 2 2 \2 


Since this last equation includes all of the perturbationa 
terms, the basic integral for 7 becomes simply 


y ¥ dP 
yo yo y 


where the zero subscript denotes the initial values of position 


;, (Ay + B) 


® In meteoritic analysis, it is convenient to include also a factor 
yap that specifies the variation of Cp which results from the 
reaction of ablated molecules, i.e., the “reverse rocket’’ effect. 
See, e.g., (6). 

10The variation of wa has been included, however, in the 
referenced works (6). 
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(y) and speed (y). Integration yields 
B B 
=| Yo + — | exp| —AD,? y(ao)ody — AP | — — [38] 
A A 
Since r? = x? + y*, we can in turn represent oy(a) by 
oy(o) = Yo + Yoy? + +... [39] 
where 
Yo 
Yo = + + +... 


Y, = Cy + 3Cexr? + +... 
Ye = Ce + 4Cgr? + + .. ., ete. 
| :mploying the summation notation, we will set 
i=n 
i=0 
anc observe that 
y | 
—AD,? ly = 41 
[vom + il [41] 


in which it is recognized that Y»; is a function of the parameter 
In practice it is found advisable to change the constants 
Co, Co, Cs, .. ., Within a few levels of the atmosphere, in order 
to .void an excessive number of terms. 
Vinally then, we can establish the equation for y as 


B |v | B 


Thus given y and P (and the initial conditions y and yo), ¥ 
can be explicitly determined in closed form. 

The computation with y as the independent variable 
proceeds by integrating 


nt dy [43] 
[44] 
y 
-f = dy [45] 
yo dy 


where the initial value of P is set equal to zero. 


Applications and Conclusions 


The perturbation and computational approaches of the 
foregoing sections are justified only if they lead to faster and 
more accurate calculations. Consequently, other comparative 
calculations were performed in which the conventional inte- 
gration of all the forces (Cowell’s method) and the perturba- 
tion scheme were employed on the same problem. The prob- 
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lems chosen for analysis are discussed in (6 and 10), and 
the Runge-Kutta integration method coupled with the Adams 
correction procedure (7) was employed in the actual computa- 
tions. 

In every computation of this particular problem, the con- 
ventional calculations (even though their equations were 
slightly simpler and more symmetric) required from 3 to 10 
times more computer time than did the perturbation computa- 
tions; furthermore, the many extra integration steps involved 
in the nonperturbation methods caused a gradual accumula- 
tion of end figure error. It is concluded, therefore, that there 
exist definite advantages to the application of astronomical 
perturbation techniques to the computations of the return 
trajectories of space vehicles. 
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dependent variable (denoted by’ ) are 


1 £ 
y Bes 
Formulation of the Drag-Perturbation Equations 


We have determined the perturbation variations in the 
parameters specifying the rectilinear gravity-free orbit in the 


dy 
dt 


preceding sections and it now remains to incorporate 
these into equations suitable for computation. Since all the 
motion of the vehicle is directed along the y-axis, the basic 
drag equation of motion is given by 


2 
] oyHay? — 
r3 
y [30] 
A —Dioypay? — 


where 


Cp, = the reference value of the drag coefficient 
Ao = the initial projected frontal area of the vehicle 
mo = the initial mass of the vehicle 


= - ( 


A+— 
y 


go = the acceleration of gravity at unit distance (979.81979 
em/sec?) 

fo = thesea-level atmospheric density (1.225 X 10-3 gm/cm) 

= p/p 

=Cd/Co, 

= m/m 

a / Ao 

y = thespeed of the vehicle with respect to atmosphere 


p, Cp, A and m represent, respectively, the true value of the 
atmospheric density, drag coefficient, projected frontal area 
and mass of the vehicle. 

Actually, y can be expected to be a function primarily of 
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the three variables 7',, o and y (where 7, is the surface tem- 
perature of the vehicle). These variables will determine the 
transitional and continuum variation of the drag coefficient. 
In certain re-entry problems in which the transitional 
variation and compressibility variation of Cp (and hence of 
y) occur at different points along the trajectory, it is possible 
to separate y into two factors® 


= [31] 


(co) specifies the variation of the drag coefficient during tlie 
transition from free molecule to continuum flow, and is 
strongly a function of the local molecular mean-free pat'i, 
and consequently, is a function of the density ratio ¢ (5). (i) 
specifies the variation of Cp with the Mach number, aid 
hence is properly a function of y (assuming that the local spec d 
is nearly constant). 

If it is assumed that wa remains constant and equal ‘o 
unity during re-entry’ and that the gross variation of oy(.) 
and y(y) (for supersonic speeds), respectively, can be repr 
sented by 


oy(a) = Co + Cor? + Curt +... [32] 
and 
Vy) =A+B/y +... [3:5] 


(where for Cpp = 0.92 and mean sonic speed = 320m/se«, 
A = 1 and B = 6.84 X 107? for a fit to the data of (12)) 


then 


B B y : 
y y 


One can proceed by either integrating with respect to time 7, or 
distance, y. For nearly circular satellite orbits (in which dy/dr 
and y are close to zero) direct integration of Equations [17, 18 
and 34] with 7 as independent variable is advantageous. In 
this case, one integrates the first term in Equation [34] 
analytically (as per Eq. [42]—with P = 0) and integrates the 
second and third perturbational terms of Equation [84] 
numerically to obtain y; a second integration of dy/dr = 

‘= y — xis then required to obtain y. On the other 
hand, for either large y? or small z? (i.e., for either high speeds, 
e.g., meteorites, or large angles of elevation, @) the employ- 
ment of y as independent variable is advantageous. In this 
latter case one proceeds as follows 


A+ 3) 
y 


where dy/dr = y¥ + y = y — 2*/yr', and we define 


B (: — EE + (+) +.. [36] 


y 


Since this last equation includes all of the perturbationa 
terms, the basic integral for y becomes simply 


dy y ¥dP 
yo (Ay + B) yo 


where the zero subscript denotes the initial values of position 


® In meteoritic analysis, it is convenient to include also a factor 
yep that specifies the variation of Cp which results from the 
reaction of ablated molecules, i.e., the “reverse rocket’’ effect. 
See, e.g., (6). 

10The variation of wa has been included, however, in the 
referenced works (6). 
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(y) and speed (y). Integration yields 


B B 
E + exp| f — P| [38] 
yo 


Since r? = xz? + y*, we can in turn represent oy(a) by 
= Yo + + +... [39] 
where 
Yo =C,+ 
Y, = C. + + +... 


Y, = Cy + 3Cor? + +... 
Ye = Cg + 4Cgxr? + +. . ., ete. 
|:mploying the summation notation, we will set 
i=n 
oy(o) = >> Yay 
i=0 
an observe that 
y ims |v 
—ADy* = —AD*)) [41] 


0 +=0 


in which it is recognized that Y2; is a function of the parameter 
In practice it is found advisable to change the constants 
(, C2, Cs, .. ., Within a few levels of the atmosphere, in order 
to avoid an excessive number of terms. 
Finally then, we can establish the equation for 7 as 


B |v ] B 


Thus given y and P (and the initial conditions y and yo), ¥ 
can be explicitly determined in closed form. 

The computation with y as the independent variable 
proceeds by integrating 


x'dy [44] 
ye 
y 
r -f = dy [45] 
ye dy 


where the initial value of P is set equal to zero. 


Applications and Conclusions 


The perturbation and computational approaches of the 
foregoing sections are justified only if they lead to faster and 
more accurate calculations. Consequently, other comparative 
calculations were performed in which the conventional inte- 
gration of all the forces (Cowell’s method) and the perturba- 
tion scheme were employed on the same problem. The prob- 
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lems chosen for analysis are discussed in (6 and 10), and 
the Runge-Kutta integration method coupled with the Adams 
correction procedure (7) was employed in the actual computa- 
tions. 

In every computation of this particular problem, the con- 
ventional calculations (even though their equations were 
slightly simpler and more symmetric) required from 3 to 10 
times more computer time than did the perturbation computa- 
tions; furthermore, the many extra integration steps involved 
in the nonperturbation methods caused a gradual accumula- 
tion of end figure error. It is concluded, therefore, that there 
exist definite advantages to the application of astronomical 
perturbation techniques to the computations of the return 
trajectories of space vehicles. 
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Technical Notes 


Thrust of a Conical Nozzle’ 


ELLIS M. LANDSBAUM? 


Jet Propulsion Laboratory, California Institute of 
Technology, Pasadena, Calif. 


The equations for calculating the thrust of a conical 
nozzle, originally presented by Malina, are reviewed. The 
correct use of the equations is demonstrated. The usual 
improper use of these equations results in a slight error. 
Although the error is negligible, the correct method should 
be more widely known. 


HE THRUST obtained from a conical nozzle is usually 
represented by 


F = Amy, + (Pe — Pa)Ac . 


where 


[2] 


2 
and a is the half-angle of the expansion cone. The other 
symbols have their usual meaning and are further identified 
in the nomenclature. Malina* correctly presented the 
original derivation of both equations. Usually, the values of 
v, and P,, for use in Equation [1], are based on the plane area 
and the assumption of axial flow at the exit. This is incorrect, 
and, although the difference in values is not large, the cor- 
rect values of v, and P, are based on the spherical area ob- 
tained from the assumption of radial flow at the exit. Since 
the above statements are not clear without reference to the 
derivation of the equations, the derivation, essentially as 
originally presented by Malina, is repeated in this paper. 


Determination of Thrust 


The thrust propelling a rocket may be related to the mo- 
mentum and pressure forces acting in the axial direction at 
the nozzle exit and assuming uniform axial flow 


F = mu + (Pe — Pa)Ae [3] 


As P, the atmospheric pressure is unaffected by the flow 
conditions, it will be neglected in the remainder of the 
derivation. 

At the exit of a conical nozzle, represented in Fig. 1, the 
flow conditions will be neither axial nor constant across nozzle 
exit A,, and 


F,= (Mprp)z + PyldAp [4] 


The variables cannot be integrated in closed form, and nu- 
merical integration is necessary. Since there is no net force 
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Fig. 1 Radial flow of a conical nozzle 


across the control surfaces A, and A,, the thrust can also »e 
determined by integrating across area A, 


= + As 5] 


If radial flow is assumed to exist at the exit, along the spherical 
area A,, the variables of Equation [5] are all constant, and 
the following equations can be written 


M=—= | 
A, Ps¥s 
(Myvs), = cos [7] 
(P,), = P; cos 8 [8] 
dA, = 2rR? sin 6d0 [9] 
F, = + sin cos 6d0 [10] 
F,, = + ps] sin? [11] 
Ay = TR? sin? a [12] 
A 
= [ i 42 P.Ay| [13] 
Ap “ (1 + cos a) [14] 
A, 2 
finally 
F, = + P,A, [15] 
or 
F, = (mv, + P;As) [16] 


It is evident that the values of v, and P, for use in Equation 
[1] are those which are related to the area A,, not the area A >. 

The calculation of the thrust coefficient, for use in the fol- 
lowing equation 


F, = Cr,P-A: [17] 
is then 


1+ cosa 


(Cry)ap 2 


Eprtor’s Note: The Technical Notes and Technical Comments sections of ARS JouRNAL are open to short manuscripts describing 
new developments or offering comments on papers previously published. Such manuscripts are published without editorial review, 
usually within two months of the date of receipt. Requirements as to style are the same as for regular contributions (see masthead page). 
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Fig. 2. Correction factor for divergence loss of a conical nozzle, 
y = 1.22 


It is readily seen that the use of equations such as 

F, = ACr,P,At [18] 
or 

F = At 


where the thrust coefficients are calculated in the usual man- 
ner, results in a slightly greater error than does the incorrect 
use of Equation [1]. 

The correction term is a function of €, a and y, since the 
change in thrust coefficient between A, and A, is a function of 
these variables. In order to use the vacuum thrust coefficient, 
related to Ap, the following can be written 


(1 + cos @) 


2 


(Cry)apPcA t 


The value of \; is presented in Fig. 2 as a function of the area 
ratio € and expansion angle @ for a gas with y = 1.22. De- 
creasing ¥ results in a slight increase in \;. The parameter ); 
represents the error resulting from the use of Equation [18]. 


Conclusions 


If nozzle optimization is redone using the equations cor- 
rectly, a nozzle of slightly smaller area ratio and slightly 
larger expansion angle will be obtained. Since the variation 
of \; with € and a is very small, the difference will probably 
not be noticeable. The error from neglecting ),; in determin- 
ing the thrust of a conical nozzle is small over the range of 
exit areas and expansion angles usually of interest. This will 
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result in a slightly more conservative design, since ), is 
greater than one. However, one should be aware of the cor- 
rect calculation for determining the thrust of a conical nozzle. 


Nomenclature 


area 
thrust coefficient 

thrust 

mass flow per unit area 

mass flow 

absolute pressure 

velocity of exhaust gas 

expansion angle 

specific heat ratio 

area ratio 

thrust coefficient correction term = (1 + cos a)/2 
thrust coefficient correction term 

density of gas 


PR WOM 


>a 


Subscripts 


atmosphere 
exit 

plane 
spherical 
vacuum 


uo 


Determination of Transient Pressure 
Flow Relationship by Momentum 
Measurements’ 


EDWARD J. CROKE? and JERRY GREY® 


Princeton University, Princeton, N. J. 


A momentum-sensitive device was used to measure the 
transient flow rate through a rocket motor injector. The 
relationship between instantaneous pressure drop and flow 
rate was measured under conditions of artificially induced 
small sinusoidal oscillations at various frequencies. The 
observed phase shift of flow rate behind pressure drop 
agreed quite well with previous theoretical predictions 
(1,2),4 but some discrepancy was observed between the 
predicted and observed amplitude ratios. 


HE PROBLEM of instantaneous liquid flow rate measure- 

ment in transient systems has appeared many times in the 
literature. One rather specialized need for this measurement 
became apparent in 1952, in connection with rocket combus- 
tion instability studies, in which it was necessary to measure 
the instantaneous efflux of propellant from the injector orifices 
of a liquid propellant rocket motor. This need resulted in the 
experimental evaluation of a number of transient flowmeters 
(3), including various forms of electromagnetic, angular 
momentum, hotwire and pressure drop devices (e.g., see (4) for 
discussion). None of these schemes proved suitable, and a 
theoretical calculation was used to determine instantaneous 
flow rate from existing instantaneous pressure drop methods 
(1,2). This paper describes an experimental technique used 
to check the validity of the theoretical calculation. 
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Apparatus 

The experimental setup, shown schematically in Fig. 1, 
consisted of a pressurized water tank supplying a rocket injec- 
tor through a cavitating venturi and a piston-type flow- 
modulating device. The water jet from the injector (actually 
a converging group of 12 individual jets) impinged on a T- 
shaped collector tube (see Fig. 1) so that the total axial 
momentum of the jet appeared as a thrust load on the collec- 
tor. The collector was mounted on a cantilever beam to which 
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Fig. 3 Pressure drop-flow rate amplitude ratio vs. frequency 
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four high sensitivity Baldwin C-7 strain gages were attached. 
Thus the strain gage signal was directly proportional to the 
instantaneous momentum of the jet. The measured natural 
frequency of the strain beam was 700 cps and hence was ace- 
quate for the 60 to 140 eps range of forcing frequencies used. 
The injector pressure drop was measured by a Li-Liu trans- 
ducer located in the injector manifold (the injector exhausted 
to atmosphere; hence the pressure drop and-manifold pres- 
sure were identical). 

The tests consisted of varying the small amplitude sinus- 
oidal forcing frequency produced by the flow modulator, and 
measuring the amplitude ratio and phase difference between 
the resulting flow and pressure drop signal oscillations. The 
fixed time displacement required for the steady-state flow 
component to traverse the distance between the injector orifice 
and the collector was obtained from turbine flowmeter and 
jet diameter measurements. Details of the apparatus ce- 
velopment, final configuration, procedures, calibration tech- 
niques, etc., are included in detail in (5). 


Discussion 


It has been shown in (1,4, etc.) that for small sinusoidal 
oscillations of an incompressible fluid, the instantaneous frac- 
tional pressure drop y and the instantaneous fractional flow 
rate uw are related by 


v 2 ia 
cos @ 
where W and yu are defined by 
Ap = Ap(1 + 
tm = m(1 + 


The phase lag a@ of flow behind pressure drop results from the 
fluid inertia, and hence is a function of orifice dimensions, 
steady-state fluid velocities, and fluid density. Values of a, 
as computed theoretically (2) for the specific orifices used in 
the present experimental study, appear as the curve of Fig 2. 
The experimental values of a obtained from the apparatus of 
Fig. 1 agree quite well with the theoretical predictions of (1) 
and (2). Comparison of the theoretical and experimental 
amplitude ratios, shown in Fig. 3, indicate that the experi- 
mental values are somewhat low. However, it must be re- 
membered that the amplitude might be expected to be rather 
sensitive to such indeterminate but consistent factors as 
location of the pressure pickup with respect to a true stagna- 
tion region, effect of oscillations on the velocity distribution in 
the jet and on entrainment of surrounding air, ete. 


Conclusions 


Measurement of transient liquid flow rates has been ac- 
complished by the use of a momentum sensitive “flowmeter” 
at frequencies up to 140 eps. 

The validity of the theoretical flow-pressure relationship 
for simple orifices (1,2) has been established by direct meas- 
urement of instantaneous flow rate and pressure for a specific 
configuration. Agreement between experimental and theoreti- 
cal phase lags was excellent, but only fair correlation with pre- 
dicted amplitudes was obtained. 
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Similitude Treatment of Hypersonic 
Stagnation Heat Transfer' 


DANIEL E. ROSNER? 


AeroChem Research Laboratories, Inc., Princeton, N. J. 


Using the “similarity”? between heat and mass transport 
an expression for the convective heat transfer rate at the 
nose of a blunt body is developed for dissociated but 
“chemically frozen”? laminar flow with arbitrary surface 
catalytic activity. The result is shown to be identical in 
structure with recent solutions to the boundary layer 
equations. 


HE LAMINAR boundary layer equations for the case of 

hypersonic stagnation point flow with no gas phase atom 
recombination have been integrated by Lees (1)* and Fay 
and Riddell (2) for a “fully” catalytic surface, and by Goulard 
(3) and Seala (4) for surfaces of arbitrary catalytic activity. 
As an example, with the assumptions pu = constant, h. > hu, 
Goulard gives the heat transfer rate as 


j = 0.664(Bpeue)' */*he {1 + — 1]a.Q/he} [1] 
where the factor 
¢ = [1 + [2] 


is a correction for the fact that the steady-state atom concen- 
tration a» at the wall need not, for all combinations of cata- 
lytic activity and flight conditions, be small compared to the 
value a, at the outer edge of the boundary layer. The pur- 
pose of this note is to indicate that a result of identical struc- 
ture may easily be obtained using an elementary chemical en- 
gineering approach to this problem. This simplified approach 
suggests applications to related problems (5) in the calcula- 
tion of energy transfer rates in partially dissociated gases, 
using existing “low-speed” and “undissociated”’ heat transfer 
data. 


Analysis 


If we write the contribution due to ordinary conduction as 
= (Str)PaoU o(he = hu) [3] 


and the contribution due to atom diffusion to the wall as the 
product of the heat of recombination and the mass transfer 
rate 


Go = Q(Stp)PaU — Aw) [4] 


then the net heat transfer g will be the sum of these two 
rates. From similitude theory, we assume that the Stanton 
number Sty for mass transport is obtainable from the Stanton 
number St, for heat transport by making the replacement: 
Pr, — Pr, (Prp is commonly referred to as the Schmidt num- 
ber). This is asymptotically exact for the low-speed, con- 
stant property case when the free stream concentration of re- 
actant is small compared to unity. 

For a laminar flow at the blunt nose of a body of revolution 


St, can be approximated, for example, by Sibulkin’s incom- 


pressible fluid formula (6), applied to the stagnation region 
behind the normal shock 


St, = 0.763(p._U [5] 
If the subscripts X are replaced by D, we obtain the cor- 
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responding Stanton number for mass transport, which is then 
introduced into the diffusive contribution gp to the net energy 
transfer. Furthermore, if the surface recombination kinetics 
are described by a first order rate law (3) of the form 


Re = ky pws [6) 


then, in the steady state, the conservation equation for atoms 
at the gas-solid interface may be written 


kwPww = (ae — (7) 


This relation may be considered to determine the “eigen- 
value” of the atom concentration at the wall. If this value for 
Q@y» is introduced into the Equation [4] and use is made of the 
fact that the enthalpy h of the partially dissociated gas is 
comprised everywhere of the sum h + a@Q, then the net heat 
transfer rate ¢ = g, + dp becomes 


G = — hw){1 + — 1] X 


where ¢ is again a correction factor of the form @/(1+ @); @ 
being the relevant ‘catalytic parameter” (7) 


= [9] 


This result should be compared to Goulard’s result (Eq. [1]) 
the latter being based upon the assumption of a highly cooled 
wall (h. >> hw). Itis observed that when the complete driving 
force (he — hw») is retained in place of just h., the catalytic 
correction factor g should modify the entire term (Le®-* — 1) 
and not just the Lewis number term. This can be veri- 
fied by checking the reasonableness of the limiting case g > 0, 
since g must approach q in the absence of surface catalysis.‘ 

By way of further comparison, the “fully” catalytic wall 
investigated by Fay and Riddell (2) corresponds to the limit- 
ing case g > 1. It is interesting then to write down the sug- 
gested correlation equation which best represented their 
numerical integrations 


{1 + — (10) 


Upon setting ¢ in Equation [8] equal to unity, and comparing 
the result with Equation [10] it will be noted that, apart from 
the correction factor for variable pu, these equations are 
virtually identical. Since detailed calculations indicate that, 
at most, pu at the wall is five times the value at the outer edge 
of the boundary layer, the fact that this ratio is raised to the 
1/,9 power in Equation [10] makes the corresponding correction 
of the order of only 5 per cent. 

It appears then that relatively simple techniques of the 
type adopted here, guided by some “exact” numerical results 
for extreme cases, can provide a realistic approach to related 
problems. Conversely, these techniques should facilitate the 
“correlation” of machine calculations. 


Nomenclature 

e€ = catalytic parameter (Eq. [9]) 

h = sensible or “frozen’’ enthalpy of the mixture 

h = enthalpy of the mixture 

ky = rate constant for surface atom recombination 

Le = “frozen’’ Lewis number (Le = Pr,/Prp), assumed con- 
stant 

Pr = “frozen’’ Prandtl number, assumed constant 

gq = heat transfer rate 

Q = heat of recombination, assumed constant 

Ry» = mass rate of atom recombination at surface 

St = Stanton number 


‘Note that in this limiting process (he — hw) — (he — hw) 


SINCE aw —> ay. 
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aQ/(he — hw)} 18} 
j 
4 
| 


U_ = fluid velocity relative to body 

a = mass fraction of atoms in mixture 
8 = inviscid velocity gradient at the nose 
p = density of fluid 

u = absolute viscosity of fluid 
Subscripts 

2 = before normal shock 

e = atouter edge of the boundary layer 
w = atthe wall (surface of body) 

= pertaining to thermal conduction 
D = pertaining to concentration diffusion 
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Motion Units to Simplify Space Travel 
Computations 
BEMROSE BOYD! 
General Motors Corp., Milwaukee, Wisc. 


NDER the special conditions of constant acceleration and 
zero initial velocity, we have 
V =at {1] 


This is a basic mathematical relationship and is valid in any 
consistent set of units. It can be nondimensionalized by use 
of a reference velocity V, 


[2] 


Furthermore, the quantities can be expressed in arbitrary 
units if proper conversion factors are incorporated 


V’ = [3a] 
a’ = Ka [3b] 
t’ = Kst [3c] 
= a’t Ka't [4] 


If V, is chosen to be the velocity of light, and a’ is measured 
in units of Earth gravity and ¢’ in years, the composite co- 
efficient K has the numerical value 1.03. Thus, within an 
accuracy of 3 per cent, the original equation (Eq. [1]) is valid 
in these units with V replaced by the ratio V’/V, 

[5] 
The ratio V’/V, might be given a designation such as “lux 
number” in analogy to Mach number. 
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Equation [5] shows that, except for relativistic effects, an 
acceleration of 1 g applied continuously for one year would 
bring a body up to the velocity of light. Similarly an ac- 
celeration of 10~* g for one year or 1 g for 10~* year (8.8 hr) 
would produce a terminal velocity of this fraction (10~%) of 
the velocity of light, or 670,000 mph. Considering the 
velocity of light as a fundamental limitation on a physical 
system, the above relationship facilitates an understanding of 
values of thrust and travel time required for various space 
travel missions. 

These units of generalized velocity, acceleration and tine 
which lead to the simplified form of Equation [5] are par- 
ticularly useful in conjunction with distances expressed in 
light-years as is done in some areas of astronomical usage. 
The mean distance of the Earth from the sun is 1.58 X 105 
light-years (8.3 light-min), and the distance to Pluto, tiie 
farthest planet, is 6.26 XX 10~‘ light-years (5.6 light-hr), 
whereas the distance to Proxima Centauri, one of the nearest 
stars, is 4.3 light-years. 

It is evident that there is a big step between the minor 
planets (Mercury, Venus, Earth, Mars) and the major 
planets (Jupiter, Saturn, Uranus, Neptune, Pluto), and a 
much greater step between the limits of the solar system and 
any stellar or extra-solar body. All space travel is presently 
acceleration-limited by virtue of thrust limitation. To the 
extent that this limitation is removed by developments in 
nuclear propulsion or utilization of radiation, travel outside 
the solar system will be velocity-limited. Travel within the 
solar system might continue to be acceleration-limited as 
a result of practical considerations of human tolerances, 
economics, etc. 


A Short Form Method for Determining 
Near-Cireular Orbit Quantities 


JOSEPH N. BENEZRA! 
Boeing Airplane Co., Seattle, Wash. 


A method is presented that allows quick and accurate de- 
termination of the precision with which a satellite must be 
injected into a “‘near-circular’’ orbit so that apogee and 
perigee remain within preselected limits. Various orbital 
quantities of interest are presented as functions of injec- 
tion velocity and attitude “‘errors,”’ given the injection al- 
titude. The sensitivities of these orbital quantities to 
velocity or attitude errors are given as functions of the 
errors themselves and the injection altitude. 


IVEN the injection altitude of a satellite, it is possible to 
determine the effect of velocity and attitude injection 
errors on various orbital quantities, such as perigee and apogee 
velocities and altitudes, period, etc. (By “errors” is meant 
simply the differences AV, y between the injection quantities 
of velocity and attitude and those corresponding to a circular 
orbit at the injection altitude.) Singer (1)? has done this 
rigorously for apogee and perigee altitudes. Jensen (2) briefly 
discusses by numerical example the sensitivity of apogee and 
perigee altitudes as well as orbit period to injection errors. 
Both Jensen (2) and Arnowitz (3) deal with the question of 
injection errors which can be tolerated in order to restrain 
the satellite to orbit within specified apogee and perigee limits. 
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The short form presented here is obtained by applying stand- 
ard perturbation techniques to basic equations. It has the 
advantage of expressing the orbit quantities listed below and 
their sensitivity to injection errors as explicit functions of the 
injection errors and injection altitude: 
Eccentricity, semimajor axis, period, energy per unit mass. 
Velocity and altitude at apogee and perigee. 


Method 


An expression for orbit eccentricity in terms of injection 
conditions is (see Fig. 1) 


V2 v2 
€= [1 - (1 - cos? {1] 


Applying standard perturbation technique by letting V = 
V. + AV and cos? y = 1 — y? it can be easily shown that 
Equation [1] becomes 


ELLIPTIC 
ORBIT 


AV\?71'4 
+4 (57) ] [2] SEML-MINOR 
AXIS 
SEMI-MAJOR AXIS 
a Fig. 1 Elliptic orbit characteristics 
Table 1 Expressions for near-circular orbit quantities 
AV ~0,¥ #0 =0,y¥ = 
Quantity General case Circular orbit 
[ AV 
€ w+ v. 0 
b 
(: 1 —7y2)(1 2 : 
R R +2 -1 R 
h h : h 
R 
+R 
Vp AV 
ha 
1+2e+2¢ 1 
Va hp 41 
R 
Vp — Va AV 
2(1+4) 
D 
+1 1 
a AV h 
> +5 (47) ] 
R\'2 \ 
T. + + 
AV AV gk 
23h (51) 
1424 
R 
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Table 2 Sensitivities of near-circular orbit quantities : | 
AV #0, #0 AV #0, #0 
Derivative General case Derivative General case 


From the polar equation of the ellipse and conservation of 


w moment of momentum we obtain for the apogee and perigee 
S altitudes 
2V? cos? 
_ w Again, letting V = V. + AV, etc., Equation [3] becomes 
hé h 
[4] 
ae 400 600 800 1000 1200 1400 1600 
are: If Y = 0, Equation [4] shows immediately that perigee occurs 
APOGEE ALTITUDE, h,, (ST. Mi.) at injection when AV > 0, and that apogee occurs at injection 
when AV <0. It is interesting to note that when AV = 0 
Fig. 2a Effect of injection velocity V and attitude y on apogee 
and perigve altitedes; injection at 300 statute miles altitude y A gama is automatically at one end of the minor axis. 
rgk? 
5 
[1490 1600 _| — (5) 
= 
6°} — 1000 \ and V? = V,? = gR?/r, it follows that r = a, which can occur 
| |__| 800 fat at only two points on the ellipse, i.e., at either end of the 
4° 600 200 minor axis. 
= 250 From Equations [2 and 4], it is possible to plot curves of 
w apogee as a function of perigee altitude showing lines of con- 
stant injection attitude and velocity. Alternately, curves 
ge || to ear | may be plotted of injection attitude as a function of injection 
Ll LG velocity (given the injection altitude) for constant values of 
di 24 _ 2\__) 26 os? 10? apogee and perigee altitude. These plots are shown in Figs. 
. 2a and 2b for an injection altitude of 300 statute miles. From 


the curves, tolerances may be established on injection attitude 
and velocity when apogee and perigee are held within desig- 
Fig. 2b Concluded—injection at 300 statute miles nated limits. 


INJECTION VELOCITY, V, (FPS) 
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Results 


Summarized in Table 1 are quantities pertaining to near- 
circular orbits expressed in dimensionless form as functions 
of injection errors in velocity and attitude. Those involving 
apogee and perigee velocities are derived from Equations 
(3and 4]. The expression for semimajor axis a may be derived 
from [5], whence follows the expressions for orbit period using 
Kepler’s Third Law. It will be noted that the expression for 
orbit period depends only upon the velocity and altitude at 
injection, as would be expected. 

By differentiating the expressions in Table 1 with respect to 
either attitude or velocity error, the sensitivities are obtained; 
these have been tabulated in Table 2. 


Accuracy 


This scheme of determining various orbit quantities should 
be sufficiently accurate for preliminary performance work. 
As the injection errors (or eccentricity) become smaller, 
the accuracy improves; this is apparent from Fig. 2a. It 
may be seen that for small positive velocity injection errors, 
perigee occurs at the injection altitude for zero attitude error. 
However, as the velocity error increases, perigee is seen to 
occur slightly above the injection altitude. Rigorously, of 
course, perigee occurs precisely at injection altitude regard- 
less of the magnitude of the injection velocity (V > V.). A 
quick indication of the relative accuracy of the method may 
be had from the following examples: 

Suppose a satellite is injected with y = 0.07 rad; V = 
26,000 fps; h = 300 statute miles. The circular velocity at 
this altitude V. = 24,972 fps. From Equation [1], « = 0.1092. 
The approximate expression, Equation [2], gives « 0.1081. 
The discrepancy between the two is —1.0 per cent. 

It is instructive to consider another example. The Naval 
Research Laboratory reported that the Army’s first Explorer 
satellite, on the basis of its first 32 runs around the Earth, 
goes out as far as 1587 statute miles and comes in as close as 
219 miles. From these figures the velocities at apogee and 
perigee may be computed using Kepler’s laws; they are 20,282 
fps and 26,938 fps, respectively. The period is 114.818 min. 
To obtain the corresponding figures from the approximate 
method, the eccentricity of the orbit is first computed. 

From Table 1, (ha/R + 1)/(he/R + 1) 14+ 2€4+ 22, 
whence « & 0.1426. Since actual injection data are not 
known, it may be assumed for computation purposes that in- 
jection takes place anywhere between apogee and perigee. 
For simplicity, assume that injection occurs at perigee. The 
circular orbit velocity at this altitude is 25,214 fps. From 
Table 1 


sare 


giving 20,078 fps and 26,756 fps for apogee and perigee ve- 
locity, respectively. Again from Table 1 


3 9 
T= = 115.16 min 


h 3/5 
(*) (: + = 91.419 min 
g R 


the period of a circular orbit at perigee altitude. Comparing 
the results of the approximate and rigorous methods we find: 

discrepancy in V4 = —1.01 per cent 

discrepancy in Vp = —0.68 per cent 

discrepancy in T = +0.30 per cent 

Higher order terms have been retained in the expression 
presented in both Tables 1 and 2 because, in general, they 
enhance the accuracy of the results in case the injection 


where 


3 The eccentricity may be determined directly in this example 
from the relationship « = (ra — rp)/(ra + rp) = 0.1409. Using 
this slightly more accurate value of ¢ the resulting discrepancies 
in Va, Vp and T are even less. 
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errors may not be, mathematically speaking, small. The ex- 
pressions may be further simplified at the discretion of the 
user. 


Nomenclature 

a = semimajor axis 

b = semiminor axis 

e = linear eccentricity [a? — b?]'/2 

E = total energy per unit mass; zero potential energy ref- 


erence at Earth’s surface 


F, F’ = foci of ellipse 

g = acceleration of gravity at Earth’s surface, 32.174 fps? 
at 45 deg latitude 

h = injection altitude, orbit altitude 

r = injection radius 

R = Earth’s radius, 20.855 < 10° ft at 45 deg latitude 

T = orbit period 

V = injection velocity, orbital velocity 

AV = velocity excess or defect, V — V, 

€ = eccentricity, [1 — (b/a)?]'/2 

0 = geocentric angle 

y = injection attitude, orbit direction angle (angle formed 
by tangent line to orbit and local horizontal) 

Subscripts 

A = apogee 

c = circular orbit corresponding to injection altitude 

= perigee 

apogee, perigee 
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General Equation for Rocket Velocity 
SYLVAN RUBIN! 


Stanford Research Institute, Menlo Park, Calif. 


T HAS been suggested that high rocket velocities can be 
attained by means of propulsion systems producing very 
high exhaust velocities, such as ion propulsion. In comparison 
with conventional chemical rocket motors, such systems in- 
troduce two new independent variables: 

1 The exhaust velocity may be varied, independently of 
the specific energy of the fuel. 

2 The fuel and the propellant may be different materials, 
with different masses and rates of consumption. 

An equation for the final velocity of a rocket using any pro- 
pulsion system is derived, in terms of the exhaust velocity, the 
specific energy of the fuel, and the mass ratios of fuel and pro- 
pellant to the rocket mass at burnout. The calculation is 
made for a single stage rocket. 

It is shown that the final rocket velocity is a maximum 
when the exhaust velocity is equal to the fuel kinetic velocity. 
Since the fuel kinetic velocity is approximately equal to the 
exhaust velocity of a conventional rocket motor using the 
same fuel, the conventional motor yields the maximum 
rocket velocity for a given fuel. In other words, the specific 
impulse of the propellant should be matched to the specific 
impulse of the fuel. 

For exhaust velocities greater than the fuel kinetic velocity, 
the final rocket velocity is inversely proportional to the exhaust 
velocity. 
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Assumptions 
1 The exhaust velocity is freely variable, independent of 
the fuel properties, but constant throughout the propulsion 


phase. 

2 All the fuel and propellant to be used are initially 
loaded in the rocket, and they are consumed in a constant 
proportion, equal to their initial loading ratio, so that both 
are completely consumed at the end of the propulsion phase. 
It should be noted that all or part of the used fuel may be 
used for propellant. Therefore, the propellant mass defined 
in the equations is only the initially inert portion of the propul- 
sive material. 

3 Used fuel, if not expelled as propellant, is jettisoned 
continuously during the propulsion phase. 

4 The exhaust jet is unidirectional, with no momentum 
going into lateral or divergent flow velocity components. 

5 Gravitational and drag forces are neglected. 


Calculation 


The equations for the differential momentum balance and 
for the conversion of fuel energy to kinetic energy are most 
readily stated in a coordinate system momentarily stationary 
with respect to the rocket. 


Momentum: —cdp =d[((M+p+f)V]) =(M+p4+f)dV [1] 


1 
Energy: —dE = edf = 3° dp [2] 


Combining the above equations 
2edf 
(M+ptf) 


Since fuel and propellant are used in a constant ratio 
p+f=af 
where a =1+(P/F) 
2e f 0 df 2e, M+aF 
c JF (M+af) ca M 


-" (75 In(1 + R, + Ry) 


dV = [3] 


[4] 


It is evident that if no further restrictions are applied, the 
burnout velocity will generally be inversely proportional to the 
exhaust velocity, and proportional to the square of the fuel 
kinetic velocity. However, for most rockets of practical in- 
terest, there will be further restrictions, so that this general 
equation must be considered to give an upper limit to the 
possible burnout velocity. 

The simplest restriction is obtained for a conventional 
chemical rocket, where the fuel and propellant are identical. 
In this case c = »,, and a = 1, giving the conventional equa- 
tion for burnout velocity. 

A very important restriction applies to all cases in which 
c < v,, arising from the necessity for supplying enough inert 
propellant mass to carry off all the energy as kinetic energy. 
The minimum additional inert propellant P is given by the 
equation for energy balance 


1 1 
= (F + P)ec? 


On substituting this restriction into Equation [4] 
V, =cln(1 +R, + R;) [5] 


which is again identical with the equation for a chemical 
rocket. This gives a smaller upper limit to burnout velocity 
than Equation [4] for c < »,. 

These two limiting curves for burnout velocity correspond to 
the two conditions of: (a) limited total energy and more than 
sufficient propellant mass; (b) limited propellant mass and 
more than sufficient energy. Because of the way », is defined, 


e¢--, 


Fig. 1 Upper limits to rocket burnout velocity V,.as a function 

of varying exhaust velocity c for a fixed fuel kinetic velocity v,. 

For any value of c, the lowest limit is the actual limit, indicated 
by the solid line 


the transition point between these two conditions occurs at 
c = v,, where the two limiting curves intersect. Jf additional 
mass or additional energy could be obtained in flight, it would 
be possible to continue on the ascending branch of the cor- 
responding curve to a higher burnout velocity, but without 
such a process, the maximum possible burnout velocity occurs 
atc = 


Conclusions 


It is evident that the greatest final velocity will be given to 
a rocket if the specific impulse of the propellant is equal to the 
specific impulse of the fuel. 

The physical reasons for this behavior are quite obvious, 
after one has the equations. In the high exhaust velocity re- 
gion, increasing exhaust velocity causes a greater fraction of 
the available energy to be expended on the propellant, and 
correspondingly less supplied to the rocket. This is also 
pointed out by Oberth? showing that impulse is inversely pro- 
portional to particle velocity for a given energy input. In the 
low velocity region, reduction of exhaust velocity is ac- 
complished essentially by diluting the fuel with inert propel- 
lant, which obviously can only diminish the rocket velocity. 

Burning a chemical fuel in a conventional rocket motor is 
the most efficient propulsion system for obtaining maximum 
rocket velocity from that particular fuel. Ion propulsion 
systems are useful for utilization of high energy fuels which 
cannot be directly burned, such as nuclear fuels. The ion 
velocity must equal the fuel kinetic velocity to obtain maxi- 
mum effectiveness. For example, the kinetic velocity corre- 
sponding to full utilization of the fission energy of pure U?* is 
1.2 X 10° cm per sec (about 0.7 Mev per nucleon). 

These results do not apply, of course, to space vehicles de- 
signed to obtain energy from their environment during pro- 
pulsion, as by means of solar batteries or solar radiation pres- 
sure. 


Nomenclature 

c = exhaust velocity 

V_ = rocket velocity (variable) 

V, = final rocket velocity 

e = specific energy of fuel (available work per unit mass) 
v, = kinetic velocity of fuel, defined by e = }2,? 
E = energy supplied by fuel (variable) 

M = rocket mass (less fuel and propellant) 

f = remaining fuel mass (variable) 

F = total fuel mass 

p = remaining propellant mass (variable) 

P = total propellant mass 

Ry; = F/M (fuel-mass ratio) 

R, = P/M (propellant-mass ratio) 


2Oberth, H., ‘“Wege Zur Raumschiffahrt,” Oldenbourg, 
Munich, 3rd ed., 1929, p. 410, Equation 237. 
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Technical Comments 


Reply to Crocco’s Criticism of the 
Zucrow-Osborn Paper 


M. J. ZUCROW! and J. R. OSBORN? 
Purdue University, Lafayette, Ind. 


N THE December 1958 issue of Jer PRopuLsIoN there ap- 
peared a Technical Note hereafter designated as Refer- 
ence (1).4 The authors of (1) presented the results they ob- 
tained from an experimental study of combustion pressure 
oscillations of the longitudinal mode. They conducted their 
experiments with an approximately 250-lb thrust liquid 
propellant rocket motor burning 95 per cent ethyl! alcohol and 
liquid oxygen. The authors of this discussion were greatly 
interested in the results presented in (1) since they had per- 
formed similar experiments with a rocket motor burning pre- 
mixed gases; their experiments were reported in (2). 

We were pleased to note that, in the ranges where the re- 
sults obtained by (1) could be compared with our own experi- 
mental results, the agreement was satisfactory. For the 
guidance of the reader, it is pointed out that most of the 
results published in (2) were obtained approximately two years 
ago and were reported at three different symposia (3,4,5) on 
liquid rocket combustion stability, and that those symposia 
were attended by one or more of the authors of (1). 

In the section “Discussion” the authors of (1) make the 
following statement with regard to the effect of combustion 
chamber length on combustion pressure oscillations. We 
quote, ‘‘Note this existence of an upper limit clearly invalidates 
oscillation-producing mechanism advanced by Zucrow and 
Osborn (2) as discussed by Crocco (6) since this mechanism 
cannot account for the observed cessation of oscillations with 
increasing chamber lengths.” We will show later in this dis- 
cussion that not only is that statement incorrect, but the 
possible existence of such an upper critical combustion chamber 
length was postulated by Osborn and Schiewe some time ago 

We now turn our attention to a Technical Comment which 
appeared in the same issue of Jer Proputsion. That article 
(6), by Prof. L. Crocco of Princeton University, is entitled 
“Comments on the Zucrow-Osborn Paper on Combustion Os- 
cillations” (2). 

In a Comment of more than 1500 words by Crocco upon a 
paper reporting experimental data and results, he did not 
make a single comment concerning the experiments them- 
selves. In his first paragraph he makes the assertion con- 
cerning his own theories that ‘in the United States and abroad 
_. . the mechanism (time lag theory) on which these theories 
are based have been accepted as correct.”’ To the senior 
author this is a surprising statement. The senior author and 
other investigators, to cite a few (8,9,10,11,12,13), have 
challenged those theories at practically every symposium held 
on combustion stability. In fact, the primary reason for 
initiating a research program on combustion pressure oscilla- 
tions at Purdue University can be attributed to the conviction 
of the senior author of this discussion that the time lag theory 
was insufficient for explaining the complex phenomena asso- 
ciated with high-frequency combustion pressure oscillations. 

In the second paragraph of his comments Crocco implies 
in a somewhat disparaging tone that the experimental tech- 
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niques used at Purdue are similar to those used at Princeton. 
He then correctly points out that the Purdue experiments were 
conducted with rocket motors burning premixed gases. 
We are of the opinion, and our progress indicates, that from 
the standpoints of obtaining basic information (1), simplifying 
the experimental procedures, increasing the speed of conducting 
the experiments, and economy, there are a great many ad- 
vantages in using premixed gaseous rocket motors. 

In the third paragraph of the aforementioned Comment, 
Crocco presents his summary of the conclusions to be reached 
from our experiments. In that summary he interprets our 
conclusions under five numbered items. Since there is no 
disagreement with Item 1 we will discuss only the remaining 
four. 

Item 2. He states that we conclude that “The tendency 
toward instability . . . , and therefore (it) increases with de- 
creasing frequency.” We made no such statement in our 
paper (2). As a matter of fact, the statement is inconsistent 
with our data on the effect of combustion chamber length. 
Apparently Prof. Crocco misunderstood our experimental 
objectives and also our physical explanation of the existence 
of a (lower) critical combustion chamber length, that is, the 
length of combustion chamber below which no combustion 
pressure oscillations were detected. Our experiments on the 
effect of combustion chamber length were conducted for the 
purpose of determining the (lower) critical combustion cham- 
ber length. Consequently, our discussion and comments are 
applicable only to the range of lengths we investigated. We 
refer Prof. Crocco to our conclusion ((2), p. 658). ““Experi- 
menters using liquid bipropellant rocket motors have also 
noted that long combustion chambers tended to initiate com- 
bustion pressure oscillations of the longitudinal mode more 
readily than do short ones” (14). The above statement is 
correct and refers to rocket motors of practical lengths and hav- 
ing the lengths we investigated. For Crocco’s comment to be 
consistent with our data and physical explanation it should 
read “‘and therefore (the tendency toward instability) increases 
with decreasing frequency for combustion chamber lengths 
near the (lower) critical length.” Any other conclusion is in- 
valid and incomplete. The authors are also of the opinion 
that the frequency of the combustion pressure oscillations 
should not be considered to be an independent variable. 

Ttem 3. This item, which is also concerned with the effect 
of combustion chamber length, is again an incorrect conclusion 
drawn by Crocco from our experiments, and not by us. His 
statement, “Oscillations with a pressure amplitude that 
steadily increases with increasing chamber length,” as in the 
case of Item 2 is inconsistent with our data and physical ex- 
planation. Table 1, reproduced from (2), summarizes the re- 
sults of our investigation of the effects of combustion chamber 
length. 


Table 1 Typical reduced data from the study on the effect 
of chamber length for superheated propane and air 


Frequency Amplitude 


Chamber Type of of oscilla- of oscilla- 

length oscillation tions, cps tions, psi 
28 shock (L)* 570 22 
22.3 shock (L) 685 22 
16.3 shock (L) 948 17 
12.6 sinusoidal (L) 1219 4 


Values for the variables: Equivalence ratio, 1.2 + 0.08; | 
propane line pressure drop, 158 + 8 psi; air line pressure 
drop, 150 + 10 psi; mean combustion pressure, 45 psia + 
2 psi. | 

*L indicates a longitudinal mode of oscillation. 


221 


~ 
— 
| 
ion 
UK. 
ted 
at 
nal 
ild 
ut 
Irs 
j 
he 
IS, 
of 
ad 
SO 
he 
y. 
Is 
+h 
n 
. 
IS 
ye 
| 
| 
| 
| 
2 
| 
= 


The table shows that the pressure amplitudes of the oscilla- 
tions were identical for a 22- and 28-in. long combustion cham- 
ber, and that the rapid change in the amplitude of the oscilla- 
tions occurred near the critical length. The error arose from 
the same considerations we have discussed under Item 2 
above. As pointed out earlier in this discussion, we have for 
some time been aware of the possible existence of an (upper) 
critical combustion chamber length, that is, a length beyond 
which there would be no combustion pressure oscillations. 
What is noteworthy, however, is that our former graduate 
student R. M. Schiewe arrived at the above conclusion by 
applying the very mechanism which Prof. Crocco asserts can- 
not explain the existence of the upper critical combustion 
length. We quote from a report by Osborn and Schiewe (7): 
“Tt is conceivable, however, that if the combustion chamber 
were made sufficiently long the pressure disturbance will be 
considerably attenuated before it can traverse twice the 
length of the combustion chamber. Thus the returning pres- 
sure wave could be too weak for raising the temperature and 
pressure of the unburned and burning mixture significantly, 
and, therefore, there would be no amplification of the pressure 
disturbance.” Because the upper critical combustion cham- 
ber length was considerably longer than the lengths of practical 
rocket motors, we did not deem it of sufficient importance to 
determine it experimentally. We are glad to learn, however, 
that the authors of (1) have determined its existence experi- 
mentally. 

Item 4. This is concerned with our prognostication that the 
mechanism by which the transverse modes of combustion 
pressure oscillation arise is probably similar to that influencing 
the longitudinal modes. Prof. Crocco presents a long discus- 
sion which completely ignores the discussions presented in 
our paper (2). 

Item 5. We used the phrase “a type of detonation wave” 
to make a distinction between a Chapman-Jouguet wave and 
a shock wave with combustion taking place behind the wave. 
We are still unconvinced that our choice of words is in- 
correct (15). 

In the closing paragraph of his comments (6) Crocco implies 
that we have no right to advance a new mechanism for ex- 
plaining combustion instability in rockets. He also implies 
that we have not made even “a limited perusal” of the 
published literature. Such remarks are irrelevant and in- 
accurate. Furthermore, they contribute nothing to improv- 
ing our understanding of the complicated phenomena asso- 
ciated with combustion pressure oscillations. Reference (16) 
of this discussion, listed as (6) in our paper (2), shows that we 
have made a rather complete study of the literature. 

We sincerely regret that we were not given the opportunity 
of reading Prof. Crocco’s comments before they appeared in 
the December 1958 issue of Jer PRopuLsion. Consequently, 
we were denied the opportunity of publishing our reply in the 
same issue. 
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And Now the Editor Has His Say 


IKE the unwary referee in a TV wrestling match, your 
editor was jostled by one of the combatants in the 
Crocco-Zucrow fray. Prof. Zucrow in the above Comment 
questions our editorial policy, implying that we should 
have shown him Prof. Crocco’s Technical Comment (JET 
PRroPpuULSION, December 1958) before we published it. So 
now we have to explain our policy and defend our decision. 
We consider it to be an established fact that intellectual 
controversy is a strong stimulus to scientific progress. The 
history of science is full of such examples. -The argument 
over the nature of the primary cosmic rays in the 1930’s 
spurred many physicists to perform searching experiments 
that ultimately added much to our knowledge, and the 
disputes at meetings of the American Physical Society at 
the time often became quite heated. No one would have 
suggested that these disputes ought to have been settled 
“behind the scenes” or that the disputants should have 
cleared their arguments with each other beforehand. 

The position of your editor is that publishing such 
Comments in the ARS Journat when they are not de- 
famatory helps to clarify the issues for the benefit of our 
readers, and in our view, the JouRNAL exists for the 
benefit of its readers. We therefore published Crocco’s 
criticism of Zucrow’s paper and also Zucrow’s defense 
against Crocco’s criticism. We earnestly recommend to 
those of our readers interested in rocket combustion that 
they read both Technical Comments: We ourselves learned 
a lot about the scientific issues at stake, and we are grate- 
ful to both authors for it. 

—Martin Summerfield 
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This is the cabin interior of Lockheed’s test Electra airplane. The 
engineer adjusts control panel for oscillographic recorder at lower right. 


“WE USE DUPONT LINO-WRIT 4 
FOR OUR MOST CRITICAL TESTS” 


—Dynamics Flight Test Laboratory, 
Lockheed Aircraft Corporation, Burbank, California 


Here’s a report from a Lockheed Group Engineer. 
“The photorecording papers we use have to deliver 
dependable performance while combining a high 
enough degree of emulsion sensitivity to handle 
special tests. We started using Lino-Writ 4 early in 
1957 and now use about 2000 feet per week. 


“The first advantage we noticed was Lino-Writ’s 
extra thinness which allows us to load almost 100% 
more footage in a 12-inch magazine than we could 
before. This meant that we could run longer tests 
without interruption and obtain a correspondingly 
larger volume of test data from each flight.” 


Other advantages of Lino-Writ 4 cited by Lock- 


REG. us pat OFF 


Better Things for Better Living . . . through Chemistry 


Marcu 1959 


* 


heed: it can be used as an intermediate in making 
diazo reproductions; it is strong and durable, capable 
of withstanding stresses incurred in processing; no 
quality problems with either exposed or unexposed 
rolls because of high temperature and humidity 
storage conditions—even the oldest records are still 


clear and thoroughly legible. 


Isn’t it time you investigated the Lino- Writ line 
for your oscillographic installation ? 


E. I. du Pont de Nemours & Co. (Inc.), Photo 
Products Department, Wilmington 98, Delaware. In 
Canada: Du Pont of Canada Limited, Toronto. 
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New Patents 


Self vaporization of liquid fuels by partial 
oxidation (2,864,233). J. G. Tschinkel, 
Huntsville, Ala., assignor to the U. 8. 
Army. 

Method of vaporizing a normally stable 
liquid fuel taken from the class consisting 
of isopropyl! alcohol and xylene. The 
process comprised vaporizing a mixture 
of the fuel and phenylhydrazine by react- 
ing the mixture with a peroxide compound. 
Buoyed rate gyro (2,865,206). T. R. 
Quermann, Huntington Station, N. Y., 
assignor to United Aircraft Corp. (Ketay 
Dept., Norden Div.) 

Unheated floated rate gyro for control- 
ling servomechanisms and _ stabilization 
systems in missiles. Stability is main- 
tained by means of a bellows, slotted 
spring steel damping cylinder and a fixed 
cam surface. The gyro is a component of 
Atlas, Polaris, Titan and Snark missiles. 
Rocket motor liquid propellant (2,865,727). 
J. W. DeDapper, Los Angeles, Calif., 
meer to North American Aviation, 

ne. 

Fuel comprised of not less than 80 per 
cent (by weight) turpentine, and the re- 
mainder an ethereal liquid in which up to 
4.5 per cent lithium aluminum hydride is 
dissolved, and which is mixed with turpen- 
tine. 

Flame holding device (2,866,314). A. 
Korti, N. Woodbury, Conn., assignor to 
United Aircraft Corp. 

Flameholder mounted in the combustion 
chamber, movable from an external posi- 
tion between a flame holding position and 
a streamlined position. 

Variable exhaust nozzle actuator (2,866,- 
315). K. K. Schakel, Cincinnati, Ohio, 
assignor to Kett Corp. 

Hydraulic device for a gas turbine 

engine afterburner. The exhaust nozzle is 
variable within a first range for normal 
operation, and movable through a second 
range to fully open condition for after- 
burner operation. 
Thrust reversing and silencing (2,866,- 
316). G. B. Towle, R. L. Vasquez and 
A. E. Wetherbee Jr., Newington, Conn., 
assignors to United Aircraft Corp. 

Exhaust walls including small nozzles 
spaced for directing gases aft. Segmental 
plates about the nozzle axis, pivoted at 
their downstream ends, swing outward to 
reverse part of the gases from the nozzles. 
Fuel turbidimeter (2,866,378). H. D. 
Warshaw and E. Nerheim, Minneapolis, 
Minn., assignors to the U.S. Navy. 

Device for determining turbidity of a 
fluid by measuring scattered light as a 
function of turbidity. The intensity of 
light, introduced through a quartz rod 
and the fluid, produces a reading on an 
ammeter. 
engine starter (2,866,385). 
K. W. iller, Granada Hills, Calif., 
assignor to Northrop Aircraft, Inc. 

Ground equipment with multiple-posi- 
tion stepping switches connected to indi- 
vidual control switches on the aircraft. 
Sequence stepping means are operable 
when a given desired condition is reached, 
and switches are turned off to a static con- 
dition as before starting the engine. 


Eprtor’s Nore: Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of the 
Heat Transfer Laboratory, University of 
Minnesota, are gratefully acknowledged. 
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George F. McLaughlin, Contributor 


Vertical takeoff aircraft with jet driven 
rotor (2,866,608). L. H. Leonard, Ingle- 
wood, Calif. 

Jet engines installed at the tips of pro- 
peller blades mounted for rotation about 


the longitudinal axis of the fuselage. 
Blades and engines may be folded into 1 
position parallel with each other by means 
of a pivotal connection between each 
engine and its associated blade. 


Data computer for converting oscillograph 
traces into direct-reading records (2,866,- 
596). E.G. Hoefs and G. T. Litton, Dal- 
las, Texas, assignor to Chance Vought 
Aircraft, Ine. 

Potentiometer mounted on a carriage 
movable longitudinally over an oscillo- 
gram, and a stylus on the carriage trans- 
versely movable for following the trace. 
A computer connected between a movable 
contact connected to the stylus is respon- 
sive to movement of the stylus. 


Aircraft jet thrust control (2,866,610). 
J. E. M. Taylor, Cleveland, Ohio. 
Deflector member located outside, and 
aft of, a jet engine to intercept the jet 
stream. The angle of impact may be 
varied between zero and 90 deg to effect 


variations in the proportion of reverse 
thrust and lift. 
Mechanical stick force producer (2,966,- 
611). H. B. Thompson, Redondo Beach, 
Calif., assignor to Northrop Aircraft, Inc. 
Pilot’s attitude control including means 
for producing a differential ratio between 
ram air pressure and static pressure on the 
airplane in flight. The resulting force 
tends to center the control in neutral posi- 
tion. 
Automatic steering system (2,866,930). 
C. Russell, Milldale, Conn., assignor to 
Bendix Aviation Corp. 
Steering means, course setting means 
and controls for automatically maintaining 
a dirigible craft on a preset course. Fol- 


(continued on page 227) 


Target indicating mechanism (2,866,967). 

L. M. Bernbaum, Philadelphia, Pa. 
Fixed target detection system for per- 

forming sequential search of successive 


sectors in discrete steps. Means respon- 
sive to target signal input to produce a 
vertical deflection on a spaced line ob- 
served on a visual indicator. 
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TEST 


Testing is a vital part of every stage in the development of 
missile and space programs at Lockheed Missiles and Space 
Division. 

The Division maintains one of the most completely 
equipped missile and space test laboratories in the world. 
Equipment includes: altitude, temperature and humidity 
chambers; shaker and vibration systems; G-accelerators; and 
apparatus capable of performing chemical, metallurgical, 
plastic, heat transfer, hydraulic, pneumatic, shock, accelera- 
tion, sinusoidal and random vibration, structural, electrical, 
and electronic tests. Static field testing; research and develop- 
ment testing on controls; testing in ordnance and hydraulics 
and high-pressure gas and propulsion systems are conducted 
at the 4,000 acre, company-owned test base in the Ben 
Lomond mountains near Santa Cruz, California. 

As weapon systems manager for such major, long-term 
projects as the Navy Polaris IRBM; Discoverer Satellite; 
Army Kingfisher; Air Force Q-5 and X-7; and other impor- 
tant research and development programs, Lockheed is 
engaged in expanding the frontiers of technology in all areas. 


(left) Pop-up test of Navy Polaris IRBM. 


(below left) One of the Santa Cruz test stands with 
dynamic thrust mount to simulate flight environment. 
Vibration oscillator functions during static firings. 


(below right) Large centrifuge for environmental 
testing has unique shaker attachment to provide 
vibration simultaneously with high G-loadings. 


Expanding the Frontiers of Space Technology 


Flight testing is conducted at Cape Canaveral, Florida; 
Alamogordo, New Mexico; and Vandenberg AFB near Santa 
Maria, California; in a unique manner. All components and 
sub-systems of a new project are initially tested on known- 
performance, production missiles. Thus, when the final sys- 
tem is ready for first flight, its individual components already 
have flight-tested reliability. This new concept of flight test- 
ing is a major contribution and has enabled Lockheed to 
produce extremely complex missile systems in record time 
and at greatly reduced expense. 

Underwater launch tests—including studies of cavitation. 
wave simulation and skip motion—are carried on at the 
Sunnyvale facility and at the Navy test base on San Clemente 
Island. In addition, structural and other tests are performed 
at Hunter’s Point Naval Shipyard, California. 

If you are experienced in any of the various phases of 
testing, we invite your inquiry. Positions also are available 
in physics, mathematics, chemistry, or one of the engineer- 
ing sciences. Write : Research and Development Staff, Dept. 
C-25, 962 W. El Camino Real, Sunnyvale, California. 


“The organization that contributed most in the past year to the advancement of the 
art of missiles and astronautics.’” NATIONAL MISSILE INDUSTRY CONFERENCE AWARD 


Lochheed y, MISSILES AND SPACE DIVISION 


SUNNYVALE, PALO ALTO, VAN NUYS, 
SANTA CRUZ, SANTA MARIA, CALIFORNIA 
CAPE CANAVERAL, FLORIDA - ALAMOGORDO, NEW MEXICO 
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You can explore new areas at IBM im 


Electronic systems engineering covers the full spectrum of applied research and development at 
IBM. With this broad, flexible range, and IBM’s outstanding resources and facilities, your own 
imagination and creativeness will find unusual opportunity for expression. Currently, IBM 
engineers and scientists concerned with broad systems problems in business, science and govern- 
ment are working toward self-optimizing computers. These computers may some day program 
themselves and arrive at the one solution to a problem. Digital computers are being applied 
to simulate all the problems of a business firm operating under specified economic and statistical 
conditions. Progress is being made in advanced studies for radically different data systems for 
terrestrial and stellar navigational problems. For problems like these, and many others, IBM 
needs people who want to convert challenges into careers. 


A NEW WORLD OF OPPORTUNITY. Both technical and administrative engineering careers offer 
parallel advancement opportunities and. rewards at IBM. You will enjoy unusual professional 
freedom, comprehensive education programs, the assistance of specialists of diverse disciplines, 
and IBM’s wealth of systems know-how. Working independently or as a member of a small team, 
your individual contributions are quickly recognized and rewarded. This is a unique opportunity 
for a career with a company that has an outstanding growth record. 
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CAREERS AVAILABLE 
IN THESE AREAS... 
APPLIED MATHEMATICS &. 
STATISTICS 
CIRCUIT DESIGN & 
DEVELOPMENT 
COMPONENT ENGINEERING 
COMPUTER ANALYSIS 
CRYOGENICS 
FLIGHT TEST ANALYSIS 
HUMAN FACTORS 
INERTIAL GUIDANCE 
INFORMATION THEORY 
LOGICAL DESIGN 
OPERATIONS RESEARCH 
PROGRAMMING 
RADAR CIRCUITRY 
THEORETICAL PHYSICS 
TRANSISTOR DESIGN 


QUALIFICATIONS: 


B.S., M.S., or Ph.D. in 
Electrical or Mechanical 
Engineering, Physics, or 


Mathematics—and proven ability 


to assume a high degree of 
technical responsibility in 
your sphere of interest. 


Marcu 1959 


GENERAL DESCRIPTIONS OF 
SOME ASSIGNMENTS: 


INERTIAL GUIDANCE ENGINEER to assume broad 
project leadership in the planning and controlling 
of development projects. Analyze relationship of in- 
ertial equipment with bombing and navigation com- 
puter. Must have experience in servo-mechanisms, 
astro-compass or similar devices. 


COMPUTER ENGINEER to perform physical, mathe- 
matical analyses for solving complex control prob- 
lems by use of digital computers. Applications in 
missile systems and special-purpose computer 
systems such as DDA plus extensive experience in 
computer analyses. 


RADAR ENGINEER to analyze ultimate limits and 
present techniques and to develop new concepts 
providing topographical sensors for advanced air- 
borne and space systems. Design airborne radar 
pulse, microwave and deflection circuitry. Analyze 
doppler radar systems to determine theoretical 
accuracy and performance limitations. Experience 
with transistor circuits or radar test equipment is 
highly desirable. 


TRANSISTOR ENGINEER to design transistor am- 
plifiers, delay lines, transistor-tube conversion cir- 
cuits, and develop specifications for transistorized 
equipment. 


ENGINEERING PHYSICIST to assist in design and 
development of advanced solid state computers. 
Must have extensive background in electronic funda- 
mentals plus knowledge of solid state phenomena 
and mathematics. Must be capable of participating 
in logic development with minimum supervision. 


For details, write, outlining 

background and interests, to: Mr. R. E. Rodgers, 
Dept. 572 C 
IBM Corp. 
590 Madison Avenue 


New York 22, N. Y. 


INTERNATIONAL BUSINESS MACHINES CORPORATION 


low-up means regulates the controls to 
change the course setting and inactivate 
the follow-up. 


Directional system sensitive to rates of 
turn (2,866,934). H. W. Whitehead, 
Marblehead, Mass., assignor to General 
Electric Co. 

Inertial directional reference means 
maintaining a stable spatial orientation. 
A slaving control corrects directional infor- 
mation upon departure of the reference 
means from a predetermined course. The 
turn signal interrupts operation of the 
slaving control when the signal exceeds a 
preset amplitude. 


Automatic tracking circuits (2,866,966). 
E. L. C. White, Iver, England, assignor to 
Electrical & Musical Industries, Ltd. 

Position defining pulses generated by 
four electrical wave forms which re- 
peatedly pass through a predetermined 
reference value. Pulses are derived from 
voltage variations set up across a load 
impedance connected to a point of fixed 
potential. 


Non-propulsion attachment for rockets 


(2,866,413). C. R. Roquemore, China 
Lakes, Calif., assignor to the U. S. Navy. 


2,866,413 
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Thrust nullifier for a rocket having a 
launcher retainer groove near its rear end. 
Gas from the rocket issues through holes in 
the wall in such radial directions that na 
transverse forces are applied to the rocket, 


Launching control device (2,867,153), 
J. Hirsch, Pacific Palisades, Calif., as- 
signor to Aerophysics Development Corp, 


System for gradually releasing a rocket 
propelled missile on a launching frame. 
Spring force on fingers gripping an adapter 
in alignment with the CG of the missile is 
retarded for a period greater than the 
period during which initial transient thrust 
eccentricities exist. 


Boundary layer control (2,867,392). W.P. 
Lear, Pacific Palisades, Calif., assignor to 
Lear, Inc. 

Hollow airfoils divided into compart- 
ments on either side of the longitudinal 
axis, and predetermined areas of the airfoil 
surfaces having apertures between the 
compartments and the external atmos- 
phere. Air pumps are arranged to supply 
equal boundary layer capacity to both 
sides of the axis. 
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Book Reviews 


Behaviour of Metals at Elevated Tempera- 
tures, Institution of Metallurgists, Lon- 
don, 1957, 122 pp. 


Reviewed by J. P. FRANKEL 
University of California 


This volumie contains the four lectures 
that comprised the Institution of Metal- 
lurgists Refresher Course given at Llan- 
dudno in October 1956. The lecture titles 
were: Engineering Properties of Metals at 
High Temperatures, Effect of Tempera- 
tures up to 450C on Metals, Non-Ferrous 
High Temperature Materials, and High 
Temperature Steels. 

With the exception of some data on S. 
A. P., the mixture of aluminum and its 
oxide, the discussions concentrate on met- 
als and alloys, excluding metal ceramic 
mixtures. 

The first lecture is an excellent summary 
of the import of materials properties on 
engineering design procedures; the re- 
maining lectures deal with specific tem- 
perature ranges and materials properties. 
This volume is of particular interest to the 
airframe and engine designer, since it dis- 
cusses thoroughly the problem of material 
selection for short time high temperature 
applications. The interests of designers 
in the long time application fields are by no 
means neglected, however. 

In short, the book is commended to the 
designer and others who have not the op- 
portunities to follow closely developments 
in the rapidly changing materials field. 
The only drawback, from the point of view 
of the reader in the U.S., is the use of 
British alloy data for representation of the 
basic principles. In many cases, as the 
book itself points out, the characteristic 
behavior of particular alloys is quite sensi- 
tive to composition; accordingly, the 
American reader must use this book pri- 
marily as an overall guide rather than a 
source of data. The book may be con- 
sidered as an excellent pocket guide to the 
exotic field of High Temperature Met- 
allurgy whose inhabitants obey strange 
laws, but tend to be somewhat friendly 
when understood, and quite amenable to 
exploitation. 


The Exterior Ballistics of Rockets, by 
L. Davis Jr., J. W. Follin Jr. and L. 
Blitzer, Van Nostrand Co., Princeton, 
1958, v + 457 pp. $8.50. 


Reviewed by A. E. Funs 
Northwestern University 


This book is a modified version of a pre- 
viously classified U.S. Navy publication of 
the same title. Most ofthe book was writ- 
ten immediately following World War II 
and considers the exterior ballistics of solid 
fuel rockets without moving control sur- 
faces. The theoretical analysis evolved 
simultaneously with a wartime ex- 
perimental program; hence the book is 


Ali Bulent Cambel, Northwestern University, Associate Editor 


based on tested techniques. The methods 
of solving ballistic problems, which are pre- 
sented by this clearly written book, are 
still pertinent today. Problems associated 
with more complex rocket systems can be 
solved by an extension of the basic treat- 
ment given in this book. 

Before we consider in detail the material 
included in the book, let us examine some 
of its features. In order to adequately de- 
velop the mathematical theory it is nec- 
essary to introduce an enormous array of 
symbols. To help the reader, a list of sym- 
bols is provided, and the section where the 
symbol is first introduced is indicated. In 
many cases the formulas are derived in al- 
ternate ways. The physical significance of 
each mathematical result is discussed suffi- 
ciently to give the reader insight into its 
meaning. Typical values for the various 
parameters are given to help fix in mind 
the order of magnitude of the terms. Pre- 
cise derivations are given from which more 
useful approximate expressions are de- 
veloped. Solutions are usually shown both 
in graphical and analytical form. The 
book is divided into two parts (fin- and 
spin-stabilized rockets) each preceded by 
an introductory chapter. 

The force system acting on the rocket 
consists of gravity, jet and aerodynamie 
forces. Assuming these forces to be 
known, the theory necessary to predict the 
trajectory is given. The trajectory can be 
considered to be composed of different seg- 
ments, namely, launching, motion during 
burning and motion after burning. Each 
segment of the tracjectory is considered 
fully, including launching from stationary 
sites and moving aircraft. 

To be able to understand the book, the 
reader should be familiar with vector 
analysis, elementary differential equations 
and the basic concepts from analytical 
mechanics. Where other mathematical 
tools are needed, the necessary background 
is presented. 


Random Vibration, edited by Stephen H. 
Crandall, The Technology Press, Cam- 
bridge, 1958, 390 pp. 


Reviewed by Rosert L. SurHERLAND 
University of Wyoming 


This book, consisting of the notes for the 
M.I.T. special summer program on Ran- 
dom Vibration, covers a range of topics in- 
cluding theory of random processes, re- 
sponse to random vibration, instrumenta- 
tion for random vibration analysis, simula- 
tion equipment and associated problems, 
and recommended design procedure for 
equipment subjected to random loading. 
There are also discussions somewhat more 
general in subject matter rather than spe- 
cifically concerned with random vibration 
theory, such as in the chapters on structural 
damping, fatigue of metals and examples of 
stochastic processes of mechanical origin. 
Hence it might well be titled (for example) 
“Considerations in the Design of Equip- 


ment Subjected co Random Loading.” 
Nevertheless, all sections contribute to the 
development and understanding of the 
overall picture of random vibration and 
associated problems. 

As pointed out by several of the con- 
tributing authors, random vibration repr«- 
sents an involved and not completely un- 
derstood subject. However, certain a-- 
pects of the theory are not entirely new, 
since electrical engineers have already et.- 
countered the complexities of randoin 
phenomena in their studies of noise in siy- 
nal transmission. The current necessity 
for the undertstanding of random me- 
chanical vibration is largely a result of the 
high level of noise and vibration energy 
generated by powerful rocket and jet en- 
gines of recent development. Thus the 
inclusion of chapters devoted to the re- 
sponse of structures to jet noise and the e-- 
timation of sound-induced missile vibra- 
tions is both timely and significant. 

The fundamental approach (including 
necessary limitations) to the designation 
and analysis of random processes in general 
and random vibration in particular is pre- 
sented in chapters 2 and 4. The most con- 
veniently handled representation of ran- 
dom excitation is that of the Gaussian 
probability distribution; the resulting 
power spectrum and probability distribu- 
tion of the response is presented. The dis- 
cussion of instrumentation for random vi- 
bration analysis is both thoroughly and 
ably handled, and the final chapter con- 
sisting of four papers on the random vibra- 
tion design problem is an interesting, ap- 
propriate and complementary conclusion 
for this collection of notes. 

In common with the usual compilation 
of papers by a number of authors, a certain 
amount of discontinuity of development 
and non-uniformity of style is unavoidable. 
This is overshadowed, however, by the 
value of the material as a reference source 
of information for engineers and scientists 
dealing with design and analytical prob- 
lems involving random processes. A help- 
ful addition would have been the inclusion 
of a good subject index, which would add 
greatly to the book’s usefulness as a ref- 
erence text. 


Landing Gear Design, by H. G. Conway, 
Chapman and Hall Ltd., London, 1958, 
viii + 342 pp. Distributed in America 
by Macmillan Co., New York. 


Reviewed by James J. KAuZLARICH 
Worcester Polytechnic Institute 


This monograph on airplane landing 
gears is the first such treatment of this 
subject in book form known to the re- 
viewer. The author has attempted to give 
a full coverage emphasizing practical de- 
sign advice. The context is principally 
historical, treating recent landing gear de- 
signs in detail. This reviewer feels that 
the material is a light treatment of the sub- 
ject having little depth in discussing topics 
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Engineers — Scientists 
Required for many aspects 
of General Electric’s 
Aircraft Nuclear Propulsion 
Program... 


APPLIED 
RESEARCH & 
DEVELOPMENT 
AT THE PhD 
LEVEL 
ON AN INTER- 
DISCIPLINARY 
BASIS 


Physics or 
Nuclear Engineering 
Physical Chemistry 
or Ceramic 


Engineering 


Opportunities now exist for a 
number of scientists and engineers 
with imaginative minds and estab- 


lished reputations to contribute to Mathematics, 
extremely sophisticated nuclear Electrical 
Engineering 


propulsion systems, where theoret- 
ical advances must be wedded to 
new standards of reliability and 
performance. 


Qualified engineers and scien- 
tists are invited to inquire about 


these opportunities to pursue orig- Physics 
inal work in a pioneering area with or Nuclear 
General Electric. Engineering 
Write in confidence, including 0) 
salary requirements, to: Cy Physical Metallurgy 
Mr. G. W. Travers, Div. 36-MC or Mechanical 
| Engineering 
AIRCRAFT NUCLEAR PROPULSION DEPT. GB)ssssess: 
Nuclear 
GENERALQQELECTRIC | 


P.O. Box 132 Cincinnati 15, Ohio 


Advise Materials Testing Laboratory on materials prob- 
lems and techniques. Direct experimental, physical, me- 
chanical and electrical tests. Analyze data. 


Carry out reactor analysis and nuclear design of power 
plants, including formulation of nuclear design specifica- 
tions and test requirements. 


Supervise applied ceramic research including fabrication 
processes, testing of ceramic bodies; development of novel 
ceramic materials. 


Conduct theoretical investigations of effect of neutrons and 
photons on matter. 


Engineering analysis of physical systems in electro-me- 
chanical areas, deriving equations associated with systems 
study. Develop generalized digital program for parametric 
study of these systems. 


Apply advanced mathematical procedures and approaches 
in resolving complex problems in areas of reactor analysis. 


Advanced Reactor Analysis. Plan objectives of experimen- 
tal reactor investigations, which adequately combine po- 
tentialities of experimental techniques and improvement 
of theoretical methods. 


Undertake responsibility for studies regarding various ma- 
terials used in moderator materials development. 


Plan, design and analyze specific nuclear shield tests — 
part of broad nuclear shield design program. 


Participate in shield physics experiments, to provide fun- 
damental data for shield design. 


A number of additional positions in areas described above are also available to candidates with a Master’s Degree. 
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of current interest, e.g., wheel shimmy, 
skidding, stability, etc. As the title im- 
plies, the book covers many disciplines of 
engineering. In particular, the landing 
gear tires, wheels, brakes and shock ab- 
sorbers are treated in separate chapters. 
In addition, there are chapters on energy 
absorption, general layout of the landing 
gear, retraction and landing gear stressing. 
In each chapter the important require- 
ments for the particular element of the 
landing gear under discussion are summa- 
rized. This is usually followed with ex- 
ample calculations based on simplified 
analysis or empirical results. In cases 


Jet and Rocket Propulsion 
Engines 


The Hot Water Rocket, by Otto Miill- 
hauser, Raketentechnik und Raumfahrt- 
forschung, vol. 2, no. 2, April 1958, pp. 
38-44 (in German). 


Altitude-test-chamber Investigation of 
Performance of a 28-inch Ram-jet Engine, 
IV: Effect of Inlet-air Temperature, Com- 
bustion-chamber-inlet Mach Number, and 
Fuel Volatility on Combustion Perform- 
ance, by Robert W. Kahn, Shigeo Nakani- 
shi and James L. Harp Jr.. NACA Res. 
Mem. E51D11, July 1951, 27 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 127, p. 13, 6/5/58.) 


Altitude Test Chamber Investigation of 
Performance of a 28 inch Ram Jet Engine, 
III: Combustion and Operational Per- 
formance of Three Flame Holders with a 
Center Pilot Burner, by Thomas B. Shil- 
lito, George G. Younger, and James G. 
Henzel Jr.. NACA Res. Mem. E50J20, 
Feb. 1951, 30 pp. (Declassified from Con- 
fidential by authority of NACA Res. 
Abstr. 127, p. 18, 6/5/58.) 


A Proposed Ram-jet Control System 
Operated by Use of Diffuser Pressure 
Recovery, by Maxime A. Faget. NACA 
Res. Mem. L52E05D, Sept. 1952, 19 pp. 
(Declassified from Confidential by au- 
thority of NACA Res Abstr. 127, p. 18, 
6/5/58.) 

Vanguard Engine, by Louis Michelson, 
Missiles and Rockets, vol. 3, May 1958, 
pp. 139, 141-142. 

Some Actuator Disc Theories for the 
Flow of Air through an Axial Turbo Ma- 
chine, by J. H. Horlock, Gt. Brit., Aeron. 


Epitor’s Note; Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of the 
Heat Transfer Laboratory, University of 
Minnesota, are gratefully acknowledged. 


Technical Literature Digest 


where empirical data are used, the refer- 
ence is given, but without comment. 
There is no separate bibliography. 

The discussion of typical landing gear 
designs includes American, British and 
Germ n design practices. However, the 
author has, in general, ignored the exten- 
sive resear h activities concerning landing 
gear problems that have been conducted in 
America. The Wright Air Development 
Center and the Langley Aeronautical Lab- 
oratory, in addition to the very extensive 
facilities of private companies, have made 
many significant contributions on the sub- 
ject. For example, the author has over- 


looked the analytical treatment of landing 
gear wheel shimmy including extensive 
testing of the theory that has been success- 
fully conducted by the Wright Air De- 
velopment Center (see J. Aeron. Sci., vol. 
21, p. 793). 

In summary, “Landing Gear Design,’ 
is a book that may ke useful to the novice 
aeronautical engineer or student, but does 
not cover the subject in sufficient depth to 
be very useful for the experienced designer 
confronted with a critical design problem. 
The book does, however, serve the useful 
purpose of giving the reader a broad pe:- 
spective of the subject. 


M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


Res. Council. Rep. & Mem. 3030 (formerly 
ARC Tech. Rep. 15491), 1958, 31 pp. 

A Study of High Energy Level, Low 
Power Output Turbines, American Ma- 
chine and Foundry Co., Turbo Div., Rep. 
1196, 1 Feb. 31 Dec. 1957, 181 pp. 63 figs. 

A Study of Liquid Boric Oxide Particle 
Growth Rates in a Gas Stream from a 
Simulated Jet Engine Combustor, by 
Paul C. Setze, NACA Res. Mem. F55120a, 
April 1957, 41 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 126, p. 31, 5/2/58.) 

Thermopropulsive Characteristics of 
High-speed Thrust Generators, by A. F 
Charwat, Aero-Space Engineering (for- 
merly Aeron. Engng. Rev.), vol. 17, June 
1958, pp. 45-48. 


Hydrogen Peroxide Injection Boosts 
Jet Thrust 120%, by L. E. Varadi, Avia- 
tion Age, vol. 29, June 1958, pp. 54-58. 


Better Turbopump Design Ups Aircraft 
Range, by Joseph Bendersky, Aviation 
Age, vol. 29, June 1958, pp. 72-77. 

Power Sources for Space Explored, by 
James A. Fusca, Aviation Week, including 
Space Technology, vol. 68 June 16, 1958, 
pp. 234-241. 


Advanced Propulsion Research Speeded 
Aviation Week, including Space Technology, 
vol. 68, June 16, 1958, pp. 277-280. 

Propellant Vaporization as a Criterion 
for Rocket-engine Design: Experimental 
Effect of Fuel Temperature on Liquid- 
oxygen-Heptane Performance, by M. F. 
Heidmann, NACA Res. Mem. E57E03, 
July 1957, 21 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 126, p. 32, 5/2/58.) 

A Study of Injection Processes for 
Liquid Oxygen and Gaseous Hydrogen in 
a 200-pound-thrust Rocket Engine, by 
Carmon M. Auble, NACA Res. Mem. 
E56125a, Jan. 1957, 32 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstr. 126, p. 32, 5/2/58.) 

Injection Principles for Liquid Oxygen 
and Heptane Using Two-element Injec- 
tors, by Marcus F. Heidmann, NACA 


Res. Mem. £56D04, June 1956, 30 pp. 
(Declassified from Confidential by au- 
thority of NACA Res. Abstr. 126, p. 31, 
5/2/58.) 

Combustion Instability in Liquid Pro- 
pellant Rocket Motors, 23d Quarterly 
Progress Report for the Period 1 Nov. 
1957 to 31 Jan. 1958, Princeton Univ., 
Dept. Aeron. Engng., Rep. 216-w, March 
1958, 45 pp., 2 tabs., 22 figs. 

Investigation of Internal Film Cooling of 
1000 - pound - thrust Liquid -ammoni2- 
Liquid-oxygen Rocket-engine Combustion 
Chamber, by Gerald Morrell, NACA Res. 
Mem. E51K04, July 1951, 42 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 115, p. 11, 5/22/57.) 

Temperature Surveys of Rocket Motor 
Combustion Chambers, by Dwight I. 
Baker, Calif. Inst. Tech., Jet Prop. Lab., 
Progr. Rep. 24-5, Nov. 1952, 23 pp. 
(Declassified from Restrieted by authority 
of CIT-JPL Public. 38, p. 18, 1/3/55.) 

Small Power Plants for Use in Space, 
by Louis Rosenblum, Aero-Space Engngq., 
vol. 17, July 1958, pp. 30-33, 51. 

Booster Propulsion for Space Vehicles, 
by Richard Wentink, Aero-Space 
Engng., vol. 17, July 1958, pp. 40-45, 51. 


How Do Small Gas Turbines Scale? 
by Dennis P. Edkins. Aviation Age, vol. 
30, July 1958, pp. 38-39, 41-45. 

Preliminary Survey of Propulsion Using 
Chemical Energy Stored in the Upper 
Atmosphere, by Lionel V. Baldwin and 
Perry L. Blackshear. Appendix D: Heat 
Transfer and Friction Drag, by James F. 
Schmidt, NACA TN 4267, May 1958, 
73 pp. 

An Analysis of Ram-jet-engine Time 
Delay Following a Fuel-flow Disturbance, 
by Fred A. Wilcox and Arthur R. Ander- 
son, NACA Res. Mem. E55D22, June 
1955, 38 pp. (Declassified from Confiden- 
tial by authority of NACA Res. Abstr. 
126, pp. 31, 5/2/58.) 

Effect of Mechanically Induced Sinus- 
oidal Air-flow Oscillations on Operation 
of a Ram-jet Engine, by E. E. Dangle, A. 
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THE MOST ADVANCED 
SOLID PROPELLANT PROPULSION SYSTEMS 


Astrodyne has the capability and experience to design, develop, and manufacture 
complete solid propellant rocket propulsion systems, extruded and cast propellants, 
rocket motors and boosters, and gas generator charges for auxiliary power units. 
Astrodyne built this capability by bringing together in one productive group an 
elite corps of scientists, engineers, and technicians. Men with experience in research, 
design, and manufacture of superior solid propellants for JATO units, gas generator 
charges, and large boosters were teamed with others who had pioneered in the 
development of the largest rocket engines in use today. 


This experienced team now is at work on the development of large-scale cast 
propellant motors, the design of advanced lightweight rocket hardware, and the 
formulation and testing of high-energy propellants with optimum physical proper- 
ties for ballistic missile applications. 


Inquiries are welcomed on all phases of the solid propellant field— from preliminary 


design to quantity production. 
ASTRODYN E. INC. mccrecor, Texas 
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J. Cervenka and Eugene Perchonok, 
NACA Res. Mem. E54D01, June 1954, 
24 pp. (Declassified from Confidential 
by authority of NACA Res. Abstr. 126, p. 
29, 5/2/58.) 

Internal Flow and Burning Characteris- 
tics of 16 Inch Ram Jet Operating in a 
Free Jet at Mach Numbers 1.35 and 1.73, 
by Eugene Perchonok and John M. Far- 
ley, NACA Res. Mem. 51C16, May 1951, 
37 pp. (Declassified from Confidential 
by authority of NACA Res. Abstr. 126, 
p. 26, 5/2/58.) 

Summary of Scale-model Thrust- 
reverser Investigation, by John H. 
Povolny, Fred. W. Steffen and Jack G. 
McArdle, NACA Rep. 1314, 1957, 14 pp. 
(Supersedes TN 3664.) 

Effect of Fuel-orifice Diameter on 
Performance of Heptane-oxygen Rocket 
Engines, by Richard. J. Priem and Martin 
Hersch, NACA Res. Mem. E57126, Feb. 
1958, 20 pp. (Declassified from Confi- 
dential by authority of NACA Res. Abstr. 
126, p. 33, 5/2/58.) 

Injection Principles for Liquid Oxygen 
and Heptane Using Nine-element In- 
jectors in an 1800-pound-thrust Rocket 

ine, by Richard F. Neu, NACA Res. 
Mem. E57E13, July 1957, 27 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 126, p. 32, 5/2/58. ) 
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Wind Tunnel Tests at Mach Numbers 
1.5 and 2.0 on Models of a 2.75 Inch 
Folding Fin Rocket, by B. F. Jaeger, F. 
H. Knemeyer and G. M. Schroedter, 


NAVORD Rep. 1284, (NOT'S 354) Feb. 
1951, 18 pp. (Declassified from Confi- 
dential by authority of NOTS ltr 7506- 
BJS: . 75-53839 dtd 19 May 58, 
Ref. OPNAV INST 5510. 1A, Art. 0408 
dtd 2 Oct. 54.) 


The Supersonic Lift on Blunt and Sharp 
Trailing e Folding Fins in Combina- 
tion with a Long Body, by Marvin L. 
Luther and G. M. aarasiee, NAVORD 
Rep. 1317, (NOTS 391) June 1951, 12 pp. 
(Declassified from Confidential by au- 
thority of NOTS Itr 7506-BJS:nh. Reg. 
75-53839 dtd 19 May 58. Ref. OPNAV 
INST 5510.1A. Art. 0408 dtd 2 Oct. 54.) 

Cross Wind Firing of Rockets: An 
Experimental Survey, by William R. 
Haseltine, NAVORD Rep. 1181, (NOTS 
237), Sept. 1949, 16 pp. (Declassified 
from Confidential by authority of NOTS 
Itr 7506-BJS:nh. Reg. 75-53839 dtd 19 
May 58. Ref. OPNAV INST 5510.14. 
Art. 0408 dtd 2 Oct. 54.) 

The Employment of Aerodynamic 
Forces to Obtain Maximum Range of a 
Rocket Missile, by D. F. Lawden, Aeron. 
Quart. vol. 9, May 1958, pp. 97-109. 

The Shock Pattern of a Wing-body 
Combination, Far from the Flight” Path, 
by F. Walkden, Aeron. Quart., vol. 9, 
May 1958, pp. 164-194. 

Exploratory Investigation of Aerody- 
namic Effects of External Combustion of 
Aluminum Borohydride in Airstream Ad- 
jacent to Flat Plate in Mach 2.46 Tunnel, 
by Robert G. Dorsch, John S. Serafini and 
Edward A. Fletcher, NACA Res. Mem. 
E57E16, July 1957, 91 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstr. 126, p. 33, 5/2/58.) 

Exploratory Investigation of Static- 


and Base-pressure Increases Resulting 
from Combustion of Aluminum Borohy- 
dride Adjacent to Body of Revolution in 
Supersonic Wind Tunnel, by John S. 
Serafini, Robert G. Dorsch and Edward 
A. Fletcher, NACA Res. Mem. E57E15, 
Oct. 1957, 49 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 126, p. 32, 5/2/58.) 

Combustion of Aluminum Borohydride 
in a Supersonic Wind Tunnel, by Edward 
A. Fletcher, Robert G. Dorsch and Melvin 
Gerstein, NACA Res. Mem. E55F07, 
Aug. 1955, 12 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 126, pp. 31, 5/2/58.) 

Note on Ballistic Trajectories and Or- 
bits, by Richard F. Hughes, J. Aeron. 
Sci., vol. 25, May 1958, pp. 330-331 

On the Character of the Instability of 
the Laminar Boundary Layer Near the 
Nose of a Blunt Body, by A. M. Kuethe. 
J. Aeron. Sci., vol. 25, May 1958, pp. 338- 
339. 


A Program of Systemic Wind Tunnei 
Tests of Long, Fin Stabilized Rocket 
Models, by G. M. Schroedter, L. T 
Jagiello and B. F. Jaeger, NAVORD 
Rep. 1049, (NOT'S 169), Nov. 1948, 28 pp. 
(Declassified from Confidential by au- 
thority of NOTS Itr 7507-BJS:nh, Reg. 
75-53839 dtd 19 May 58, Ref. OPNAV 
INST 5510.1A. Art. 0408 dtd 2 Oct. 54. 


Experimental Studies of Forces, Pres- 
sure Distributions and Viscous Effects on 
Long Inclined Bodies of Revolution at 
Mach 2.96, by C. S. Brown, M. L. Luther 
and G. M. Schroedter, NAVORD Rep. 
1281, (NOTS 351), Feb. 1951, 60 pp. 
(Declassified from Confidential by aw- 
thority of NOTS Itr 7506-BJS:nh Reg. 
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Investigation of the Laminar Aerody- 
namic Heat-transfer Characteristics of a 
Hemisphere-cylinder in the Langley 11- 
inch Hypersonic Tunnel at a Mach Num- 
ber of 6.8, by Davis H. Crawford and Wil- 
liam D. McCauley, NACA Rep. 1323, 
1957, 21 pp. (Supersedes 7'N 3706.) 

Forces and Moments on Oscillating 
Slender Wing Body Combinations at 
Supersonic Speed, Part II: Applications 
and Comparison with Experiment, by 
Haim Kennet, Holt Ashley and Robert L. 
Stapleford, MIT, Fluid Dynamics Res. 
Group. Rep. 57-5 (AFOSR-TN-114; 
ASTIA AD 152022), Dec. 1957, 55 pp. 

Effect of Wing Camber and Twist at 
Mach Numbers from 1.4 to 2.1 on the Lift, 
Drag, and Longitudinal Stability of a 
Rocket-powered Model Having a 52.5° 
Sweptback Wing of Aspect Ratio 3 and 
Inline Tail Surfaces, by Warren Gillespie 
J:., NACA Res. Mem. L56C16, May 1956, 
20 pp. (Declassified from Confidential 
by authority of NACA Res. Abstr. 126, 
p. 35, 5/2/58.) 

Low Amplitude Damping-in-pitch Char- 
acteristics of Tailless Delta-wing-body 
Combinations at Mach Numbers from 
0.80 to 1.35 as Obtained with Rocket- 
pewered Models, by Charles T. D’Aiutolo, 
NACA Res. Mem. L54D29, June 1954, 
34 pp. (Declassified from Confidential 
by authority of NACA Res. Abstr. 126, 
p. 34, 5/2/58.) 

Transonic Flight Test of a Rocket-pow- 
ered Model to Determine Propulsive Jet 
Influence on the Configuration Drag, by 
Carlos A. deMoraes, NACA Res. Mem. 
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L54D27, June 1954 16 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstr. 126, p. 34, 5/2/58.) 

Experimental Effects of Propulsive Jets 
and Afterbody Configurations on the 
Zero-lift Drag of Bodies of Revolution at 
a Mach Number of 1.59, by Carlos A. 
deMoraes and Albin M. Nowitzky, NACA 
Res. Mem. L54C16, Apr. 1954, 32 pp. 
(Declassified from Confidential by au- 
thority of NACA Res. Abstr. 126, p. 34, 
5/2/58.) 

Aerodynamics of Slender Bodies at 
Mach Number of 3.12 and Reynolds Num- 
bers from 2 X 10° to 15 X 10°, II: Aero- 
dynamic Load Distributions of Series of 
Five Bodies Having Conical Noses and 
Cylindrical Afterbodies, by John R. Jack 
and Lawrence I. Gould, NACA Res. 
Mem. E52C10, May 1952, 28 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 126, p. 57, 5/2/58.) 

Low-speed Measurement of Static 
Stability and Damping Derivatives of a 60° 
Delta-wing Model for Angles of Attack of 
0° to 90°, by Donald E. Hewes, NACA 
Res. Mem. L54G22a, Sept. 1954, 30 pp. 
diagrs., photo., tab. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 127, p. 20, 6/5/58.) 

Wind-tunnel Investigation of the Use of 
Spoilers for Obtaining Static Longitudinal 
Stability of a Canard-missile Model in 
Reverse Flight, by Herman §. Fletcher, 
NACA Res. Mem. L54K05, June 1954, 15 
pp., diagrs., tab. (Declassified from Con- 
fidential by authority of NACA Res. 
Abstr. 127, p. 19, 6/5/58.) 

Rocket Trajectory Equations and Dif- 
ferential Corrections in Spherical Coor- 
dinates, by E. F. Dobies, Calif. Inst. of 


Techn., Jet Propulsion Lab., Progr. Rep. 
20-233, March 1958, 12 pp. 

The Interaction of Shock Waves and 
Turbulent Boundary Layers, by Andrew 
G. Hammitt, J. Aeron. Sci., vol. 25, June 
1958, pp. 345-356. 

Note on the Lift of Slender Nose Shapes 
According to Newtonian Theory, by J. D. 
Cole, J. Aeron. Sct., vol. 25, June 1958, 
pp. 399. 

The Wave Drag of Nonlifting Combina- 
tions of Thin Wings and ‘‘Non Slender’”’ 
Bodies, by L. M. Sheppard, Gt. Brit., 
Aeron. Res. Council. Rep. & Mem. 3-76 
(formerly ARC Tech. Rep. 19319; Royal 
Aircr. Estab., TN Aero 2496), 1958, 17 pp. 

Investigation of Factors Affecting the 
Control and Stability of High Speed Air- 
craft and Missiles. A Review of Work 
Conducted in This Field (9/1/52 to 2/28/ 
1958), by Phrixos Theodorides, Univ. of 
Maryland, Inst. for Fluid Dynamics and 
Appld. Math., Tech. Rep. BN-129 (AFOS- 
R-T R-58-52; ASTIA AD 154215), April 
1958, 12 pp. 

The Pressure Drag Due to Turbulent 
Separation on Bodies of Revolution with 
Varying Boundary Layer Thickness, by 
Erik G. M. Petersohn, Sweden, Aeron. 
Res. Inst., Rep. 75, 1957, 9 pp. 

Aerodynamic Characteristics at Super- 
sonic Speeds of a Series of Wing-body 
Combinations Having Cambered Wings 
with an Aspect Ratio of 3.5 and a Taper 
Ratio of 0.2. Effects of Sweep Angle and 
Thickness Ratio on the Aerodynamic 
Characteristics in Pitch at M-2.-1, by 
Ross B. Robinson, NACA Res. Mem. 
L52E09, July 1952, 27 pp., diagrs., pho- 
tos., 2 tabs. (Declassified from Confiden- 
tial by authority of NACA Res. Abstr. 
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Aerodynamic Characteristics at Super- 
sonic Speeds of a Series of Wing-body 
Combinations Having Cambered Wings 
with an Aspect Ratio of 3.5 and a Taper 
Ratio of 0.2. Effects of Sweep Angle 
and Thickness Ratio on the Static Lateral 
Stability Characteristics at M-1.60, by M. 
Leroy Spearman and John H. Hilton Jr., 
NACA Res. Mem. L51K15a, Jan. 1952, 
31 pp., diagrs., 2 tabs. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 127, p. 17, 6/5/58.) 

Effect of Camber on the Drag of a Body 
of Revolution, by Robert R. Dickey, 
mga: Res. Mem. A56E23, Oct. 1956, 8 

p., diagrs., tabs. (Declassified from 
Confidential by authoriy of NACA Res. 
Abstr. 127, p. 12, 6/5/58.) 

Forces and Moments on Pointed and 
Blunt-nosed Bodies of Revolution at Mach 
Numbers from 2.75 to 5.00, by David H. 
Dennis and Bernard E. Cunningham, 
NACA Res. Mem. A52E22, Aug. 1952, 47 

pp., diagrs., photos, 11 tabs. (Declassi- 

from Confidential by authority of 
NACA Res. Abstr. 127, p. 12, 6/5/58.) 

Exact Solution of the Neumann Prob- 
lem. Calculation of Non-circulatory 
Plane and Axially Symmetric Flows about 
or within Arbitrary Boundaries, by A. M. 
O. Smith and Jesse Pierce, Douglas Air- 
craft Co., Inc., Rep. ES 26988, April 1958, 
106 pp. 


A Numerical Method for Evaluation 
Wave Drag, by Maurice 8. Cahn and Wal- 
ter B. Olstad, NACA TN 4258, June 
1958, 13 pp., diagrs., tabs. 

Low-speed Measurements of Rolling 
and Yawing Stability Derivatives of a 60° 
Delta-wing Model, by Joseph L. Johnson 
Jr., NACA Res. Mem. L54G27, Dec. 1954, 
17 pp., diagrs., tab. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 127, p. 20, 6/5/58.) 

Conditions Governing the Initial Slope 
of a Sonic Line on a Blunt Body, by Andrew 
Hammitt and John J. Sullivan, Princeton 
Univ., Dept. Aeron. Engny., R 
(Wright Air Devel. Center, 
—" AD 151066), Feb. 1958, 6 pp., 3 

gs. 

An Experimental Investigation at a 
Mach Number of 2.01 of the Effects of 
Body Cross-section Shape on the Aero- 
dynamic Characteristics of Bodies and 
Wing-body Combinations, by Harry W. 
Carlson and John P. Gapeynski, NACA 
Res. Mem. L55E23, July 1955, 29 pp., 
diagrs., 2 tabs. (Declassified from Con- 
fidential by authority of NACA Res. 
Abstr. 127, p. 21, 6/5/58.) 

Investigation of the Effects of Body 
Indentation and of Wing-plan-form Modi- 
fication on the Longitudinal Chaacteristics 
of a 60° Swept-wing-body Combination at 
Mach Numbers of 1.41, 1.61 and 2.01, by 
John R. Sevier Jr., NACA Res. Mem. 
L55E17, July 1955, '37 pp., diagrs., tabs. 
(Declassified from’ Confidential by au- 
thority of NACA Res. Abstr. 127, p. 20, 
6/5/58.) 

General Variational Theory of the Flight 
Paths of Rocket Powered Aircraft, Mis- 
siles and Satellite Carriers, by Angelo 
Miele, Purdue Univ., School of Aeron. 
Engng., Rep A-58-2 (AFOSR-T N-58-246; 
eg AD 154148), Jan. 1958, 57 pp., 11 

gs. 

Aerodynamic Characteristics at Super- 
sonic Speeds of a Series of Wing-body 
Combinations Having Cambered Wings 
with an Aspect Ratio of 3.5 and a Taper 
Ratio of 0.2. Effect at M-2.-1 of Nacelle 
Shape and Position on the Aerodynamic 
Characteristics in Pitch of Two Wing- 
body Combinations with 47° Sweptback 
Wings, by Cornelius Driver, NACA Res. 
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Mem. L52F03, July 1952, 27 pp., diagrs., 
(Declassified from Con- 
fidential by authority of NACA Res. 
Abstr. 127, p. 18, 6/5/58.) 

Aerodynamic Characteristics at Super- 
sonic Speeds of a Series of Wing-body 
Combinations Having Cambered Wings 
with an Aspect Ratio of 3.5 and a Taper 
Ratio of 0.2. Effects of Sweep Angle and 
Thickness Ratio on the Static Lateral 
Stability Characteristics at M-2.01, by 
Clyde V. Hamilton, NACA Res. Mem. 
L52E23, Aug. 1952, 37 pp., diagrs., 2 
tabs. (Declassified from Confidential by 
ene of NACA Res. Abstr. 127, p. 18, 
6/5/58 

An Experiment on Flat Plate Turbulent 
Boundary Layer Flow—the Effect of Local 
Fluid Addition on Friction and Velocity 
Distribution, by H. Barrow, J. Roy. Aeron. 
Soc., vol. 62, Feb. 1958, pp. 135-138. 

Heat Transfer Through the Laminar 
Boundary Layer on a Circular Cylinder in 
Axial Incompressible Flow, by D. E. 
Bourne and D. R. Davies, Quart. J. Mech. 
and Appl. Math., Feb. 1958, pp. 52-66. 

Recovery Temperatures and Heat 
Transfer Near Two-dimensional Rough- 
ness Elements at Mach 3.1; by P. F. 
Brinich, NACA TN 4213, Feb, 1958, 20 pp. 

On the Diffusion of a Chemically Reac- 
tive Species in a Laminary Boundary 
Layer Flow, by P. L. Chambre and J. D. 
Young, Phys. of Fluids, vol. 1, Jan.-Feb. 
1958, pp. 48-56. 

Combustion in the Laminar Boundary 
Layer of Chemically Active Sublimating 
Surfaces, by M. R. Denison and D. A. 
Dooley, J. Aeron. Sci., vol. 25, April 
1958, pp. 271-272. 

A Method for Computing the Turbulent 
Boundary Layer When a Longitudinal 
Pressure Gradient is Present, by K. K. 
Fediaevskii and A. 8. Ginevskii. Sov. 
Phys.-Tech. Phys., Feb. 1957, pp. 275- 
291. 


Transient Temperature Distribution in 
Aircraft Structures, by I. Frank., J. Aeron. 
Sci., vol. 25, April 1958, pp. 265-267. 

On the Heat Transfer to Constant Prop- 
erty Laminar Boundary Layer with 
Power-Function Free-stream Velocity and 
Wall-temperature Distribution, by I. 
Imai, Quart. J. Appl. Math., April 1958, pp. 
33-45 

Some Effects of Bluntness on Bound- 
ary-layer Transition and Heat Transfer at 
Supersonic Speeds, by W. E. Mocckel, 
NACA Rep. 1312, 1957, 14 pp. 

An Extension of the Linearized Theory 
of Supersonic Flow Past Quasi-cylindrical 
Bodies, with to igi 
Interference, by R. C. Lock, J. Fluid 
Mech., vol. 4, May 1958, pp. 33-63. 
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The Three Dimensional Flow Field of an 
Axial Compressor Blade Row-Subsonic, 
Transonic and Supersonic, by James F lliot 
McCune, Cornell Univ., Grad. School of 
Aeron. Engng. (AFOSR-T N-58-72; AS- 
TIA AD 148116 and 148117), Feb. 1958, 
259 pp. 

Measurement of the Effect of an Axial 
Magnetic Field on the Reynolds Number 
of Transition in Mercury Flowing through a 
Glass Tube, by Michel Bader and William 
C. A. Carlson, NACA T'N 4274, May 1958, 
8 pp. 

A Performance Analysis of Methods for 
Handling Excess Inlet Flow at Supersonic 
Speeds, by Donald P. Hearth and James 
F. Connors, NACA TN 4270, May 1958, 
26 pp. 


Droplet Impingement and Ingestion by 
Supersonic Nose Inlet in Subsonic Tunnel! 
Conditions, by Thomas F. Golder, NACA 
TN 4268, May 1958, 56 pp. 


Summary of Experimental Heat-trans- 
fer Measurements in Turbulent Flow for a 
Mach Number Range from 0.87 to 5.05, by 
Maurice J. Brevoort and Barbara D. Ara- 
bian, NACA TN 4248, May 1958, 43 pp 

A Preliminary Investigation of Methods 
for Improving the Pressure-recovery Char- 
acteristics of Variable-geometry Super- 
sonic-subsonic Diffuser Systems, by Lo- 
well E. Hasel and Archibald R. Sinclair, 
NACA Res. Mem. L57H02, Oct. 1957, 57 
pp. (Declassified from Confidential b) 
authority of NACA Res. Abstr. 126, p 
37, 5/2/58.) 

Internal Performance Characteristics of 
Short Convergent-divergent Exhaust Noz- 
zles Designed by the Method of Charac- 
teristics, by H. George Krull and Willian: 
T. Beale, NACA Res. Mem. E56D27a 
July 1956, 17 pp. (Declassified from Con 
fidential by authority of NACA Res. 
Abstr. 126, p. 32, 5/2/58.) 


Matching of Auxiliary Inlets to Second- 
ary-air Requirements of Aircraft Ejector 
Exhaust Nozzles, by Donald P. Hearth 
Gerald W. Englert and Kenneth L. Kowal- 
ski, NACA Res. Mem. E55D21, Aug 
1955, 39 pp. (Declassified from Confi- 
dential by authority of NACA Res. Abstr 
126, p. 30, 5/2/58.) 

Heat Transfer from Surfaces on Non- 
uniform Temperature, by J. T. Stuart, 
J. Fluid Mech., vol. 4, May 1958, pp. 1-21. 


On the Theory of Vortex Flow, by K. 
H. Miller, ZAMM, vol. 38, May-June 
1958, pp. 179-186 (in German). 


On the Calculation of Radial Turbine 
Cascade Flow, by W. H. Isay, ZAMM, 
vol. 38, May-June 1958 pp. 209-219 (in 
German). 


Laminar Heat Transfer Around Blunt 
Bodies in Dissociated Air, by N. H. Kemp, 
P. H. Rose and R. W. Detra, AVCO Mfg. 
Corp., AVCO Res. Lab., Res. Rep. 15, May 
1958, 31 pp., 5 figs. 

Experimental Investigation of Heat 
Transfer by Forced Convection from the 
Hot Gas to a Cooled Stainless Steel Radial 
in Flow Gas Turbine Rotor, by Thomas R. 
Swartwout, Purdue Univ., Jet Prop. 
Center, Rep. F-58-2, Jan. 1958, 69 pp. 


On the Strength Distribution of Noise 
Sources Along a Jet, by H. S. Ribner, 
Univ. of Toronto, Inst. Aerophys., Rep. 51, 
April 1958, 19 pp., 7 figs. 

UTIA Air Duct Facility for Investigation 
of Vibration Noise Inducted by Turbulent 
Flow Past a Panel (Boundary Layer 
Noise), by L. Maestrello, Univ. of Toronto, 
Inst. Aerophys., TN 20, April 1958, 5 pp., 2 
tabs., 11 figs. 

Note on Acoustic Energy Flow in a 
Moving Medium, by H. 8. Ribner, Univ. 
of Toronto, Inst. Aerophys., TN 21, April 
1958, 6 pp., 2 figs. 

The Interaction of a Plane Strong Shock 
Wave with a Steady Magnetic Field, by 
J. H. deLeeuw, Univ. of Toronto, Inst. 
Aerophys., Rep. 49, March 1958, 47 pp., 
3 tabs., 29 figs. 


Laminar Flow in an Annulus with Porous 
Walls, by A. S. Berman, Atomic Energy 
Comm., K-1324, June 1957, 31 pp. 


Relaxation Phenomena in Air Between 
3000 and 8000° K, by M. Camac,J. Camm, 
J. Keck and C. Petty, AVCO Mfg. Corp., 
AVCO Res. Lab., Res. Rep. 22, March 
1958, 30 pp. 

Forced-convection Heat-transfer Char- 
acteristics of Molten Flinak Flowing in an 
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Inconel X System, by Milton D. Grele and 
Louis Gedeon, NACA Res. Mem. E53L18, 
Feb. 1954, 23 pp. (Declassified from 
Secret by authority of NACA Res. Abstr. 
127, p. 24, 6/5/58.) 

On the Thermal Conduction Equation 
for Self-heating Materials with Surface 
Cooling, by P. H. Thomas, TJ’rans. Fara- 
day Soc., vol. 54, Jan. 1958, pp. 60-65. 

Some Shock Tube Experiments on the 
Chemical Kinetics of Air at High Tempera- 
tures, by 8. Feldman, J. Fluid Mech., Dec. 
1957, pp. 225-242. 

Effect of Oxygen Recombination on 
One-dimensional Flow at High Mach 
Numbers, by S. P. Heims, NACA TN 
4144, Jan. 1958, 52 pp. 

Radiant Heat Exchange in a Gas-filled 
Enclosure, by H. C. Hottel and E. 8 
Cohen, AIChE J., vol. 4, no. 1, March 
1958, pp. 3-14. 

Dynamics of a Dissociating Gas, (III) 
Non-equilibrium Theory, by N. C. Free- 
man, NATO AGARD Rep. 133, July 
1957, 21 pp. 

On the Relation of Surface Friction and 
Heat Transfer, by E. A. Sidorov, Sov. 
Phys.-Tech. Phys., March 1957, pp. 499- 
504. 

An Analysis of the Turbulent Bound- 
ary-layer Characteristics on a Flat Plate 
with Distributed Light-gas Injection, by 
\M. W. Rubesin and C. C. Pappas, NACA 
7 N 4149, Feb. 1958, 43 pp. 

The Influence of Surface Injection on 
Heat Transfer and Skin Friction Associ- 
ated with the —o Turbulent 
Boundary Layer, by M. W. Rubesin 
NACA RM A55L13, Feb. 20, 1956, 16 pp. 

Heat Transfer in Isotropic Turbulence 
During the Final Period of Decay, by D. 
W. Dunn and W. H. Reid, NACA TN 
4186, June 1958, 68 pp. 

Prandtl Number Effects on Unsteady 
Forced-convection Heat Transfer, by E. 
M. Sparrow and J. L. Gregg, NACA TN 
4311, June 1958, 14 pp. 

Thermodynamic Properties and Com- 
position of Air by Dissociation Equilib- 
rium, by H. J. Kaeppeler and H. G., L. 
Krause, Forschugsinst. f. Physik der Strahl- 
antriebe, Mitteil. 1, 1954, 40 pp., 22 figs. 
(in German.) 

A Graphical Method for Solution of 
Freezing Problems, by P. 8S. Longwell, 
AIChE J., vol. 4, no. 1, March 1958, pp. 
53-57. 

Heat and Momentum Transfer in the 
Flow of Gases Through Packed Beds, by 
M. B. Glaser and G. Thodos, AIChE J., 
vol. 4, no. 1, March 1958, pp. 63-67. 


Heat Transfer from Horizontal Cylin- 
ders to a Turbulent Air Flow, by B. G. 
van der Hogge Zijnen, Appl. Sci. Res., 
Sec. A., no. 2-3, 1958, pp. 205-223. 

Heat Transmission to Fluids with Low 
Prandtl Numbers for Flow Through Tube 
Banks, by R. D. Cess and R. J. Grosh, 
ASME Trans., vol. 80, no. 3, April 1958, 
pp. 677-682. 


Heat Transfer to Fluids with Low Prandtl 
Numbers for Flow Across Plates and 
Cylinders of Various Cross Section, by R. 
J. Grosh and R. D. Cess, ASME Trans., 
vol. 80, no. 3, April 1958, pp. 667-671. 

Heat Transfer Rates to Cross-flowing 
Mercury in a Staggered Tube Bank (II), 
by C. L. Rickard, O. E. Dwyer and D. 
Dropkin, ASME Trans., vol. 80, no. 3, 
\pril 1958, pp. 646-651. 

Shock Wave Photogra ag of Large 
Subjects in Daylight, by EK. Edgerton, 
Rev. Sci. Instr., vol. 29, no. 2, Feb. 1958, 
pp. 171-172. 

Measurement of Metal Heat Conduc- 
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tivity in a Wide Range of Temperature 

a Single Experiment, by O. A. 
Kraev, Teploenergetika, vol. 4, no. 12, 
Dec. 1957, pp. 69-72 (in Russian). 

Unsteady Heat Convection Near a 
Horizontal Cylinder, by G. A. Ostroumov, 
Sov. Phys-Tech. Phys., Oct. 1957, pp. 
2627-2641. 

Natural Convection Heat Transfer in 
—— of Maximum Fluid Density, by 

8. Schechter and H. 8. Isbin, AIChE J., 
me 4, no. 1, March 1958, pp. 81-89. 

Heat Transfer in Nucleate (Il) 
New Correlating Equation, by K. Nishi- 
kawa, Trans. Japan Soc. Mech. "Engrs. " 
vol. 22, Aug. 1956, pp. 557-561. 

Heat Transfer in Nucleate Boiling (II) 
Effect of Pressure, by K. Nishikawa, 
Trans. Japan Soc. Mech. Engrs., vol. 22, 
Aug. 1956, pp. 562-570. 

On the Stability of Boiling Heat Trans- 
fer, by N. Zuber, ASME Trans., vol. 80, 
April 1958, pp. 711-718. 

Heat Transfer to Supercritical Water, 
by N. L. Dickinson and C. P. Welch, 
ASME Trans., vol. 80, April 1958, pp. 
746-752. 

Experimental Velocity and Temperature 
Profiles for Air in Turbulent Pipe Flow, by 
C. A. Sleicher Jr.,. ASME Trans., vol. 80, 
no. 3, April 1958, pp. 693-698. 

Turbulent Heat Transfer in the Thermal 
Entrance Region of a Pipe with Uniform 
Heat Flux, by E. M. Sparrow, T. M. Hall- 
man and R. Siegel, Appl. Sci. Res., sect. 
A, no. 1, 1958, pp. 37-52. 

Fluid Friction in Noncircular Ducts, by 
J. E. Walker, G. A. Whan and R. R. Roth- 
oo AIChE J., vol. 3, Dec. 1957, pp. 484- 
489. 
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Heat Transfer to Flow in a Round Tube 
with Arbitrary Velocity Distribution, b 
I. R. Whiteman and W. B. Drak, ASM. 
— vol. 80, no. 3, April 1958, pp. 728- 


Convective Heat Transfer from High- 
temperature Air Inside a Tube, by H. E. 
Zellnik and 8. W. Churchill, AJChE J., 
vol. 4, no. 1, March 1958, pp. 37-42. 

Quasilinear Heat Flow, by N. E. Fried- 
mann, ASME Trans., vol. 80, no. 3, 
April 1958, pp. 635-641. 

Graphical and Numerical Methods for 
Heat Conduction with Unsteady-state 
Ambient Temperature or Heat Genera- 


tion, by Y. Katto, 7’rans. Japan Soc. 
Mech. Engrs., vol. 22, Aug. 1956, pp. 
525-529. 


Formulae and Approximations for Aero- 
dynamic Heating Rates in High Speed 
Flight, by R. J. Monaghan, Gt. Brit., 
Aeron. Res. Council, Curr. Pap. 360, Oct. 
1955, 1957, 53 pp. 


On Heat Transfer in Laminar Boundary- 
layer Flows of Liquids Having a Very 
Small Prandtl Number, by G. W. Morgan, 
A. C. Pipkin and W. H. Warner, J. Aeron, 
Sci., vol. 25, March 1958, pp. 173-180. 


Incompressible Two-dimensional Stag- 
nation-point Flow of an Electricaly Con- 
ducting Viscous Fluid in the Presence of 
a Magnetic Field, by J. L. Neuringer and 
W. McIlroy, J. Aeron. Sci., vol. 25, March 
1958, pp. 194-198. 


Measurement of Turbulent Heat Trans- 
fer Rates on the Aft Portion and Blunt 
Base of a Hemisphere-cylinder in the 
Shock Tube, by J. Rabinowixz, Calif. 
Inst. Tech., Guggenheim Aeron. Lab., 
Mem. 41, Nov. 1957, 24 pp. 


Hypersonic Heat Transfer to Catalytic 
Surfaces, by S. M. Scala, J. Aeron. Sci., 
vol. 25, April 1958, pp. 273-275. 


Steady Temperature of a Plane Plate in 
a Reacting Gas Flow, by E. P. Vaulin, 
Sov. Phys.-Dokl., March-April 1957, pp. 
187-191. 


Film Boiling from Vertical Tubes, by 
Y. Y. Hsu and J. M. Westwater, AIChE 
J., vol. 4, March 1958, pp. 58-62. 

Effect of Superatmospheric Pressures on 
Nucleate Boiling of Organic Liquids, by 
R. B. Mesler and J. T. Banchero, AIChE 
J., vol. 4, March 1958, pp. 102-113. 

Heat Transfer in Nucleate Boiling (1) 
Effective Stirring Length of Bubbles, by 
K. Nishikawa and M. Kawano, Japan 
Soc. Mech. Engrs. Trans., vol. 22, Aug. 
1956, pp. 551-556. 

Far Noise Field of Air Jets and Jet 
Engines, by Edmund E. Callaghan and 
Willard D. Coles, NACA Rep. 1329, 
pp. (Supersedes 3590; 

591. 


Some Measurements Bearing on the 
Principle of Operation of Jet Silencing 
Devices, by G. M. Corcos, Douglas Aircr. 
Co., Rep. SM-23114, March 1958, 39 pp. 

The Compressible Flow through Cas- 
cade Actuator Discs, by J. H. Horlock, 
Aeron. Quart., vol. 9, May 1958, pp. 110- 
130. 


Turbulent Mixing of Two Streams of 
Different Gases, by L. G. Napolitano, 
L’ Aerotecnica, vol. 37, Dec. 1957, pp. 301- 
310 (in Italian). 

Application of the Integral Methods to 
Study of the Laminar Motion of a Non- 
reacting Binary Gas Mixture, by A. Pozzi, 
L’ Aerotecnica, vol. 37, Dec. 1957, pp. 311- 
320 (in Italian). 

Experimental Data on Turbulent Wall 
- by A. Sigalla, Aircr. Engng., vol. 30, 

ay 1958, pp. 131-134. 


The Contra-flow Gas-to-gas Heat Ex- 
changer, by J. S. Turton, Aircr. Engng., 
vol. 30, May 1958, pp. 135-141. 

The Effect of Shock Waves on the Isen- 
tropic Efficiency of Convergent-divergent 
Nozzles, by B. W. Martin and F. J. Bay- 
ley, Roy. Aeron. Soc. J., vol. 62, May 1958, 
pp. 377-382. 

On Heat Transfer in Slip Flow, by 
Stephen H. Maslen, a Aeron. Sci., vol. 25, 
June 1958, pp. 400-40 

On Hypersonic Flow 
with a Magnetic Field, by Nelson H. 
Kemp, J. Aeron. Sci., vol. 25, June 195s, 
pp. 405-407. 

Macroscopic Foundation of Plasma 
Dynamics, by Taro Kihara, Phys. Sov. 
Japan J., vol. 13, May 1958, pp. 473-48 

Externally Pinched Plasma, by Tihiro 
Ohkawa and Taro Kihara, Phys. Sov. 
Japan J., vol. 13, May 1958, pp. 482-48::. 

The Instability of a Pinched Fluid with 
a Longitudinal Magnetic Field, by MM. 
Kruskal and J. L. Tuck., Roy. Soc., Proc., 
vol. A245, June 3, 1958, pp. 222-237. 

A Radial Heat-flow Apparatus for Liquid 
Thermal Conductivity Determinations, by 
A. R. Challoner, H. A. Gundry and R. W. 
Powell, Roy. Soc. Proc., vol. A245, June 3, 
1958, pp. 259-267. 

On Some Fundamentals in Magneto- 
fluid-mechanics, by Eugene E. Covert, 
MIT, Naval Supersonic Lab., Tech. Rep. 
247, March 1958, 300 pp. 

Magneto-gasdynamics of Hypersonic 
Couette Flow, by Z. O. Bleviss, Douglas 
Aircr. Co., Rep. SM 23098, Feb. 1958, 56 
pp. 

Attenuation in a Shock Tube Due to 
Unsteady-boundary-layer Action, by Har- 
old Mirels, NACA Rep. 1333, 1957, 19 pp. 
(Supersedes 7'N 3278.) 


Measurements of the Effect of Surface 
Cooling on Boundary Layer Transition on 
a 15° Cone, PartI: Tests at M = 2 and3 
in an 8’’x 9’’ Wind Tunnel at RAE-Bed- 
ford, by A. C. Browning, J. F. W. Crane 
and R. J. Monaghan, Gt. Brit., Aeron. 
Res. Council, Curr. Paper 381 (formerly 
ARC Tech. Rep. 19836; Roy. Airer. 
Estab., TN Aero 2527). 


An Experimental Study of the Turbulent 
Boundary Layer on a Shock-tube Wall, by 
Paul B. Gooderum, NACA TN 4243, June 
1958, 63 pp., diagrs., photos., tab. 

Stall Propagation in a Cascade of Air- 
foils, by Anthony R. Kriebel, Harry S. 
Seidel and Richard G. Schwind, NACA 
TN 4134, June 1958, 122 pp. 


Cooling Characteristics of a Transpira- 
tion-cooled Afterburner with a Porous 
Wall of Brazed and Rolled Wire Cloth, by 
William K. Koffel, NACA Res. Mem. 
E54E25, Aug. 1954, 68 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstr. 127, p. 15, 6/5/58.) 


The Structure of a Hydromagnetic 
Shock in Steady Plane Motion, by G.S.5 
Ludford, Univ. of Maryland, Inst. for 
Fluid Dynam. and Appl. Math., TN BN- 
131 (AFOSR-TN 58-418; ASTIA AD 
158221), April 1958, 23 pp. 


Design and Operation of a Sodium to 
Lithium to Air Heat Transfer System, by 
A. R. Crocker, R. L. Potter, R. J. Spera, 
T. D. McLay and S. H. Esleeck, Atomic 
Energy Comm., APEX 327, Aug. 1957, 
195 pp. 

The Analytical Determination of the 
Local Heat Transfer Characteristics of 
Gases Flowing Turbulently in the Thermal 
Entrance Region of a Circular Duct with 
Large Temperature Differences between 
the Gas and the Duct Wall, by Helmut 
Wolf and M. J. Zucrow, Project Squid, 
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Tech. Rep. PUR-35-M (available only on 
microcards). 


Induced Air Flows in Fuel Sprays, by 
Hikmet Binark and W. E. Ranz, Project 
Squid, Tech. Rep. PSU-1-P (ASTIA AD 
156507), April 1958, 25 pp. (available 
only on microcards). 


Internal Characteristics and Perform- 
ance of Several Jet Deflectors at Primary- 
nozzle Pressure Ratios up to 3.0, by Jack 
G. MeArdle, NACA TN 4264, June 1958, 
107 pp. 

Measurements on the Velocity of Sound 
in Oxygen Gas under High Pressure, by 
A. van Itterbeek and J. Zink, Appl. Scz. 
Res., vol. A7, no. 5, 1958, pp. 375-385. 


On Simple Waves in Perfect and Dis- 
sociating Gases: Theory and Application, 
by Haim Kennet, MIT, Fluid Dynam. 
Ris. Group. Rep. 58.1 (AFOSR-TN-58- 
302), Feb. 1958, 49 pp. 

Measurements of the Velocity Distri- 
bution in a Plane Turbulent Jet of Air, by 
B. G. van der Hegge Zijnen, Appl. Sci. 
Rs., sect. A, vol. 7, no. 4, 1958, pp. 256- 
275. 

Measurements of the Distribution of 
Heat and Matter in a Plane Turbulent 
Jet of Air, by B. G. van der Hegge Zijnen, 
Appl. Sct. Res., sect. A, vol. 7, no. 4, 
1958, pp. 277-292. 


Measurements of Turbulence in a Plane 
Jet of Air by the Diffusion Method and by 
the Hot-wire Method, by B. G. van der 
Hegge Zijnen, Appl. Sci. Res., sect. A, vol. 
7, no. 4, 1958, pp. 293-313. 


Magnetohydrodynamics: Hope for 
Space, by J. S. Butz Jr., Aviation Week, 
vol. 68, May 12, 1958, pp. 48-50. 


Time Delay between High-speed Pellets 
and Associated Luminosity and Ionization, 
by P. E. Tucker, R. A. Davidson, and W. 
S. Patridge, J. Appl. Phys., vol. 29, May 
1958, pp. 868-870. 


Hydromagnetic Effects on Stagnation- 
point Heat Transfer, by Joseph L. Neur- 
inger and William Mellroy, J. Aeron. Sci. 
vol. 25, May 1958, pp. 332-334. 


Magnetohydrodynamic Analysis of Heat 
Transfer Near a Stagnation Point, by 
Vernon J. Rossow, J. Aeron. Sci., vol. 25, 
May 1958, pp. 334-335. 

Calculation of Shock Hugoniots and Re- 
lated Quantities for Nitrogen and Oxygen, 
by R. A. Alpher and H. D. Greyber, 
Physics of Fluids, vol. 1, March-April 
1958, pp. 160-161. 

Turbulence Measurement in Electri- 
cally Conducting Fluids, by Salamon Es- 
kinazi, Physics of Fluids, vol. 1, March- 
April 1958, pp. 161-162. 

Pumping and Drag Characteristics of an 
Aircraft Ejector at Subsonic and Super- 
sonic Speeds, by Gerald G. Gorton, 
NACA Res. Mem. F54D06, June 1954, 19 
pp. (Declassified from Confidential by au- 
thority of NACA Res. Abstr. 126, p. 29, 
5/2/58.) 

Investigation of Performance of Several 
Double-shroud Ejectors and Effect of 
Variable-area Exhaust Nozzle on Single 
Ejector Performance, by C. W. Ellis, D. 
P. Hollister and H. D. Wilsted, NACA 
Res. Mem. E52D25, July 1952, 25 pp. 
(Declassified from Confidential by au- 
thority of NACA Res. Abstr. 126, p. 27, 
5/2/58.) 

Internal Performance Characteristics of 
Variable-throat Plug- and a 
Convergent-divergent Nozzles, by 
George Krull, Fred W. Steffen and Carl C. 
Ciepluch, NACA Res. Mem. E53D09, 
June 1953, 34 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 126, p. 28, 5/2/58.) 


Marcu 1959 


MISSILE HARDWARE 


NEWBROOK 


Specializing in 
eMOTOR CASES 


Solid and Liquid 
Propellants 


JATO CASES 

NOZZLES 

e PLENUM CHAMBERS 
eBLAST TUBES 

eFUEL INJECTORS 


We have the 
“KNOW-HOW” 
We have developed new 
techniques, new methods, 
new processes that effect 
production economy so 
necessary to a successful 
missile program. 

Here at Newbrook you 
will find men with experi- 
ence gained from doing... 
a modern plant with up- 
to-date equipment... pre- 
cision inspection to meet 
your most exacting quality 
control requirements. 
And most important, New- 
brook specialization re- 
sults in strict reliability! Let 
us help you with your Mis- 
sile Hardware problems. 


Finest Welding Facilities 
Certified Welders 


MACHINE CORPORATION 


45 MECHANIC ST. PHONE 45 SILVER CREEK, N. Y. 
237 


Isen- 
1958, 
1958, 
Ex 
nd 3 Pe 
i 
CA Below: Machining Motor Cases 
the 
\ 
rmal 
yuid 
NAL 


AEROJET... 
for engineering careers 


STRESS ANALYSTS 
1. Experience and/or training in 
analytical stress computations, with 
emphasis on structural and pressure 
vessel applications. B.S. required; 
graduate work desirable. 


2. Experience and/or training in 
experimental stress analysis. Familiar- 
ity with advanced experimental tech- 
niques and expedient methods of stress 
analysis. B.S. required; graduate work 
desirable. 


GAS DYNAMICIST 

M.S. or Ph.D. in engineering, physics 
or applied mechanics. Five to ten years’ 
experience in analytical and experimen- 
tal work in the fields of boundary layer 
studies, supersonic and hypersonic 
flows, or application of non-ideal-gas 
theories. 


AERODYNAMICIST 

M.S. or Ph.D. in aeronautical or 
mechanical engineering or applied 
mechanics. Five to ten years’ experi- 
ence in the field of advanced aero- 
thermodynamics, missile aero- 
dynamics, heat transfer in rockets, 
missiles, etc. 


MISSILE SYSTEMS ANALYST 
B.S., M.S. or Ph.D. in engineering, 
physics or applied mechanics. Five to 
ten years’ experience in analysis and 
evaluation of interactions between pro- 
pellant systems and other missile com- 
ponents. Capable of developing 
optimum vehicle configurations from 
the standpoint of the propulsion sys- 
tem. 

U.S. Citizenship Required 

Resumes cordially invited. Write: 

E. P. JAMES 
AEROJET-GENERAL CORPORATION 
P.O. BOX 1947 
SACRAMENTO, CALIFORNIA 


Preliminary Experimental! Investigation 
of Transpiration Cooling for an After- 
burner with a Sintered, Porous Stainless- 
steel Combustion-chamber Wall, by Wil- 
liam K. Koffel, NACA Res. Mem. E53- 
DO8, June 1953, 47 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstr. 126, p. 28, 5/2/58.) 

Exploratory Tests of Transpiration 
Cooling of a Porous 8° Cone at M = 2.05 
Using Nitrogen Gas, Helium Gas, and 
Water as the Coolants, by Leo T. Chauvin 
and Howard S. Carter, NACA Res. Mem. 
L55C29, June 1955, 22 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstr. 126, p. 35, 5/2/58.) 

Theoretical Study of the Effect of the 
Non-ideality of a Dense Shocktube 
Driver Gas with Special Reference to 
Non-uniform Shocktubes, by A. E. Seigel, 
NAVORD Rep. 5707 (Naval Ord. Lab., 
Ball. Res. Rep. 5), Nov. 1957, 35 pp. 

Schlieren Evaluation of a Shock Tube 
for Pressure Gage Performance, by 
Stuart M. Hauser, Calif. Inst. Techn., Jet 
Prop. Lab., Rep. 20-115, Nov. 1957, 49 pp. 


Combustion, Fuels and 
Propellants 


Observations on the Kinetics of the 
Ionization of Alkali Metals in Flame 
Gases, by P. F. Knewstubb and T. M. 
Sugden, Faraday Soc. Trans., vol. 54, 
March 1958, pp. 372-380. 

The Thermal Reaction between Hydro- 
gen and Nitrogen Dioxide, PartIV: Ther- 
mal Ignitions of Hydrogen and Nitrogen 
Dioxide, by P. G. Ashmore and B. P. 
Levitt, Faraday Soc. Trans., vol. 54, March 
1958 pp. 390-393. 

The Dielectric Constant of a Compres- 
sed Gas Mixture, by A. D. Buckingham 
and R. E. Raab, Faraday Soc. Trans., vol. 
54, May 1958, pp. 623-628. 

Elementary Processes in the Radio- 
chemical Excitation of Gaseous Oxygen, by 
K. J. Laidler and E. Kerry Gill, Faraday 
Soc. Trans., vol. 54, pp. 633-639. 

Mixing and Chemical Reaction in the 
Laminar Wake of a Flat Plate, by S. I. 
Cheng and A. A. Kovitz, J. Fluid Mech., 
vol. 4, May 1958, pp. 64-80. 

Applications of Thermodynamic Con- 
stants and Data to Gaseous Mixing at 
Elevated Temperatures (Flames, Engines 
and Propulsion Systems), by G. Ribard 
and N. Manson, France, Ministere del’ Air, 
Publications Scientifiques et Techniques, 
no. 341, 1958, 194 pp. (in French). 

Performance Estimation of the Coolant 
Fuel System, by L. Greiner, J. Rabino- 
wicz and A. Zuckerman, Odin Assoc., Tech. 
Rep. 100-1 (AFOSR-T R-58-43; ASTIA 
AD 154150), Dec. 1957, 33 pp. 

Evaluation of Some Analytical Proce- 
dures for Solid Rocket Oxidants, Part I, 
by E. A. Burns and R. F. Muraca, Calif. 
Inst. Tech., Jet Prop. Lab., Mem. 20-158, 
June 1958, 17 pp. 

Ignition Limits of Hydrogen Peroxide 
Vapor at Pressures above Atmospheric, by 
C. N. Satterfield, F. Feakes and N. Sekler, 
Mass. Inst. Tech., Dept. Chem. Engngq., 
Rep. 53, June 1958, 9 pp. 

I. Macromolecular Symposium, Prague, 
Sept. 1957; II. Soviet Polymer Seience, 
1957, Condensed Report, by H. F. Mark, 
Solid Propellant Info. Agency, July 1958, 
50 pp. 

Relation between Gas Absorptivities and 
Emissivities, by D. Olfe and S. S. Penner, 
Calif. Inst. Tech., Daniel and Florence 
Guggenheim Jet Prop. Center, Tech. Rep. 
27, June 1958, 4 pp. 


High Energy Solid Fuels May Be 
Hydrids, by Michael Yaffee, Aviation 
Week, vol. 68, June 23, 1958, pp. 47, 50-51. 

Some Aspects of Detonation, Part I: 
Detonation Velocity and Chemical Con- 
stitution by A. R. Martin and H. J. Yallop; 
Part II: Detonation Velocity as a Func- 
tion of Oxygen Balance and Heat of For- 
mation, by A. R. Martin and H. J. Yallop, 
Faraday Soc. Trans., vol. 54, pt. 2, Feb. 
1958, pp. 257-271. 

New Method for Reducing Data on 
Interior Ballistic Trajectories, by Tadeusz 
Leser and John Lanahan, Aberdeen Pro»- 
ing Ground, Ballistic Res. Labs., Re). 
1024, Aug. 1957, 41 pp. 

Theoretical Performance of Liquid Am- 
monia with Liquid Oxygen as a Rocket 
Propellant, by Sanford Grodon and Al:in 
R. Glueck, NACA Res. Mem. E58A2\, 
May 1958, 84 pp. 

Theoretical Rocket Performance of 
Liquid Methane with Several Fluorine- 
Oxygen Mixtures Assuming Frozen Com- 
position, by Sanford Gordon and Micha»! 
E. Kastner, NACA Res. Mem. E58B2), 
May 1958, 44 pp. 

New Formulas for Rapid Calculation of 
Linear Burning Rates of Solid Propellants, 
by William F. Wallace, Picatinny Arsenal, 
Samuel Feltman Ammunition Labs,. Tec’. 
Rep. 2488, April 1958, 40 pp. 

Studies of OH, CO, CH, and C, Radia- 
tion from Laminar and Turbulent Pro- 

e-air and Ethylene-air Flames, by 
homas P. Clark, NACA TN 4266, June 
1958, 23 pp., diagrs., photos. 

Flame Structure Studies, IV: Premixed 
Carbon Monoxide Combustion, by Ray- 
mond Friedman and Robert G. Nugent, 
Project Squid, Tech. Rep. ARC-3-P (AS- 
TIA AD 156505), 18 pp., April 1958 
(available only on microcards). 

Considerations Pertaining to Spherical 
Vessel Combustion, by J. Grummer, E. B. 
Cook and T. A. Kubala, Project Squid, 
Tech. Rep. BUM-25-P (ASTIA AD 156- 
947), April 1958, 18 pp. (available only 
on microcards). 

Composition Traverses in Curved Lam- 
inar Fronts, by George H. Markstein, 
Project Squid, Tech. Rep. CAL-75-P 
(ASTIA AD 156503), April 1958, 21 pp. 
(available only on microcards). 


The Decay of Radical Radiation, by 
Edward T. Child and Kurt Wohl, Project 
Squid, Tech. Rep. DEL-9-P (ASTIA AD 
156506), April 1958, 23 pp. (available only 
on microcards). 

es Studies of Laminar 
Flames, II, by Raymond Kushida and 
Kurt Wohl, Project Squid, Tech. Rep. 
DEL-10-P (ASTIA AD 156504), April 
1958 (available only on microcards). 


Optical Studies of Structure of Turbu- 
lent Flames, by Kurt Wohl and Rajah H. 
Atalla, Project Squid, Tech. Rep. DEL-11- 
P (ASTIA AD 156948), April 1958, 20 
pp. (available only on microcards). 

Two Dimensional Effects in Gaseous 
Detonation Waves, by James A. Fay and 
Gerhard Opel, Project Squid, Tech. Rep. 
MIT-16-P (ASTIA AD 153326), March 
1958, 5 pp. (available only on microcards). 

Ionization in Detonation Waves, by 
Satyaki Basu and James A. Fay, Project 
Squid, Tech. Rep. MIT-17-P (ASTIA AD 
156502), April 1958, 16 pp. (available only 
on microcards). 

On the Structure of Detonation Waves, 
by T. C. Anderson, Project Squid. Tech. 
Rep. MICH-4-M (ASTIA AD 156509), 
April 1958, 19: pp. (available only on 
microcards). 
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Solid Propellant Burning, by Donald L. 
Reid, Kimball P. Hall and Martin Summer- 
field, Project Squid, Tech. Rep. PR-77-P 
(ASTIA AD 156501), April 1958, 19 pp. 
(available only on microcards). 


Investigation of Liquid Fluorine-Liquid 

Propellant Combination in a 

und-thrust Rocket Engine, by Ed- 

oe A. Rothenberg and Howard W. 

Douglass, NACA Res. Mem. E53E08, 

July 1953, 31 pp. (Declassified from 

Confidential by authority of NACA Res. 
Abstr. 127, p. 14, 6/5/58.) 

The Analysis of Red Fuming Nitric 
Acid, II: Determination of Water by the 
Karl Fischer Titration, by E. A. Burns 
and R. F. Muraca, Calif. Inst. Tech., Jet 
Prop. Lab., Progr. Rep. 20-344, Jan. 1958, 
8} 

A Statistical Theory of Chemical Reac- 
tion Rates, by James Keck, AVCO Mfg. 
Corp., AVO Res. Lab., Res. Rep. 20, April 
1958, 18 pp. 

The Reactions of Hydrogen Atoms with 
Oxygen, by B. deB. Darwent, and V. J. 
Kiasnansky, Project Squid, Tech. Rep 
CUA-1-P (ASTIA AD 157195), Apri 
1958, 18 pp. (available only on micro- 
cards). 

Diffusion and Heterogeneous Reaction, 
I: The Dynamics of Radical Reactions, by 
Henry Wise and Clarence M. Ablow, 
Project Squid, Tech. Rep. SRI-1-P (ASTIA 
AD 156508), March 1958, 17 pp. (avail- 
able only on microcards). 


Infrared Absorptivities and Emissivities 
of Gases, by A. Thomson, W. J. Hooker 
and S. S. Penner, Calif. Inst. Tech., Daniel 
and Florence Guggenheim Jet Prop. Center, 
Tech. Rep. 26, June 1958, 6 pp. 


Micro Determination of Nitrogen in 
Organic Compounds, by M. Bens, 
NAVORD Rep. 5865 *(WOTS 1988), 
April 1958, 26 pp. 


The Determination of Unsaturation in 
Glycols, by E. A. Burns, 

F. Muraca and F. Chang, Calif. Inst. 
Techn. Jet Prop. Lab., Progr. Rep. 20- 
345, Jan. 1958, 12 pp. 


Internal Ballistics of a High Velocity 
Light Gas Gun, by M. S. Sodha and V. K. 
Jain, Appl. Sci. Res., vol. A7, no. 5, 
1958, pp. 351-357. 


On the Internal Ballistics of a Supergun 
(Hochdrukpumpe), by V. A. Jain and M. 
$. Sodha, Appl. Sci. Res., vol. A7, no. 5, 
1958, pp. 369-374. 

Liquids vs. Solids, by Kurt R. Stehling, 
Aviation Age, vol. 29, June 1958, pp. 22- 
23, 90-91, 94-95. 


Advances Made in Conventional Fuels, 
by Michael Yaffee, Aviation Week, in- 
cluding Space Technology, vol. 58, June 
16, 1958, pp. 267-272. 


High Temperature Radiation Chemistry 
of Hydrocarbons, by P. J. Lucchesi, B. L. 
Tarmy, R. B. Long, D. L. Baeder and J. 
P. Longwell, Ind. Engng. Chem., vol. 50, 
June 1958, pp. 879-884. 


Gaseous Detonations, X: Study of 
Reaction Zones, by John P. Chesick and 
G. B. Kistiakowsky, J. Chem. Phys., vol. 
28, May 1958, pp. 956-961. 


Basic Considerations in the Combustion 
of Hydrocarbon Fuels with Air, Propul- 
sion Chemistry Division, edited by Henry 
C. Barnett and Robert R. Hibbard. 
Chap. I: Atomization and Evaporation of 
Liquid Fuels, by Charles C. Graves and 
Donald W. Bahr. Chap. II: Flow and 
Mixing Processes in Combustion Cham- 
bers, by Wilfred E. Scull and William R. 
Mickelsen. Chap. III: Ignition and Flam- 
mability of Hydrocarbon Fuels, by Frank 
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E. Belles and Clyde C. Swett. Chap. IV: 
Laminar Flame Propagation, by aden 
L. Dugger, a . Simon and Melvin 
Gerstein. Chap. V: Turbulent Flames, 
by Melvin Gerstein and Gordon L. Dug- 
er. Chap. VI: Flame Stabilization, by 

ordon L. Dugger and Melvin Gerstein. 
y- VII: Diffusion Flames, by Richard 
. Brokaw and Melvin Gerstein. Chap. 
VIII: Oscillations in Combustors, by 
Perry L. Blackshear Jr. and Warren D. 
Rayle. Chap. IX: Smoke and Coke For- 
mation in the Combustion of Hydrocar- 
bon-air Mixtures, by Rose L. Schalla and 
Robert R. Hibbard. NACA Rep. 1300, 
1957, 259 pp. 

Double Waves in Gaseous Detonations, 
by Mark Cher and G. B. Kistiakowsky, 
Harvard Univ., Dept. Chem., March 1958, 
15 pp., 3 figs. 

Chemically Frozen Boundary Layers 
with Surface Reaction, by D. E. Rosner, 
Princeton Univ., Dept. Aeron. Engngq., 
Rep. 419, March 1958, 17 pp. 

The Spark Ignition of Dust Clouds, by 
Lloyd E. Line Jr., Harold A. Rhodes and 
Thomas E. Gilmer Jr., Experiment Inc., 
TP-129, March 1958, 12 pp., 10 figs. 

Shock Tube Studies of Acetylene 
Compounds, by W. J. Hooker, Calif. Inst. 
Techn., Daniel and Florence ‘Guggenheim 
Jet Prop. Center, Tech. Rep. 7 (AFOSR- 
TN-58-293: ASTIA AD 154202), April 
1958, 10 pp., 5 figs. 

The Reaction Rate Law in Non-Uniform 
Flows, by D. E. Rosner, Princeton Univ., 
Dept. Aeron. Engng., Rep. 420, Nov. 
1957, 10 pp. 

On the Effect of Preliminary Vaporiza- 
tion on the Completeness and Stability of 
Combustion of a Fuel Spray, by V. Y. 
Basevich, Akademiia Nauk SSSR, Izves- 
tia, Otdelenie Tekhnicheskikh Nauk, no. 2, 
Feb. 1958, pp. 26-32 (in Russian). 

Effect of Pressure on the Normal Veloc- 
ity of the Propagation of Flames, by S. 
A. Gol’denberg and V. S. Pelevin, Aka- 
demiia Nauk SSSR, Izvestia, Otdelenie 
Tekhnicheskikh Nauk, no. 2, Feb. 1958, 
pp. 33-41 (in Russian). 

The Pressure - Volume - temperature 
Properties of Fluorobenzene, by D. R. 
Douslin, R. T. Moore, J. P. Dawson and 
Guy Waddington, Am. Chem. Soc. J., vol. 
80, May 5, 1958, pp. 2031-2038. 

Flow Field of Bunsen Flame, by M. S. 
Uberoi, A. M. Kuethe and H. R. Menkes, 
Physics of Fluids, vol. 1, March-April 
1958, pp. 150-158. 

NACA Probes Fluorine Uses, by Alfred 
J. Zaehringer, Missiles and Rockets, vol. 
3, May 1958, pp. 104-105. 

Stopping Time—the Photographv of 
Explosives, by Charles H. Bagley, SRI //., 
vol. 2, no. 1, 1958, pp. 27-32. 
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Outlook for 600 Deg. F. Silicone Rub- 
bers, by J. F. Dellaria, Aviation Age, vol. 
29, June 1958, pp. 60-64. 

Sheet Steels for High-speed Aircraft 
and Missiles, by A. L. Feild and M. E. 
Carruthers, Aero-Space Engng. (formerly 
Aeron. Engng. Rev.), vol. 17, June 1958, pp. 
41-44. 

Regularities in Creep and Hot Fatigue 
Data, Part II, by K. F. A. Walles and A. 
Graham, Gt. Brit., Aeron. Res. Council, 
Curr. Paper 380 (formerly Nat. Gas Tur- 
bine Estab. Rep. R190), Dec. 1956, 46 pp., 
97 figs. 

Elastomeric Mountings for High Tem- 
peratures, by R. P. Thorn, Aviation Age, 
vol. 30, July 1958, pp. 66-69. 
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SCIENTISTS 


Here is your opportunity to grow 
with a young, expanding subsidiary 
of the Ford Motor Company. Out- 
standing career opportunities are 
open in Aeronutronic’s new RE- 
SEARCH CENTER, overlooking the 
Pacific at Newport Beach, and the 
facility in Glendale, California. 
You will have all the advantages of 
a stimulating mental environment, 
working with advanced equipment 
in a new facility, located where you 
- enjoy California living at its 
inest. 


PhD and MS RESEARCH SPECIAL- 
ISTS with 5 to 7 years’ experience in heat 
transfer, fluid mechanics, thermodynamics, 
combustion and chemical kinetics, and 
thermoelasticity. To work on theoretical 
and experimental programs related to re- 
entry technology and advanced rocket pro- 
pulsion. Specific assignments are open 
in re-entry body design, high temperature 
materials studies, boundary layer heat 
t fer with ti thermal 
stress analysis, and high temperature 
thermodynamics. 

PROPULSION ENGINEERS with 5 years’ 
experience in liq id and solid rocket design 
and test. iliarity with _ heat transfer 
To work 
on pone col wide scope in R & D of ad- 
vanced concepts in rocket engine com- 
Fonents, and for missile project work. 


APPLIED MATHEMATICIANS, 3-5 
years’ recent experience required and A.B. 
or M.A. degree. Experience in numerical 
analysis and ter work in con 

with rockets and rocket propellants. 


STRESS HIGH TEMPERATURE MA- 
TERIALS. Mechanical or Metallurgical 
Engineer, must be familiar with thermal 
stress and shock as well as elasticity and 
plasticity. Application to nose cones and 
thrust chambers. 


FLIGHT TEST & INSTRUMENTATION 
ENGINEERS with 5 to 10 years’ experi- 
ence in laboratory and flight test instru- 
mentation techniques. Will develop tech- 
niques utilizing advanced instrumentation 
associated with space vehicles. 


THEORETICAL AEROTHERMODY- 
NAMICIST. Advanced degree and at 
least 5 years’ experience in high-speed 
aerodynamics. Knowledge of viscid and 
gas flows required. To work on program 
leading to advanced missile configurations. 
Work involves analysis of the re-entry of 
ypersonic missiles and space craft for 
determining optimum configuration. 


CERAMICIST. M.S. or Ph.D. required 
and 3-5 years’ recent experience with high 
temperature materials, structures and cer- 
mets. 


ENGINEER or PHYSICIST. With ex- 
perience in the use of scientific instruments 
for making physical measurement. Work 
related to flight test and facility instrumenta- 
tion. Advanced degree desired with mini- 
mum of 3 years of related experience. 


Qualitied a ates are invited to send 
and in to Mr. R. W. Speich, 
ic. 


AERONUTRONIC 


a subsidiary of Ford Motor Company 


2234 Air Way Bldg. 19, Glendale, Calif. 
CHapman 5-6651 


Newport Beach, Glendale, Santa Ana, 
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SELF-ALIGNING BEARINGS 


ROD 


PLAIN TYPES 


PATENTED U. S. A. 
World Rights Reserved 


CHARACTERISTICS 
ANALYSIS RECOMMENDED USE 


; For s operating under high 
Stainless Steel Ball and Race 1200 degrees F.). 
Chrome Alloy Steel Ball { For types operating under yt radial 
and Race ultimate loads ( -893,000 Ibs.). 


Bronze Race and Chrome { For types operating under normal loads 
Steel Ball with minimum friction requirements. 


DOO 


Thousands in use. Backed by years of service life. Wide variety of 
Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in similar 
size range with externally or internally threaded shanks. Our Engi- 
neers welcome an opportunity of studying individual requirements 
and prescribing a type or types which will serve under your demand- 
ing conditions. Southwest can design special types to fit individual 
specifications. As a result of thorough study of different operating 
conditions, various steel alloys have been used to meet specific 
needs. Write for Engineering Manual No. 551.Address Dept. JP-59. 


SOUTHWEST PRODUCTS CO. 


1705 SO. MOUNTAIN AVE., MONROVIA, CALIFORNIA 


RESEARCH ENGINEERS 


Stimulating and creative research positions are open 
at ARMOUR RESEARCH FOUNDATION in the fields of: 


Hypersonic Flight Problems 
Hydrodynamics 

Cavitation Problems 

Gas Dynamics 


Boundary layer investigations and wake flow prob- 
lems of high speed flight are some of the programs 
and concepts you will work on. 

M.S. to Ph.D. in Aerodynamics or Mechanical Engineer- 
ing required with good theoretical background 
compressible flow, heat transfer and thermodynamics. 
We offer outstanding employment benefits including 
liberal vacation policy. If you would like to work 
in a well known research organization with some of 
the leading engineers in this field, please write to: 


E. P. Bloch 


ARMOUR RESEARCH 
FOUNDATION 


Illinois Institute of Technology 
10 West 35th Street 
Chicago 16, Illinois 
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